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AND THEIR IMPACT ON A NOMINAL SHIP DESIGN 


by JAMES CLINTON DAVIS 

Submitted to the departments of Ocean Engineering and 
Electrical Engineering in partial fulfillment of the 
requirements of the degrees of NAVAL ENGINEER and MASTER OF 
SCIENCE in ELECTRICAL ENGINEERING AND COMPUTER SCIENCE. 

ABSTRACT 

Synchronous, permanent magnet, and induction machines 
are modeled using computer programs. The computer programs 
incorporate an optimization algorithm which converges on 
lowest weight, volume, and inefficiency. Machine designs 
for high and low rpms are performed, with a varying number 
of pole-pairs. The machine designs are analyzed to find the 
optimum combination of generator and motor for inclusion in 
a naval ship propulsion system. 

Three ships are used for the system study: a baseline 
mechanical transmission ship, a ship retaining the same sub¬ 
division as the baseline but with the electric machinery, 
and an electric transmission ship with subdivision and 
machinery box arrangement chosen to benefit from the in¬ 
herent arrangeability of electric transmissions. 

Two generator/motor combinations are used in the final 
ship analysis. Both employ a 3600 rpm, six-pole synchronous 
generator, which turns at the shaft speed of the prime 
mover. One combination uses a 180 rpm, direct-drive, 16- 
pole synchronous motor, and the other uses an 1800 rpm, 
geared, 8-pole synchronous motor. Power converters are used 
in both combinations to control motor speed. 

The geared combination in the rearranged ship 
demonstrated the best endurance speed efficiency, reducing 
the endurance fuel load by 18%, while maintaining the maxi¬ 
mum and sustained speed of the baseline ship. The savings 
in ship volume translated to an additional twenty Tomahawk 
missile cells in the rearranged ship. When the fuel load 
was held at the tonnage of the baseline ship, endurance 
range increased as much as 25%. 

Permanent magnet machines were not competitive in this 
study due to their high weight and volume, even though their 
individual machine efficiency was the highest of all types. 
Induction machines were not used as propulsion generators 
because of the inherent difficulties in control. The induc¬ 
tion machine motor candidates were not competitive because 
of off-design-point inefficiency. 
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Chapter One: Introduction 


The use of electric drive as a propulsion method for 
naval ships brings to the ship design process improved ar- 
rangeability and efficiency, though electric machines may 
increase the weight of the plant. Water-cooled electric 
machines are being studied today for ship transmissions 1 ; 
these are smaller and lighter, for the same power output, 
than air-cooled machines. They promise reduced overall 
weight for the ship through more economic prime mover load¬ 
ing, as less fuel will be needed on board. A review of the 
literature has found no work comparing various types of con¬ 
ventional motors and the effect of each type on the overall 
ship system when used as a propulsion method. St. John [1] 
showed the effects of a superconducting generator/motor 
transmission on the design of a DD963 destroyer hull. Many 
simulations of electric motors and their transients have 
been done. There have been several papers written on naval 
ship integrated electric propulsion systems 2 . Also, much 
effort has been expended in the area of electric motor 
design and optimization. Herein, various kinds of conven¬ 
tional electric machines are modeled. Those machines were 
used in ship designs to find the sensitivity of the designs 
to their use. 


1. Greene, Mole, Welch, and Seng, "Analysis of a High-Power 
Water-Cooled Electric Propulsion System," SNAME Trans., 

Vol 86, 1978, pp 140-162. 

2. Ames, "Marine AC Generation Systems," LSE Journal, 

Vol 12(1), pp 13-29. 


9 



1.1. Review of electric drive 


Navy ships operate in almost every salt water location 
in the world, including the Black Sea and the Indian Ocean. 
Regular deployments are made to the Mediterranean Sea, North 
and South Atlantic Ocean, and all areas of the Pacific. 

Ships transit the Suez and Panama Canals, the Saint Lawrence 
Seaway, and operate in the Arctic and Antarctic Oceans. 

These operations are made under greatly varying environmen¬ 
tal conditions, ranging from the sub-freezing temperatures 
of the high latitudes to the hot, dusty conditions of the 
Middle East. All this variety requires ships (and men) 
capable of sustained and efficient operation under any known 
condition. To that end, naval ships are tremendously com¬ 
plicated systems. 

It is impractical to equip naval ships for every con¬ 
tingency, but ship designers try to include as much as pos¬ 
sible when deciding what systems, equipment, and spares to 
put aboard a ship. Once a ship is designed, or as part of 
the design process, political acceptance of the product and 
its purchase is required. Since ships cost tax dollars, 
there are usually limitations on the size and complexity of 
the design. Still, the designers try to work within the 
given constraints and produce an acceptable and survivable 
(in both the battle and political senses) design. Usually 
this results in a ship that has very small margins in avail¬ 
able weight and volume. 3 

Weight (displacement) is a semi-direct measure of ship 
cost. Volume is required to place desired systems aboard a 
ship. Therefore, any design change that results in less 
required weight or volume with no decrease in ship 


3. To have small margins is to be limited in the quantity of 
additional system weight and volume that can be added to the 
ship over its lifetime. If a ship is limited in this 
fashion, the flexibility one has in backfitting new systems 
as the ship ages is significantly decreased. 
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effectiveness, is usually a welcome one. 

There are several different kinds of ship propulsion 
systems now in use. They include steam, nuclear, diesel, 
and gas turbine, and there are two principle drive systems: 
mechanical and electric drive. Below is a crude comparison 
of the various propulsion systems, for the purpose of plac¬ 
ing the thesis in perspective. 

Steam plants burning coal or oil have been in use for 
over a hundred years. They require large amounts of prime 
ship volume, in the center part of the ship. 4 Steam is 
produced in boilers and used to power turbines that rotate 
the shafts and propellers through reduction gears. Steam 
plants are very reliable mechanically, but are not terribly 
economic. The large size of the system components demands 
that the boilers and turbine machinery be placed in the cen¬ 
ter of the ship. This necessitates long runs of shafting 
from the center to the stern of a ship. Shafts typically 
are 18“-24" hollow steel cylinders of two to four inch 
thickness; they are heavy and their required placement and 
length makes valuable volume unavailable for other use. All 
propulsion plants except for electric drive have this ar¬ 
rangement and shafting disadvantage. Steam plants are used 
in all sizes of ships, from the 3000-ton displacement 
frigates to the 50,000-ton battleships to the 80,000-ton 
aircraft carriers. See Table 1 for a list of ship types and 
principle dimensions. 

Nuclear plants are steam plants with a different heat 
source. They also produce steam to power turbines and suf¬ 
fer the same volume disadvantages as conventional steam 
plants. They are also very heavy and very costly. Manning 
requirements are more stringent, since personnel levels are 


4. The center part of the ship is the most useful ship 
region to place and arrange systems. Ship designers try to 
keep free as much center ship volume as possible. This al¬ 
lows much greater flexibility in arranging systems that have 
large objects, such as boilers and turbines. 
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Table 1. Typical principle dimensions for various ship types 


Ship types Total 

shaft horsepower 

Length 

Tonnage 

Patrol Hydrofoil 

18,000 shp 

145 ft. 

240 

Frigate 

37,417 

441 

3880 

Destroyer 

78,555 

510 

6395 

Cruiser (non-nuclear) 

82,462 

556 

8872 

Cruiser (nuclear) 

77,460 

590 

9487 

Carrier (non-nuclear) 

280,000 

1050 

80,000 

Carrier (nuclear) 

270,000 

1100 

73,000 

Battleship 

212,000 

887 

58,000 

Submarine (nuclear, attack) 

15,000 

292 

4640 

Submarine (nuclear, strategic missiles) 

15,000 

410 

7880 


Source: Jane’s Fighting Shins 1985-86 . Jane’s Publishing 

Company, Ltd., London, edited by John Moore. 

These numbers represent the geometric mean of several 
classes of ships within a type and should not be taken to be 
those of a particular ship. Their value lies in the ap¬ 
preciation of the differences between various ship types. 

An apt weight comparison would be that a typical forty foot 
sailboat might displace fifteen tons. 

The U. S. Navy has other ship types besides those listed 
above. They include amphibious warfare ships, submarine and 
destroyer tenders, and fleet oilers, and supply ships. 
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rigidly controlled and the operators of the propulsion plant 
must be nuclear trained. Nuclear plants are used on 
cruisers, submarines, and aircraft carriers. 

Diesel propulsion plants have high weight to volume 
ratios, making their use costly in weight dollars, but 
cheaper in volume dollars. 5 They are very noisy, ruling out 
their use in antisubmarine warfare ships. Although very 
reliable mechanically, they require much tinkering and 
tuning. Their specific fuel consumption® is among the 
lowest of all the plants. Medium speed diesel engines are 
not commonly manufactured in the 25,000 hp range, which 
means diesel plants cannot be used in ships requiring high 
shaft horsepower. They are typically used in smaller ships 
as cruise engines and in amphibious warfare ships as main 
propulsion. (Amphibious warfare ships typically have lower 
top speeds than frigates or destroyers.) 

Gas turbine propulsion plants seem to have many good 
points. They are reliable, quiet, relatively low weight, 
and come in power ranges that are useful in ships ranging 
from 300 ton hydrofoils to 8000 ton destroyers. They 
require large amounts of volume for intake and exhaust duct¬ 
ing, but this is not a great disadvantage. Their fuel 
economy is not as good as other plants, but this is not 
intrinsic to the gas turbine engine. It is a fault of the 
operating method; gas turbine plants have mostly been built 
with mechanical transmissions. Usually one or two engines 
are coupled to each shaft. If the ship is proceeding at 
high speed, the gas turbines are operating at their full 


5. When designing a ship, total ship cost is monitored by 
the use of marginal cost factors. Every additional cubic 
foot or ton of weight added has a marginal cost associated 
with it. When a ship's total cost is constrained, design 
changes that add cost are discouraged or must be offset by 
the reduction of other systems’ weight or volume. 

6. Specific fuel consumption is the ratio of pounds of fuel 

burned per horsepower-hour. SFC=lbf/hp-hr. 
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load design point and are relatively fuel economic. Good 
fuel economy is not usually realized, however. In the main, 
the ship proceeds at a cruise speed, using one gas turbine 
for each shaft, and the gas turbines operate at about half 
power. Specific fuel consumption rises rapidly as gas tur¬ 
bine power level drops, which makes for inefficient opera¬ 
tion. 

Usually, the above propulsion plants have mechanical 
transmissions. This means the main engines, whether they 
are diesels, gas turbines or steam turbines, are mechani¬ 
cally connected to the shafts and propellers. There is 
usually a reduction gear between the engine and shafting. 
These gears are large, very heavy, and expensive. To 
provide an idea of size, the largest, or "bull” gear in a 
typical locked-train double-reduction gear is about seven 
feet in diameter. The reduction gear must be placed in-line 
with the shafting, thereby using more of that prime ship 
volume. Some mechanical transmissions have cross- 
connections between shafts, but this is not common. 

Electrical transmissions are characterized by prime 
movers of any type providing power to generators. The out¬ 
put electricity is conditioned and sent to propulsion motors 
via a distribution network. Cross-connection is done with 
switches and breakers. There can be a mechanical reduction 
gear if it is desired to operate the propulsion motors at 
higher than propeller rotational speeds. Direct-drive 
motors may also provide the desired propeller rpm, e.g., by 
controlling the field current in a DC motor. The propulsion 
motors can be very near the propeller, i.e. aft, eliminating 
the long runs of shafting associated with a mechanical 
transmission. 

Naval ship propellers are of two types, controllable or 
fixed pitch. The pitch of a propeller is the distance the 
ship moves forward in the water for one turn of the propel¬ 
ler. A fixed pitch propeller has this characteristic dis¬ 
tance the same at all times. A controllable pitch propeller 
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can vary this distance by changing the angle of attack (the 
angle at which the blade slices through the water) of the 
propeller blades, including reversing the blade so that the 
ship moves astern. Controllable pitch propellers are prac¬ 
tically required for propulsion plants that have non¬ 
reversing main engines, such as diesel and gas turbine 
plants. ^ Steam plants can reverse their propellers and 
shafts by use of an astern turbine, albeit much more slowly 
than a ship with a controllable pitch propeller system. 

Fixed pitch propellers have a slight advantage in efficiency 
(1-3%) over the controllable ones. This is mostly due to 
the large propeller hub required for varying the blade angle 
of a controllable pitch propeller. Quick reversal of shaft 
direction or propeller pitch means quick ship braking and/or 
ship reversal. This ship quickness is mandatory for an¬ 
tisubmarine operations and safe navigation. For example, 
the ability to stop "on a dime" may be important in a 
crowded sea lane, where a small wooden sailboat has the 
right of way over a powered naval vessel. 

Electric drive seems to combine the best of all the 
propulsion plants. It has all the advantages of a conven¬ 
tional prime mover plus the advantage of electrical cross 
connection and better arrangements. In the cruise scenario 
above, the electric drive ship could have both shafts 
operating from one gas turbine engine. That engine would be 
coupled to an electric generator which would produce enough 
power to run the motors that turn each shaft. The drive 
motors would be placed at the stern of the ship, near the 
propellers, on the same level. The long runs of shafting 
would be replaced by electric cable, which is smaller, 
weighs less, and can be placed in non-prime real estate. 


7. Ships with non-reversing prime movers can also have a re¬ 
versible reduction gear with a fixed pitch propeller instead 
of a controllable reversible pitch propeller. This is new 
technology for the United States and only the latest naval 
ship design, the DDG51, has a reversible reduction gear. 
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Cable is also in many cases cheaper than shafting, espe¬ 
cially to repair. Since one gas turbine would provide power 
to both shafts, it would operate at a higher power level and 
would be therefore more fuel economic. A typical propulsion 
plant might consist of three gas turbines with three gener¬ 
ators. Most ships have an even number of prime movers be¬ 
cause mechanical shaft cross connection between shafts is 
not often used and each shaft in a mechanical transmission 
ship requires the same number of prime movers to balance 
loading at high power levels. The extra prime mover 
requires a lot of weight and volume. An advantage of 
electric drive is that it becomes possible and perhaps 
desirable to use an odd number of prime movers. Each of the 
two shafts would have one propulsion motor. The heavy 
reduction gears could be replaced by the motors, which would 
have an infinitely variable reduction ratio. The control¬ 
lable pitch propeller system so far required by this gas 
turbine ship would be replaced by the cheaper, slightly 
smaller, and far less complicated fixed pitch propeller, ad¬ 
ding a small efficiency gain. The hydraulic system used to 
vary blade angle would be eliminated. Ship braking and 
reversal would be accomplished by electrically controlling 
the motor rotation direction, combined with energy dissipa¬ 
tion through the use of resistor banks. 

A disadvantage of this arrangement would be the high 
weight of the propulsion motors. They would be special 
designs and have a high capital cost. Hopefully, the high 
weight of the motors would be offset by the reduction in 
shafting and fuel weight and the possible elimination of the 
reduction gears. The high cost would be made palatable by 
the savings in fuel over the life of the ship. A Life Cycle 
Cost comparison of various propulsion plants, including 
electric drive, is available in reference two. 

The change to electric drive would likely be accom¬ 
panied by an overall decrease in propulsion plant weight and 
volume. The ship could be smaller and lighter, and would 
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require less onboard fuel for the same endurance range (the 
distance the ship can travel without refueling). Since less 
volume would be required for the fuel, the ship could be 
smaller and lighter. Since the ship would then be smaller 
and lighter, less horsepower would be required to achieve 
the same top speed. Since less horsepower would be 
required, the ship could be smaller and lighter. This is an 
example of the design spiral that would result in a smaller, 
lighter, cheaper, more risk-free ship. An example of this 
type of ship improvement is given in reference one. There 
is a limit on ship improvement, usually due to the non¬ 
propulsion systems or payload. One cannot make an ocean 
crossing missile ship the size of a small yacht. 

So why are not all Navy ships electric drive? They are 
not largely because the technology has not existed in a 
usable, fully developed, and manageable form. Because of 
the high cost of naval ships (a small one may cost $350 
million) and the lack of experience with current electric 
drive technology, the Defense Department is reluctant to 
build large electric drive ships. There have been electric 
drive ships, including five battleships with 21 MW shaft 
output and two aircraft carriers with 135 MW shaft output. 
Over 160 escort vessels were built during the Second World 
War with turboelectric or diesel-electric drives ranging 
from about 4.5 to 9 MW. a A new class of ocean surveillance 
ship, the T-AGOS 19, is being built with diesel-electric 
propulsion, but it is only a 3500 shaft horsepower (shp) 
ship. 

Electric drive was replaced by conventional mechanical 
transmission plants after World War Two because of the com¬ 
petition afforded by improved gear cutting methods. Double¬ 
reduction locked-train gear transmissions became the 
standard. Since the electric drive ships all had non- 


8. Doyle, T. J. and Harrison, J. H., "Navy Superconductive 
Machinery Program," Trans. SNAME, 1978, p. 20-1. 


17 
















































superconducting, air-cooled motors and generators, they had 
higher weight and increased space requirements and suffered 
in comparison with the mechanical transmission ships.9 The 
importance of the improvement in power electronics must be 
mentioned. World War Two ships did not have the advantages 
afforded by those electronics. 

Integrated electric drive propulsion must also be men¬ 
tioned. This propulsion plant is the same as any of those 
discussed above, except that ship service electrical power 
is derived from the main propulsion plant, usually by taking 
power off the reduction gears or main propulsion generator. 
Power conditioning equipment, such as a cycloconverter, is 
needed to "clean up" the power and change it to fixed 
frequency for use in other equipment. Variable speed con¬ 
stant frequency equipment and concepts embody the integrated 
electric drive concept. The U. S. Navy has investigated 
this in some detail .10 

Some requirements of electric drive may be viewed as 
disadvantages. The power from the electric generators has 
to be conditioned to provide frequency control of the 
propulsion motors. The power conditioners add weight and 
volume to the overall system, as well as reducing the 
efficiency of the transmission. Braking resistors, used to 
dynamically and quickly slow the propulsion motors, add more 
weight and volume to the system. There also may be a high- 
frequency radiated noise signature associated with alternat¬ 
ing current systems that may be deleterious to the mission 
of the ship. 

The research done to date has not explored specific 
motor types in detail. How can it be decided whether to put 
a synchronous, inductive, permanent magnet or other motor in 


9. Ibid. 

10. Robey, Stevens, and Page, "Application of Variable Speed 
Constant Frequency Generators to Propulsion-Derived Ship 
Service," Naval Engineers Journal, May 1985. 
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the electric drive system? What makes the ship system 
"best"? The effect of each motor type on the ship system 
has not been analyzed. The electric transmissions used in 
the current version of the Advanced Surface Ship Evaluation 
Tool (ASSET), a ship design computer program written for the 
United States Navy by Boeing Computer Systems, Inc., are 
generic combinations of AC and DC motors and generators, 
using rough estimates of weight, volume, and efficiency. 

Ship designs more involved than the feasibility level need 
detail on just those items. 

Motor design is a well known subject and there are many 
texts on the subject. The use of motors and pertinent tech¬ 
nologies in a ship as a propulsion method is discussed in 
Greene, Powell, and Gripp [3]. The advantages of electric 
drive include flexibility of arrangement, controllability, 
variable reduction ratio, reliability, and provision of ship 
service power from the main bus. Jolliff and Greene [4] go 
on to propose a specific water-cooled Advanced Integrated 
Electric Propulsion Plant (AIEPP) for a frigate/destroyer¬ 
sized ship. They discuss the essential characteristics of 
such a plant, establish the feasibility of the drive system 
and identify the method to technically demonstrate the sys¬ 
tem. Acker, Greene, and Jolliff [5] present several model¬ 
ing techniques and scaling relationships that allow estima¬ 
tion of volumes and weights of propulsion motors and gener¬ 
ators, solid state power conditioners, electrical 
switchgear, and associated electrical propulsion systems 
components as functions of propeller shaft power. A case 
study of AIEPP is given in the paper by Kastner, Davidson, 
and Hills [6]. 

Simulation of electric motors and associated systems is 
a popular topic. Many persons have done work in this area, 
from the micro-consideration of high frequency inductance 
changes to the more macro-consideration of hunting tran¬ 
sients, etc. Smith, Stronach, and Tsao [7] model a complete 
electromechanical marine drive system while Smith, Stronach, 
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Tsao, and Goodman [8] concern themselves more with a marine 
power system, including pump drives. Nonlinearities and 
operational transients are addressed. 

1.2. Optimization 

Optimization is the process of making a system, sub¬ 
system, or idea the best it can be. "Best" is defined by an 
"objective function," a measure of what is optimum. For ex¬ 
ample, an optimal manufacturing process may produce the max¬ 
imum number of units at the lowest cost. The objective 
function would combine units-produced with cost in an equa¬ 
tion that could be analyzed to find the proper production 
level. The output of the objective function is a single 
scalar measure of "goodness." It may be difficult to repre¬ 
sent complicated processes with only one number. 

Optimization can be performed on a global or subor¬ 
dinate basis. The optimum motor might be the one that has 
the highest efficiency, even if that efficiency was achieved 
by designing a very large and heavy motor. The sub-system 
(the motor) has efficiency as its objective function. The 
ship in which the motor is to be placed may be optimum when 
its overall weight and volume are the lowest (ignoring cost, 
for example). A large, heavy motor, then, may not be op¬ 
timum for the ship, even if it is very efficient. A good 
case study of motor optimization is the EPRI report authored 
by Fuchs, et al. [9] 

"Optimization" can be an ill-defined term but there are 
fairly well defined methods of achieving it. Linear 
programming, Markov modeling, and Monte Carlo schemes are 
examples of these methods. The accessibility of high-speed 
digital computers has made multiple random excursions in a 
multi-variable space a much easier way of finding the 
"optimum" solution, provided an objective function and con¬ 
straints can be devised to describe the problem. This 
method of random excursions (Monte Carlo scheme) and ex- 
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amples of it are among the methods described in references 
ten, eleven, and twelve. 

Monte Carlo schemes take their name from the action of 
the roulette wheel in the gambling casinos of Monte Carlo. 
Around and around the wheel goes, stopping on random num¬ 
bers. If the computing power is available, this is an ac¬ 
ceptable method of exploring a large variable space. It can 
be much quicker than looking at every possible permutation 
of all variables. 

The steepest-descent scheme is so named because op¬ 
timization moves down the sharpest gradient of the objective 
function. From a valid design point, random steps are taken 
in every variable and the design point is moved "downhill" 
toward the objective function over the steepest slope. This 
is different from the "drunkard’s walk," where the random 
steps are only evaluated on whether or not the objective 
function’s output has improved, not if it is improving at 
the fastest possible rate. 

Optimization is almost always subject to constraints. 

In the previous manufacturing process, warehouse space may 
be limited, so only a certain number of units may be stored. 
This could act to limit production. For motor design, con¬ 
straints include maximum rotor tip speed, maximum current 
density, minimum power output, etc. All constraints should 
be combined with the objective function to yield a 
"constrained optimization." 

1.3. The objective function 

Optimization is not possible without an objective func¬ 
tion. It may be very difficult to devise a good objective 
function for a complicated system such as a ship. It may be 
even harder to find one for a sub-system of that complicated 
system. There are very many characteristics that could be 
optimizing variables, and assembling them into one objective 
function with all the constraints is not easy. 
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Even with a properly defined objective function, it may 
be difficult to choose among designs that result from the 
constrained optimization. For example, a low-volume and 
low-weight motor may have poor efficiency. A very efficient 
motor may also be large and heavy. If all three elements 
are important, which is the best motor? Deciding between 
competing designs has been the subject of various papers, 
one of which is by Schweppe and Merrill [13]. In that 
paper, the authors suggest the use of "knee curves," saying 
that the essential characteristics of a multiple attribute 
tradeoff can be plotted on a series of x-y graphs. Uncom¬ 
petitive designs are easily discerned and discarded. The 
decision process can be limited to only those designs that 
are competitive. 

Table 2. Optimizing characteristics for ship and transmis¬ 
sion 


Weight 

Volume 

Efficiency 

Cost 

Reliability 

Maintainability 

Commonality 

Manning 


The above table lists many of the possible optimizing 
characteristics for the ship and its propulsion sub-system. 
Manning estimates are typically based on historical data and 
do not indicate that the baseline and variant ships will 
require a significantly different number of personnel. Com¬ 
monality measures the use of the same equipment in other 
ships. Since there are no other electric drive ships at the 
power levels used in this thesis, commonality is not an 
issue. It is very difficult to discern maintainability and 
reliability differences between designs that are as close as 
the machine designs of this thesis, so these two charac- 
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teristics were not made part of the objective function. 

Cost is a measure that should be part of this sort of op¬ 
timization. The cost of the ship system, quantified in the 
ship displacement, was used in the final recommendations for 
the transmission sub-system. In the case of permanent mag¬ 
net machines, the relative costs of magnet material and mag¬ 
net steel were included in the objective function. 

Weight, volume, and efficiency were made an explicit 
part of the objective function for the computer design of 
the machines. An Effective Weight was calculated for every 
machine design. The design with the lowest effective weight 
was the "best" within its class. The generic objective 
function is 

Effective Weight = weight + ke*(1-efficiency) + kv*volume 

where ke and kv are weighting factors for efficiency and 
volume, respectively. The weighting factors were obtained 
from changes in ship displacement for marginal changes in 
efficiency and volume of the transmission. They were 
modified to reflect the actual designs resulting from the 
process. 
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Chapter Two. General Considerations 


Only steady-state behavior of electric machines was 
modeled. The modeling of dynamic behavior is very difficult 
and was not viewed as being within the design problems posed 
by this thesis. The changes in machine design necessary to 
solve dynamic instabilities, etc., are much smaller than the 
approximate nature of the algorithms used here. 

All derivation work was performed without specifying 
the number of winding turns or the number of rotor or stator 
slots. The only exception to this was the case of induction 
machines, where an arbitrary number of rotor slots was 
selected. This selection was necessary for the calculation 
of the equivalent circuit components. The number of rotor 
slots chosen, 71, was a number designed not to induce pole 
harmonics. Since no turn numbers were specified, units in¬ 
clude volts-per-turn, ampere-turns, and impedance-per-turns- 
squared. Power is measured in watts. 

All machines used as their synchronous frequency the 
maximum allowed by the particular combination of pole pairs 
and shaft rpm. Developmental work showed that the optimiza¬ 
tion algorithm converged to the highest frequencies, so the 
algorithm now starts at the highest possible frequencies. 

Up to six pole pairs were used in the higher rpm 
machines and up to 25 pole pairs in the 180 rpm machines. 
Diminishing improvements in volume, weight, and efficiency 
show up at half these limits. 

The random number generator was taken from Kelley and 
Pohl [14], with one change. After every run of each 
program, the random number generator seed is stored. This 
means the sequence of pseudo random numbers is not repeated 
until the full range possible has been used. For the 
machines, it gives differences at every run and means the 
multidimensional variable space is more fully explored. 
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2.1. Optimization method 


The chosen optimization method is a combination of the 
Monte Carlo and steepest-descent schemes. A design point is 
established by randomly selecting machine geometric 
parameters, subject to constraints. Ten random steps are 
taken around the design point, in all variable directions. 
The effective weight of each random step is evaluated and 
compared with that of the design point. The best of the 
eleven is designated the new design point. More random 
steps are taken, and the process continues until no more im¬ 
provement is seen in effective weight. At that point, the 
size of the random steps is halved, and the process repeats 
itself, with the step size continually halved (up to ten 
times). The best effective weight is the index to the best 
design. The number of original design points used in any 
particular run of a program is under user control. If ten 
original design points are desired, the algorithm will look 
at over a thousand designs. 

The purpose of having original design points is to 
start the optimization process in different sectors of the 
multidimensional variable space. In this fashion, the op¬ 
timization process zeroes in on several local "best" points. 

The variables that are randomly selected include stator 
current density, rotor radius, air gap dimension, stator 
slot space factor, and rotor slot space factor. The back 
iron dimension (the iron behind the stator teeth) is sized 
to handle a saturation level of flux. The stator slot depth 
is originally sized as a random fraction of the back iron 
dimension, and the rotor slot depth is originally a random 
fraction of the rotor radius. 

Only steady-state behavior was modeled in this algo¬ 
rithm. Dynamic modeling may or may not show different op¬ 
timum configurations for machine types. 
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2.2. Constraints 


The constraints placed on the optimization process are 
listed in the following table. The most difficult con¬ 
straint to satisfy while still achieving a valid design was 
the rotor current density constraint in induction machines. 
Only a few valid designs were achieved in induction machines 
using the above algorithm, leaving some question about the 
application of the algorithm in the case of induction 
machines. Only those induction machine designs in which 
there was reasonable confidence were included in the thesis 
analysis. 


Table 3. Optimization constraints during machine design 


Minimum air gap flux density 
Maximum (saturation) flux density 
Maximum rotor radius 
Maximum rotor tip speed 
Stator and rotor space factor 
Maximum rotor slot depth 
Maximum synchronous reactance: 
synchronous machines 
permanent magnet machines 
Power factor 


1.05 tesla rms 
1.5 tesla rms 
2.0 meters 
200.0 meters/sec 
0.35 

33% of rotor radius 

2.0 per unit 
3.0 per unit 
0.8 


The magnet steel chosen was 26 gauge M19. Its magnetic 
properties were found in USX technical data [15). It has 
been observed that saturation flux levels in electric 
machines occur first where the area perpendicular to the 
flux path is the smallest. If the back iron dimension is 
made appropriately large, this saturation will first occur 
in the teeth, as is desirable. Accordingly, the back iron 
dimension was set to 

Br r 

dcore = - 

Bsat p 

where Br is the radial air gap flux density, r is the rotor 
radius, Bsat is the saturation flux density, and p is the 
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number of pole pairs in the machine. This equation is 
derived from Gauss’ Law. 

If electric machines are to be installed in a ship, th 
78 ey obviously will need to fit into the space designated 
for them. An electric machine with a two meter rotor radius 
is at least thirteen feet in envelope diameter. This is a 
very large machine to install in a machinery space where 
volume is already at a premium. The rotor radius limit de¬ 
scends from the physical ability to fit an electric machine 
in a ship. 

The tip speed limitation represents the physical limit 
on material strength with regard to the rotor conductors. 
Rotor conductors may break free from the rotor at higher 
tangential velocities than this limit. The magnitude of the 
limit was taken from a tip-speed-limited, 3600 rpm, two-pole 
turbogenerator, and was verified against standard Navy 
design practices. Several as-built electric machines were 
analyzed and this number seemed to fit their characteristics 
well. The tip speed limit arises when choosing a rotor 
radius (given a particular frequency and number of poles), 
and is less stringent a constraint than maximum rotor 
radius. 

Thermal considerations are often extremely important in 
machine design. The heat build-up in electric machines has 
led to many cooling schemes over the years, including 
natural convection, forced air cooling, and hydrogen cool¬ 
ing. One of the latest methods is liquid cooling of the 
stator and rotor conductors through cooling passages through 
the copper itself. This has been made possible by better 
de-ionizing methods for cooling fluid and better rotating 
seals for the rotor. Naturally, the cooling passages and 
insulation limit the amount of copper area in a slot cross 
section. The copper area in a typical conductor bar was 
measured and found to be about thirty-five percent of the 
bar cross section. This number was used for the stator and 
rotor slot space factors. 
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Other thermal considerations must be made for permanent 
magnet materials, which suffer a degradation in flux as tem¬ 
peratures rise. Flux loss rises slowly with increasing tem¬ 
perature until about 100°C. Above 100°C, flux loss is more 
rapid. An assumption in the design of these machines is 
that there will be sufficient cooling in the operating space 
to limit ambient temperature to about 80°C. This, combined 
with the machine liquid cooling, should keep flux loss to a 
minimum. Transmission lines were assumed to function satis¬ 
factorily at the same temperature. 

Insulation also has a thermal rating. No insulation 
class was specified in this thesis but a typical insulation 
used in electric generators by the Navy is Class F. For 
this class, a permissible rise of 100°C over an ambient tem¬ 
perature of 50°C is standard, but lesser insulation classes 
must run cooler. If a machine must be designed with a les¬ 
ser insulation class, the consequent lower temperatures may 
result in a quieter machine and longer machine and insula¬ 
tion life. It probably will be larger than a machine with a 
greater class of insulation. The temperatures quoted above 
are at hot spots, not in the bulk of the machine. 11 

Along with the reduction in copper area for liquid 
cooling, a maximum current density was imposed. The copper 
losses, in the form of heat, have to be removed by the cool¬ 
ing fluid. There is a tradeoff between the size of the 
cooling passages, the allowable current density, and the 
rating of the machine. Twelve million amperes per square 
meter equates to forty amperes in a twelve gauge copper 
wire. 

Rotor slots were constrained to no more than one third 
of the rotor radius. Some reasonable shaft diameter is 
required to transmit the mechanical torque. Stator slots 
were allowed to grow as needed to meet the stator current 


11. Private communication, D.F. Schmucker, Naval Sea Systems 
Command, Code 56Z31, 11 December 1986. 
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density limit. 

Synchronous reactance limits were taken from as-built 
machines. A power factor of 0.8 was used for all machines, 
though one researcher has indicated a power factor of 1.0 
might be best for permanent magnet machines. 12 

2.3. Geometric considerations 

End turns were modeled as described in Appendix B. 

While not exactly as machines are constructed, this model 
gives reasonable results. A length allowance equal to one 
rotor diameter on either end of the active length of the 
machine was made to allow for containment of the end turns. 

Fractional slot pitches were not considered. A stator 
winding pitch of 0.8 was assumed, resulting in the elimina¬ 
tion of the fifth harmonic from the steady-state output 
waveforms. 

Once weight and volume were calculated, an extra ten 
percent was added to allow for the frame and foundation of 
the machine. The calculated weight included an additional 
three percent of the rotor weight to allow for bearings. It 
is the final envelope weight and volume that were used in 
the decision process. 

2.4. Efficiency and losses 

The general equation for efficiency is 

minpwr 

efficiency = - 

minpwr + ph + pe + i2r + i2rr 

where minpwr is the minimum mechanical power expected of the 
machine, ph is the hysteresis loss, pe is the eddy current 


12. Robey, H.N., "Permanent Magnet Machine Technology 
Assessment,” DTNSRDC Report TM-27-80-87, September 1980. 
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loss, i2r is the stator copper loss, and i2rr is the rotor 
copper loss. This formulation is for a motor, but the ef¬ 
ficiency calculated will not be significantly different if 
the machine is a generator. 

Hysteresis and eddy current losses arise from currents 
circulating within the magnet steel that forms the rotor and 
stator. They are two different mechanisms and depend on the 
metallurgy of the steel. 

Eddy currents are a result of the time-varying magnetic 
fields within the machine, and they oppose the change in 
flux density within the machine. Eddy current losses in¬ 
crease as the square of the electrical frequency of the 
machine and also as the square of the peak flux density. 

One method of lessening eddy current losses is to use thin 
laminations to build up the rotor and stator. If the var¬ 
nish used on the laminations is sufficiently insulating, the 
eddy currents are limited to azimuthal circulation. Axial 
circulation is practically zero because of the small lamina¬ 
tion thickness. 

Hysteresis losses are inherent to magnetic materials, 
and are proportional to the total volume of the material, 
the area of the hysteresis loop, and the machine electrical 
frequency. 

USX has developed equations to calculate eddy current 
and hysteresis losses. They are 


ph 


0^01445j3_f_Br_Hc 


and pe 


0.4818 H Em* t 2 f* 
rho D 


In these equations, (i is the hysteresis loss factor (the 
ratio of the actual hysteresis losses to the area of a 
square hysteresis loop passing through Br and He), f is the 
frequency in Hertz, Br is the residual induction in 
kilogauss, He is the coercive force in oersteds, D is the 
density in grams per cubic centimeters, rho is the electri¬ 
cal resistivity of M19 in microhm-cm, fl is the anomalous 
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loss factor. The losses are in watts per pound of material. 
The numerical factors at the beginning of both equations 
were altered to reflect the use of SI units. The factors <i , 
fl, Br, He, rho, and D change with the type of magnetic 
material used. The thickness of the laminations, t, is 
0.014 inches for 26 gauge steel. These equations do not 
reflect some variations caused by differences in silicon 
content and differences in processing treatments leading to 
variations in grain size and crystallographic texture. 
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Chapter Three. Synchronous machines 


Synchronous machines operate because of an interaction 
between stator and rotor flux waves. The rotor flux wave is 
developed by a field winding. The stator flux wave is 
developed by an armature winding. These two waves try to 
align themselves, which is how the machine action is 
produced. In the case of a motor, the armature wave is 
"rotating” around the periphery of the stator bore because 
of the 120° separation between the three phases. A rotating 
action ensues. In the case of a generator, the rotation is 
provided by a prime mover, such as a gas turbine, and volt¬ 
age is induced in the stator phases. 

Synchronous machines operate at a steady-state shaft 
speed specified by the number of poles and the electrical, 
or synchronous, frequency. This synchronous speed is main¬ 
tained despite changes in load. This feature makes 
synchronous machines attractive for applications where speed 
control is important. Shaft rpm = (120*frequency)/(poles). 

A derivation of the equations of synchronous machines 
and the computer modeling program are presented in 
Appendix B. 

3.1. Assumptions 

The rotors of synchronous machines may exhibit 
saliency, or may be smooth cylindrical rotors. The dif¬ 
ferences in properties and parameters among salient and 
round-rotor machines amount to only a few percent* *. The 
approximate nature of the modeling means the saliency ef¬ 
fects will not be important. Therefore, no special provi¬ 
sion for salient rotors were made. 

The ships used in this thesis have displacements of 


13. Fitzgerald, Kingsley, and Umans, Electric Machinery , 
McGraw-Hill, New York, 1983. 
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about 5000 LT ( 1 LT = 2240 lbs). Their fuel load is deter¬ 
mined by the required range and transmission efficiency at 
endurance speed. It was first estimated that a one-percent 
increase in transmission efficiency would reduce total ship 
displacement by about 89 LT. The efficiency factor (ke = 
90,000 kg/percent) corresponds to this 89 LT change in ship 
full load displacement. (89 LT x 2240 lbs/LT x 2.205 kg/lb) 
When this produced machines with efficiencies about 95%, it 
was doubled to 180,000. Obviously, this factor may be ad¬ 
justed to any level. The volume efficiency factor (kv = 
1286.1 kg/m 3 ) corresponds to the density of a LM-2500 gas 
turbine module. These factors were used throughout the 
thesis. 

3.2. Machine description 

Machines with shaft speeds of 1800, 2400, 3000, 3600, 
and 7200 rpm, with the number of pole pairs varying from one 
to six, were modeled. Also, 180 rpm machines using from one 
to twenty-five pole pairs were modeled. This provided a 
good coverage of the variable space. 

3.2.1. Efficiency 

Synchronous machine efficiency at full load was about 
98.5% for the higher rpm machines, while the 180 rpm 
machines hovered around 93% to 94% efficiency. The number 
of pole-pairs seemed to have little effect in the high rpm 
machines, but there was an "arch" in the efficiency curve of 
the 180 rpm machines, peaking at 97% with 36 poles. Though 
not fully understood, the 24- and 26-pole 180 rpm machines 
had very low efficiencies. Generally speaking, machine ef¬ 
ficiency was higher when rpm was higher. This was an ex¬ 
pected result. Off-design-point efficiency was good for the 
higher rpm machines, but bad for the 180 rpm machines. (See 
the discussion in Chapter Seven.) 
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3.2.2. Weight and volume 


Weight and volume increased almost linearly with the 
number of pole-pairs in the 180 rpm machine. This is almost 
certainly a function of tip-speed limitations, as the maxi¬ 
mum tip speed in a machine is a function of the number of 
poles in the machine. When the rotor radius is limited, the 
machine must grow in length to develop enough torque. The 
weight and volume of these machines were much higher than 
for the higher rpm machines. 

The higher rpm machines saw significant decreases in 
weight and volume when the number of pole-pairs increased 
from one to two. There also was an observable increase in 
weight and volume as the number of pole-pairs further in¬ 
creased. As rpm increased, the machines became smaller and 
lighter. 
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EFFICIENCY 


Figure 1. High rpm synchronous machine efficiency 


SYNCHRONOUS MOTOR EFFICIENCY 
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EFFICIENCY 


Figure 2. 180 rpm synchronous machine efficiency 


SYNCHRONOUS MOTOR EFFICIENCY 


DIRECT-DRIVE, 180 RPM. 25775 HP 
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VOLUME (cubic meters) 


Figure 3. High rpm synchronous machine volume 


SYNCHRONOUS MACHINE VOLUME 
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VOLUME (M~3) 


Figure 4. 180 rpm synchronous machine volume 


SYNCHRONOUS MOTOR VOLUME 


DIRECT-DRIVE, 180 RPM, 25775 HP 
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WEIGHT 

(Thousands) 


Figure 5. High rpm synchronous machine weight 


SYNCHRONOUS MOTOR WEIGHT 


25775 HP 
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WEIGHT (KG) 
(Thousands) 


Figure 6. 180 rpm synchronous machine weight 


SYNCHRONOUS MOTOR WEIGHT 


DIRECT-DRIVE, 180 RPM, 25775 HP 
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3.3. Machine characteristics 


The following tables give machine characteristics for 
many of the higher rpm and the 180 rpm synchronous machines 
designed for this thesis. 

The stator slot factor tended to increase to the limit 
of 0.75, while the rotor slot factor moved around the value 
0.58 quite a bit. This demonstrates a partial limit on the 
depth of the stator slots (to the same dimension as the back 
iron depth). The overall diameter was limited and more 
stator slot area was needed to develop the required power. 
Similarly, the stator current density converged to its cook¬ 
ing limit. 

The longest higher rpm machine was 5.42 meters, while 
the largest overall diameter was 0.85 meters. Machines of 
this size will cause no difficulties when placed in the 
machinery spaces of most ships. The 180 rpm machines are 
larger, typically less than 6 meters long (discounting one 
16 meter machine) and 1.8 meters in diameter. They are also 
good candidates for ship systems. 
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Table 4. Characteristics of 1800 rpm synchronous machines 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
shaft rpa 
stator current density 
synchronous frequency 
rotor radius 
gap dimension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall length 
overall diaieter 


1 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

1800 

1800 

1800 

8.00E+06 

1.20E+07 

1.20E+07 

30 

60 

90 

0.1720 

0.2077 

0.2603 

0.0158 

0.0440 

0.0250 

0.1204 

0.0727 

0.0607 

0.0772 

0.0722 

0.0607 

0.0401 

0.1859 

0.0747 

0.653 

0.750 

0.750 

0.650 

0.447 

0.687 

2.781 

1.522 

1.509 

20835.40 

9827.96 

8766.09 

16926.49 

10677.39 

17111.27 

3033.2 

3826.7 

9198.8 

275834.4 

374646.5 

293335.5 

0.985 

0.980 

0.984 

4.670 

1.794 

1.498 

1.19E+07 

1.06E+07 

1.17E+07 

4.40E+06 

6.39E+06 

7.00E+06 

1.99 

0.86 

0.87 

2.71 

1.66 

1.67 

5.42 

2.80 

2.64 

0.77 

0.79 

0.81 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

1800 

1800 

1800 

1.20E+07 

1.20E+07 

1.20E+07 

120 

150 

180 

0.2905 

0.3267 

0.3189 

0.0124 

0.0090 

0.0092 

0.0508 

0.0457 

0.0372 

0.0508 

0.0457 

0.0372 

0.1333 

0.0907 

0.0289 

0.749 

0.750 

0.748 

0.659 

0.449 

0.628 

1.541 

1.697 

1.709 

9186.03 

8661.17 

8853.21 

19469.83 

31370.98 

43892.02 

13955.6 

28107.7 

47191.5 

164289.8 

173116.1 

269901.7 

0.990 

0.988 

0.982 

1.513 

1.348 

1.741 

4.58E+-06 

8.34E+06 

1.89E+07 

2.43E+06 

4.24E+06 

1.19E+07 

1.11 

1.20 

0.78 

1.89 

1.97 

1.60 

2.72 

2.69 

3.05 

0.81 

0.85 

0.80 
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Table 5. Characteristics of 2400 rpm synchronous machines 


nuiber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
shaft rpa 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 
overall length 
overall diaaeter 


1 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

2400 

2400 

2400 

1.15E+07 

1.20E+07 

1.19E+07 

40 

80 

120 

0.1563 

0.1986 

0.2388 

0.0273 

0.0164 

0.0089 

0.1094 

0.0695 

0.0557 

0.0929 

0.0695 

0.0557 

0.0416 

0.0215 

0.0088 

0.621 

0.750 

0.750 

0.676 

0.749 

0.750 

1.611 

1.112 

1.133 

11619.57 

7204.42 

6677.15 

119B1.17 

14511.75 

21941.41 

2862.6 

6934.5 

15727.2 

444605.4 

427167.9 

482231.4 

0.977 

0.977 

0.974 

2.397 

1.709 

1.552 

2.11E+07 

3.11E+07 

5.39E+07 

7.75E+06 

1.28E+07 

2.11E+07 

2.00 

1.69 

1.83 

2.72 

2.43 

2.56 

3.13 

2.57 

2.54 

0.77 

0.71 

0.72 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

2400 

2400 

2400 

1.20E+07 

1.20E+07 

1.07E+07 

160 

200 

240 

0.2890 

0.3141 

0.3397 

0.0135 

0.0109 

0.0066 

0.0506 

0.0440 

0.0396 

0.0505 

0.0434 

0.0396 

0.0563 

0.1047 

0.1040 

0.748 

0.750 

0.750 

0.432 

0.453 

0.327 

1.329 

1.436 

1.633 

6733.44 

7106.43 

7727.01 

27615.01 

32769.75 

47103.39 

26391.9 

39147.9 

67525.7 

232958.2 

143185.6 

114790.8 

0.985 

0.989 

0.988 

1.152 

1.144 

1.223 

1.73E+07 

7.83E+06 

7.99E+06 

9.64E+06 

4.55E+06 

4.32E+06 

1.01 

0.92 

1.07 

1.80 

1.72 

1.85 

2.36 

2.44 

2.61 

0.81 

0.82 

0.85 
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Table 6. Characteristics of 3000 rpm synchronous machines 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
shaft rpa 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 

overall length 
overall diaieter 


1 

n 

L. 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

3000 

3000 

3000 

1.20E+07 

1.20E+07 

1.20E+07 

50 

100 

150 

0.1494 

0.1890 

0.2304 

0.0390 

0.0172 

0.0101 

0.1046 

0.0662 

0.0538 

0.0862 

0.062S 

0.0537 

0.0904 

0.0879 

0.0372 

0.750 

0.750 

0.750 

0.598 

0.494 

0.724 

1.205 

0.963 

0.972 

8815.49 

6498.77 

5801.99 

9574.79 

14072.41 

21089.10 

2859.6 

8405.7 

18895.4 

380285.7 

212641.9 

183124.5 

0.980 

0.988 

0.989 

1.670 

1.657 

1.363 

1.53E+07 

1.13E+07 

1.38E+07 

6.03E+06 

5.23E+06 

6.09E+06 

1.80 

1.41 

1.53 

2.53 

2.16 

2.27 

2.42 

2.48 

2.32 

0.76 

0.67 

0.70 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

3000 

3000 

3000 

1.20E+07 

1.10E+07 

1.18E+07 

200 

250 

300 

0.2432 

0.2872 

0.2974 

0.0052 

0.0068 

0.0088 

0.0426 

0.0402 

0.0347 

0.0422 

0.0402 

0.0347 

0.0151 

0.0784 

0.0427 

0.750 

0.750 

0.750 

0.730 

0.255 

0.750 

0.991 

1.203 

1.236 

5725.40 

6176.39 

6143.64 

31331.61 

41076.93 

46178.46 

37429.9 

61340.0 

82749.7 

204027.7 

140446.5 

134484.2 

0.986 

0.988 

0.986 

1.590 

1.306 

1.310 

2.42E+07 

1.40E+07 

1.06E+07 

9.63E+06 

7.44E+06 

6.96E+06 

1.78 

1.11 

0.69 

2.51 

1.89 

1.52 

2.58 

2.48 

2.53 

0.67 

0.75 

0.75 
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Table 7. Characteristics of 3600 rpm synchronous machines 


nuAber of pole pairs 
power, hp 
efficiency factor 
volute factor 
shaft rp# 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
full load current density 
no load current density 
xs/turns-squared, p.u. 
internal volts/turn, p.u. 

overall length 
overall diaaeter 


i 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

3600 

3600 

3600 

1.16E+07 

1.20E+07 

1.20E+07 

60 

120 

180 

0.1413 

0.1824 

0.2191 

0.0139 

0.0120 

0.0123 

0.0989 

0.0638 

0.0511 

0.0492 

0.0638 

0.0486 

0.0496 

0.0541 

0.0572 

0.748 

0.750 

0.750 

0.733 

0.409 

0.528 

1.233 

0.814 

0.870 

9385.14 

5377.73 

5216.22 

14030.95 

15587.02 

22251.57 

5028.6 

11172.5 

23924.3 

190131.9 

228605.1 

165930.4 

0.989 

0.9B7 

0.989 

3.258 

1.493 

1.372 

9.15E+06 

1.94E+07 

1.31E+07 

3.38E+06 

7.41E+06 

6.99E+06 

1.99 

1.89 

1.10 

2.71 

2.62 

1.88 

3.88 

2.27 

2.30 

0.61 

0.64 

0.66 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

3600 

3600 

3600 

1.19E+07 

1.20E+07 

1.20E+07 

240 

300 

360 

0.2736 

0.2846 

0.2877 

0.0178 

0.0209 

0.0124 

0.0479 

0.0398 

0.0336 

0.0479 

0.0398 

0.0336 

0.0298 

0.0595 

0.0385 

0.750 

0.750 

0.750 

0.679 

0.691 

0.654 

1.065 

1.127 

1.102 

4979.26 

5177.33 

5190.73 

31042.24 

36284.98 

48590.52 

44501.0 

65021.1 

104486.5 

251863.9 

196473.3 

184008.6 

0.983 

0.985 

0.983 

0.890 

0.975 

1.166 

2.47E+07 

1.45E+07 

1.73E+07 

1.61E+07 

1.12E+07 

1.28E+07 

0.71 

0.43 

0.48 

1.54 

1.30 

1.35 

2.06 

2.20 

2.37 

0.77 

0.77 

0.73 



Table 8. Characteristics of 7200 rpm synchronous machines 


nunber of pole pairs 

i 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

180000 

voluae factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

1236.1 

shaft rp* 

7200 

7200 

7200 

7200 

7200 

7200 

stator current density 

1.20E+07 

1.20E+07 

1.16E+07 

1.20E+07 

1.20E+07 

1.20E+07 

synchronous frequency 

120 

240 

360 

480 

600 

720 

rotor radius 

0.1159 

0.1710 

0.1958 

0.2044 

0.2219 

0.1932 

gap dimension 

0.0096 

0.0093 

0.0053 

0.0093 

0.0046 

0.0055 

back iron depth 

0.0812 

0.0599 

0.0457 

0.0358 

0.0311 

0.0225 

stator slot depth 

0.0405 

0.0452 

0.0457 

0.0358 

0.0311 

0.0225 

rotor slot depth 

0.0329 

0.0189 

0.0105 

0.0714 

0.0138 

0.0292 

stator slot factor 

0.750 

0.750 

0.749 

0.750 

0.750 

0.750 

rotor slot factor 

0.511 

0.617 

0.661 

0.633 

0.468 

0.570 

envelope voluie 

0.711 

0.544 

0.525 

0.544 

0.572 

0.564 

envelope weight 

5341.21 

3476.08 

2971.92 

3136.32 

2950.31 

3427.07 

hysteresis loss 

17305.B6 

22097.70 

28944.03 

30808.73 

50465.33 

64447.34 

eddy current loss 

12404.5 

31678.5 

62239.7 

88332.6 

180863.2 

277168.2 

stator copper loss 

149625.0 

159505.1 

162532.1 

74543.2 

143253.5 

117812.4 

full load efficiency 

0.991 

0.989 

0.987 

0.990 

0.981 

0.977 

active length 

2.863 

1.212 

0.972 

1.081 

1.082 

1.952 

full load current density 

1.65E+07 

2.71E+07 

3.92E+07 

7.81E+06 

3.70E+07 

1.72E+07 

no load current density 

6.14E+06 

1.16E+07 

1.46E+07 

5.00E+06 

1.99E+07 

1.28E+07 

xs/turns-squared, p.u. 

1.97 

1.59 

1.97 

0.74 

1.08 

0.48 

internal volts/turn, p.u. 

2.69 

2.33 

2.69 

1.56 

1.86 

1.34 

overall length 

3.36 

1.93 

1.78 

1.94 

1.99 

2.75 

overall diaieter 

0.49 

0.57 

0.58 

0.57 

0.58 

0.49 
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Table 9. Characteristics of 180 rpm synchronous machines 


nuiber of pole pairs 

i 

2 

3 

4 

5 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

voiuae factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

180 

180 

180 

180 

180 

stator current density 

1.18E+07 

1.12E+07 

1.20E+07 

9.25E+06 

1.20E+07 

synchronous frequency 

3 

6 

9 

12 

15 

rotor radius 

0.2641 

0.3738 

0.4588 

0.5869 

0.5507 

gap disension 

0.0331 

0.0978 

0.0723 

0.0317 

0.0198 

back iron depth 

0.1849 

0.1308 

0.1071 

0.1027 

0.0771 

stator slot depth 

0.0592 

0.1308 

0.1071 

0.1027 

0.0580 

rotor slot depth 

0.0500 

0.1867 

0.1934 

0.0565 

0.1079 

stator slot factor 

0.741 

0.750 

0.750 

0.742 

0.750 

rotor slot factor 

0.558 

0.571 

0.422 

0.535 

0.287 

envelope voiuae 

16.305 

9.329 

9.048 

11.378 

10.306 

envelope weight 

119110.5 

56300.0 

48896.8 

56398.0 

56899.5 

hysteresis loss 

10465.3 

6915.1 

9833.5 

18267.4 

23999.1 

eddy current loss 

187.5 

247.8 

528.6 

1309.4 

2150.3 

stator copper loss 

2.25E+06 2.06E+06 

1.76E+06 

1.34E+06 

1.18E+06 

full load efficiency 

0.895 

0.903 

0.916 

0.934 

0.941 

active length 

14.913 

3.135 

2.590 

2.376 

3.708 

full load current density 

1.53E+07 1.25E+07 

1.33E+07 

1.91E+07 

1.42E+07 

no load current density 

5.61E+06 6.14E+06 

7.25E+06 

8.93E+06 

7.27E+06 

xs/turns-squared, p.u. 

2.00 

1.28 

1.04 

1.38 

1.19 

internal volts/turn, p.u. 

2.72 

2.04 

1.83 

2.14 

1.96 

overall length 

16.10 

5.02 

4.71 

4.85 

5.99 

overall diaaeter 

1.08 

1.47 

1.49 

1.65 

1.41 
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Table 10. Characteristics of more 180 rpm synchronous 
machines 


nusber of pole pairs 

b 

7 

8 

9 

10 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

volute factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

180 

180 

180 

180 

180 

stator current density 

1.20E+07 

1.20E+07 

1.20E+07 

1.20E+07 

1.20E+07 

synchronous frequency 

18 

21 

24 

27 

30 

rotor radius 

0.6176 

0.5729 

0.7152 

0.6943 

0.7296 

gap diiension 

0.0647 

0.0335 

0.0229 

0.0138 

0.0396 

back iron depth 

0.0720 

0.0573 

0.0626 

0.0540 

0.0511 

stator slot depth 

0.0720 

0.0569 

0.0626 

0.0540 

0.0509 

rotor slot depth 

0.1652 

0.0855 

0.2046 

0.0740 

0.1420 

stator slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

rotor slot factor 

0.368 

0.750 

0.374 

* 0.2B3 

0.615 

envelope volute 

11.693 

10.486 

12.870 

12.385 

14.273 

envelope weight 

48739.8 

53234.1 

52286.9 

53189.2 

56440.5 

hysteresis loss 

22385.0 

27947.6 

31793.9 

42718.3 

42086.1 

eddy current loss 

2406.8 

3505.7 

4557.9 

6889.4 

7541.6 

stator copper loss 

1.B6E+06 

1.24E+06 

7.69E+05 

1.14E+06 

1.11E+06 

full load efficiency 

0.911 

0.938 

0.960 

0.942 

0.943 

active length 

2.225 

3.418 

2.045 

2.549 

2.365 

full load current density 

1.74E+07 

1.11E+07 

7.09E+06 

1.98E+07 

9.69E+06 

no load current density 

1.29E+07 

7.98E+06 

4.18E+06 

1.07E+07 

7.76E+06 

xs/turns-squared, p.u. 

0.48 

0.53 

0.89 

1.08 

0.36 

internal volts/turn, p.u. 

1.35 

1.39 

1.69 

1.86 

1.25 

overall length 

4.95 

5.84 

5.00 

5.38 

5.44 

overall diateter 

1.65 

1.44 

1.73 

1.63 

1.74 
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3.4. Verification 


Data on a large turbogenerator, Big Sandy Unit Two, was 
available in reference sixteen; it was used to verify the 
synchronous machine design program. Big Sandy is rated at 
907,000 kVA, power factor 0.9, at a rated voltage of 26 kV. 
When the machine parameters were input to the synchronous 
design program, it produced a machine very close to Big 
Sandy. It was judged that the design program would yield 
good results. 
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Chapter Four. Permanent Magnet Machines 


Permanent magnet machines are very similar to 
synchronous machines. The main difference lies in the 
method used to produce the field flux wave. Instead of a 
field current causing the wave, permanent magnets provide 
the flux. No exciter is used with these machines and there 
are no brushes. 

4.1. Magnet material 

Many different elements may be used to manufacture per¬ 
manent magnets. Past designs used ceramics, aluminum- 
nickel-cobalt-iron-titanium (AlNiCo), and samarium-cobalt 
(SmCo). However, ceramic magnets do not produce sufficient 
residual flux (see Appendix C for an explanation of terms), 
and any magnet based on cobalt is high in cost and may be in 
limited supply. Recently, magnets of neodymium-iron-boron 
(NdFeB), have entered the marketplace. None of the con¬ 
stituents of the NdFeB magnets are strategic materials; it 
is expected that availability and cost will improve. 

NdFeB magnets have a high Maximum Energy Product (MEP) 
that may be used advantageously by the machine designer. 

Data from Sumitomo Special Metals [17] indicates their 
NEOMAX line have MEPs as high as 37 MGOe, which is higher 
than the 30 MGOe of the SmCo magnets marketed as REC-30 by 
TDK Corporation [18]. High MEP is not the only criteria for 
magnet selection; flux stability, cost-to-performance ratio, 
ease of machine assembly, and other characteristics may en¬ 
ter the decision process. This study needed the best per¬ 
formance of its machines, so NdFeB was selected on the basis 
of its high MEP. Thermal stability was assumed to be satis¬ 
factory if the thermal considerations in Chapter Two were 
met. 

Cost has been mentioned above. Magnet material is sig¬ 
nificantly more expensive than magnet steel and copper. 
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Magnets in the quantity used by a large production run of 
25,000 hp machines might cost as much as $120 per pound, 14 
compared with the 58<|j per pound of M19 steel. 15 Obviously, 
permanent magnet machines will be more expensive, but the 
cost of magnet materials may be made part of the optimiza¬ 
tion process. The degree of magnet overhang, discussed in 
Appendix C, also affects cost. 

4.1.1. Magnet cost factor 

A change to the objective function was made to incor¬ 
porate the cost of magnet material relative to magnet steel. 
The ratio of the above costs was taken and the result called 
the magnet cost factor, km. The objective function was 
modified to 

Effective weight = (weight + kalaagnet weight)) ♦ ke(l - effcy) + kv(voluae + ka(aagnet voluae)) 

An initial value for km of 170 was used, and several 
machines designed. Then a value of km = 25 was tried. The 
machines with km = 170 indeed had less magnet material in 
them, but at a cost. The change to km = 25 resulted In a 
larger machine (23.5%) with a lower stator current density, 
20% more magnet material, and about a 1.5% increase in 
machine efficiency. That 1.5% translates to a lot of fuel 
aboard a ship, so it was decided to use km = 25. The extra 
magnet material will add about $22,000 per machine. 

4.2. Assumptions 

The largest obstacle to assembling high-power permanent 
magnet machines is their inherent residual magnetism. If 

14. Estimate by Mr. Yokokura, President of Sumitomo Special 
Metals of America. 

15. Book price for 26 gauge M19 steel from Mr. Dagg of USX. 


51 




















the magnets possess all of their properties at assembly, it 
will be extremely difficult to place the rotor (which con¬ 
tains the magnets) inside the stator. The rotor would be 
strongly attracted to the iron of the stator. Of course, 
the magnets may be magnetized after machine assembly, but it 
may be difficult to achieve MEP without elevated tempera¬ 
tures inside the machine. The assumption is made here that 
the magnets will be magnetized prior to assembly. The 
detailed design of the machine will have to include con¬ 
sideration of the jigs and fixtures necessary for assembly. 

The only other assumption worthy of mention is that the 
load line may be modeled as described in Appendix C. 

4.3. Machine description 

The same rpm and pole-pair combinations were used as in 
the synchronous machines. The magnets on the rotor are ar¬ 
ranged in a cylindrical-wedge configuration, as shown in the 
figure below. 

Figure 7. End view of a permanent magnet machine 



M/lchNET 


MAGiNET STE.EU 

binder. 


The rotor slot factor (lr) used in the other types of 
machines is called here the magnet slot factor and refers to 
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the width of magnet per pole pitch. It may vary between 25% 
and 75% of the pole pitch, the same as for rotor slots in 
other machines. 

The rotor slot depth, ds, does not exist in this 
machine. In this case, the rotor radius, r, is added to lm, 
the magnet radial dimension, to find the actual width of the 
rotating core. (In the synchronous and induction machines, 
dr was included in the rotor radius.) 

4.3.1. Efficiency 

For the higher rpm machines, efficiency within a par¬ 
ticular rpm group decreased with an increasing number of 
poles. The most efficient machines, at about 99%, had four 
poles. This is higher than the synchronous or induction 
machines, largely because there are no rotor copper losses. 
The twelve-pole efficiencies were about 98%, which is not 
too large a spread. 

The 180 rpm machines had a fairly flat efficiency curve 
(excepting one anomaly) up to about a 28-pole machine, where 
efficiency started to vary widely. There, the conflict be¬ 
tween the number of pole-pairs and maximum rotor radius 
started to become significant. The flat efficiency was 
about 96.5%, which is less efficient than the higher rpm 
machines. 

4.3.2. Weight and volume 

The higher rpm machines had a general tendency toward 
lower weight and volume as rpm increased. Within an rpm 
group, the four- and six-pole machines had the lowest weight 
and volume. The smallest machines were larger than the 
synchronous machines. 
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EFFICIENCY 


Figure 8. High rpm permanent magnet machine efficiency 


PERMANENT MAGNET MOTOR EFFICIENCY 
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EFFICIENCY 


Figure 9. 180 rpm permanent magnet efficiency 


PERMANENT MAGNET MOTOR EFFICIENCY 


DIRECT-DRIVE, 180 RPM, 25775 HP 
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VOLUME (cubic meters) 


Figure 10. High rpm permanent magnet machine volume 


PERMANENT MAGNET MACHINE VOLUME 


25775 HP 



NUMBER OF POLE-PAIRS 
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VOLUME (M^3) 


Figure 11. 180 rpm permanent magnet machine volume 


PERMANENT MAGNET MOTOR VOLUME 


DIRECT-DRIVE. 180 RPM, 25775 HP 



57 



























































































WEIGHT (KG) 
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Figure 12. High rpm permanent magnet machine weight 


PERMANENT MAGNET MOTOR WEIGHT 


25775 HP 



NUMBER OF POLE-PAIRS 
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Figure 13. 180 rpm permanent magnet machine weight 


PERMANENT MAGNET MOTOR WEIGHT 


DIRECT-DRIVE, 180 RPM, 25775 HP 
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The weight and volume of the 180 rpm machines increased 
with the number of pole-pairs. Again, the tip speed limita¬ 
tion is linked to this increase, as with the synchronous 
machines. There was a wide variation around the general in¬ 
crease, and these machines are large relative to all others. 
They are not competitive as ship propulsion motors because 
of their size. 

4.4. Machine characteristics 

The following tables give the machine characteristics 
of the permanent magnet machines designed for this study. 

The stator slot factor, Is, tended to the maximum of 0.75 
for the higher rpm and 180 rpm machines. This occurred as 
the optimization algorithm tried to minimize envelope volume 
and weight. The rotor slot factor, lr, was a constant 0.378 
for all machines. The rotor slot factor was calculated to 
produce load-line (MEP) operating flux, as derived in Appen¬ 
dix B. With a different magnet material selection (and a 
consequent change in operating point flux), a different lr 
would have resulted. 

The magnet overhang tended toward the maximum limit, in 
an attempt to achieve the highest flux levels. Permanent 
magnet machines cannot rival the flux level produced by the 
field winding of a synchronous machine, but the optimization 
algorithm did its best. 

The per-unit synchronous reactance-per-turns-squared 
was limited to 3.0 in these machines. This reactance tended 
toward the limit, but was lower with an increasing number of 
pole-pairs in the 180 rpm machines. It was very difficult 
to achieve valid designs with xsmax = 2.0, as in the 
synchronous machines. If an xsmax greater than 2.0 is unac¬ 
ceptable, these machines will be less competitive. 

The amount of magnet material varied from 50 kg to a 
few hundred kg in the higher rpm machines to 700-4000 kg in 
the 180 rpm machines. The cost of 700 kg of NdFeB is about 
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Table 11. Characteristics o-f 1800 rpm magnet machine® 


nuiber of pole pairs 
power, hp 
efficiency factor 
volume factor 
•agnet factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
■agnet radial dimension 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
■agnet weight 
■agnet voluse 
xs/turns-squared, p.u. 
■agnet overhang 
air gap flux density 
overall length 
overall diameter 


1 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

1800 

1800 

1800 

7.49E+06 

1.20E+07 

7.32E+06 

30 

60 

90 

0.2887 

0.3396 

0.3643 

0.0326 

0.0382 

0.0067 

0.1332 

0.0870 

0.0530 

0.0687 

0.0742 

0.0456 

0.0342 

0.0401 

0.0070 

0.750 

0.750 

0.750 

0.378 

0.378 

0.378 

3.548 

2.558 

3.309 

24280.7 

12528.8 

19366.5 

20581.5 

20298.5 

52315.8 

3688.1 

7274.8 

28124.3 

171150.0 

265272.5 

93554.8 

0.990 

0.985 

0.991 

2.328 

0.876 

2.833 

405.184 

262.376 

157.355 

0.055 

0.035 

0.021 

2.994 

2.339 

2.996 

0.001 

0.298 

0.315 

0.619 

0.687 

0.642 

3.750 

2.548 

4.345 

1.047 

1.078 

0.939 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

1800 

1800 

1800 

9.14E+06 

8.00E+06 

1.04E+07 

120 

150 

180 

0.5964 

0.5544 

0.6068 

0.0188 

0.0067 

0.0123 

0.0742 

0.0502 

0.0482 

0.0353 

0.0263 

0.0277 

0.0198 

0.0071 

0.0130 

0.750 

0.750 

0.646 

0.378 

0.378 

0.378 

6.304 

5.847 

6.319 

20910.3 

24826.4 

20948.1 

80505.5 

121495.4 

124368.0 

57704.9 

108857.1 

133717.2 

98867.0 

69323.3 

85961.4 

0.988 

0.985 

0.982 

0.933 

1.888 

1.257 

273.763 

178.613 

246.460 

0.037 

0.024 

0.033 

1.165 

1.611 

0.920 

0.326 

0.331 

0.329 

0.722 

0.671 

0.700 

3.473 

4.160 

3.786 

1.450 

1.275 

1.390 
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Table 12. Characteristics of 2400 rpm magnet machines 


nuiber of pole pairs 
power, hp 
efficiency factor 
volume factor 
iagnet factor 
shaft rpi 
stator current density 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
•agnet radial diiension 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
iagnet weight 
iagnet voluse 
xs/turns-squared, p.u. 
iagnet overhang 
air gap flux density 
overall length 
overall diaieter 


1 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

2400 

2400 

2400 

12000000 

12000000 

12000000 

40 

80 

120 

0.3594 

0.0005 

0.4018 

0.0371 

0.1883 

0.0136 

0.1640 

0.0420 

0.0625 

0.0697 

0.0249 

0.0341 

0.0390 

0.1978 

0.0143 

0.348 

0.750 

0.750 

0.378 

0.378 

0.378 

4.497 

1.342 

2.727 

26548.1 

8797.0 

12335.3 

34028.8 

18160.1 

46136.0 

8130.4 

8677.9 

33069.4 

234324.1 

217110.2 

127316.4 

0.986 

0.987 

0.989 

1.536 

4.392 

1.291 

384.207 

10.204 

167.444 

0.052 

0.001 

0.023 

2.417 

0.282 

1.973 

0.018 

0.330 

0.331 

0.617 

0.635 

0.676 

3.278 

5.938 

3.010 

1.260 

0.511 

1.024 


4 

5 

6 

25775 

25775 

25775 

1B0000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

2400 

2400 

2400 

10202760 

11959030 

11997423 

160 

200 

240 

0.4408 

0.3361 

0.4760 

0.0088 

0.0032 

0.0165 

0.0504 

0.0290 

0.0379 

0.0325 

0.0234 

0.0300 

0.0092 

0.0033 

0.0173 

0.750 

0.749 

0.678 

0.378 

0.378 

0.378 

3.157 

2.349 

3.438 

13510.2 

13660.4 

11877.1 

68103.9 

86724.5 

91478.9 

65087.5 

103604.3 

131141.0 

90617.2 

112885.8 

91721.4 

0.988 

0.984 

0.984 

1.387 

3.061 

1.129 

128.106 

73.723 

217.676 

0.017 

0.010 

0.029 

1.965 

2.879 

0.751 

0.315 

0.338 

0.324 

0.672 

0.640 

0.692 

3.222 

4.431 

3.168 

1.065 

0.783 

1.121 
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Table 13. Characteristics of 3000 rpm magnet machines 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
aagnet factor 
shaft rpa 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
aagnet radial diaension 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
aagnet weight 
aagnet voluae 
xs/turns-squared, p.u. 
aagnet overhang 
air gap flux density 
overall length 
overall diaaeter 


1 

2 

3 

4 

5 

6 

25775 

25775 

25775 

25775 

25775 

25775 

180000 

180000 

180000 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

25 

25 

25 

25 

25 

25 

3000 

3000 

3000 

3000 

3000 

3000 

10639824 

9557320. 

12000000 

11999567 

12000000 

11400662 

50 

100 

150 

200 

250 

300 

0.2480 

0.2394 

0.2239 

0.0003 

0.3116 

0.4837 

0.0227 

0.0075 

0.0035 

0.1603 

0.0027 

0.0375 

0.1121 

0.0510 

0.0317 

0.0175 

0.0267 

0.0428 

0.0405 

0.0369 

0.0235 

0.0074 

0.0224 

0.0312 

0.0238 

0.0079 

0.0037 

0.1683 

0.0028 

0.0394 

0.750 

0.750 

0.750 

0.750 

0.747 

0.681 

0.378 

0.378 

0.37B 

0.378 

0.378 

0.378 

2.213 

1.834 

1.771 

2.143 

1.950 

3.643 

15051.7 

12394.7 

12381.0 

8254.8 

11507.2 

9968.3 

22692.8 

36848.2 

57339.4 

45284.3 

91118.4 

95620.6 

6777.4 

22010.1 

51374.8 

54098.3 136066.7 

171347.8 

157512.6 

108787.2 

132086.7 

170766.7 

97680.0 

74849.0 

0.990 

0.991 

0.988 

0.986 

0.983 

0.983 

2.395 

3.719 

5.490 

16.723 

3.005 

0.732 

250.152 

123.693 

91.981 

15.685 

55.765 

338.030 

0.034 

0.017 

0.012 

0.002 

0.008 

0.046 

2.934 

2.993 

2.997 

0.043 

2.991 

0.351 

0.007 

0.009 

0.326 

0.325 

0.311 

0.336 

0.618 

0.619 

0.627 

0.623 

0.637 

0.736 

3.574 

4.738 

6.415 

18.039 

4.274 

2.975 

0.847 

0.669 

0.565 

0.371 

0.727 

1.191 





































































Table 14. Characteristics of 3600 rpm magnet machines 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
■agnet factor 
shaft rp« 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
aagnet radial diaension 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
aagnet weight 
aagnet voluae 
xs/turns-squared, p.u. 
aagnet overhang 
air gap flux density 
overall length 
overall diaaeter 


1 

2 

3 

25775 

25775 

25775 

1B0000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

3600 

3600 

3600 

12000000 

12000000 

11966619 

60 

120 

180 

0.3044 

0.2596 

0.3006 

0.0228 

0.0088 

0.0121 

0.1454 

0.0574 

0.0462 

0.0332 

0.0442 

0.0411 

0.0240 

0.0092 

0.0127 

0.750 

0.553 

0.750 

0.378 

0.378 

0.378 

2.475 

1.602 

1.427 

15468.5 

9853.6 

7072.7 

29048.1 

36290.0 

37383.1 

10410.6 

26012.0 

40193.3 

163951.5 

114166.3 

104037.5 

0.990 

0.991 

0.991 

1.394 

2.276 

1.279 

221.175 

113.409 

104.478 

0.030 

0.015 

0.014 

2.960 

2.977 

2.067 

0.336 

0.324 

0.327 

0.664 

0.641 

0.663 

2.799 

3.386 

2.581 

1.012 

0.740 

0.800 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

3600 

3600 

3600 

11677983 

11998747 11975558 

240 

300 

360 

0.2526 

0.3700 

0.4536 

0.0025 

0.0063 

0.0173 

0.0268 

0.0332 

0.0371 

0.0208 

0.0223 

0.0240 

0.0026 

0.0067 

0.0182 

0.735 

0.750 

0.749 

0.378 

0.378 

0.378 

1.694 

2.063 

2.801 

11238.1 

9644.1 

8977.4 

84962.0 

92749.9 104238.5 

121798.6 

166203.6 224148.8 

95474.6 

73859.8 

66008.1 

0.985 

0.983 

0.980 

4.322 

1.669 

0.907 

58.252 

91.121 

175.220 

0.008 

0.012 

0.024 

2.981 

1.490 

0.593 

0.333 

0.340 

0.339 

0.630 

0.661 

0.708 

5.352 

3.202 

2.863 

0.605 

0.864 

1.064 
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Table 15. Characteristics of 7200 rpm magnet machines 


nuaber of pole pairs 
power, hp 
efficiency factor 
voluie factor 
iagnet factor 
shaft rp* 
stator current density 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
iagnet radial dimension 
stator slot factor 
rotor slot factor 
envelope voluie 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
iagnet weight 
aagnet voluie 
xs/turns-squared, p.u. 
•agnet overhang 
air gap flux density 
overall length 
overall diaieter 


1 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

7200 

7200 

7200 

12000000 12000000 8393050. 

120 

240 

360 

0.2266 

0.1938 

0.2622 

0.0H6 

0.0059 

0.0095 

0.1045 

0.0423 

0.0392 

0.0285 

0.0284 

0.0315 

0.0153 

0.0062 

0.0099 

0.749 

0.750 

0.750 

0.378 

0.378 

0.378 

1.247 

0.811 

1.113 

8186.9 

5319.9 

6052.8 

30513.8 

38176.1 

65438.7 

21871.7 

54727.9 140716.0 

90777.3 

71500.7 

38515.8 

0.993 

0.992 

0.987 

1.552 

2.384 

1.621 

110.276 

57.403 

87.115 

0.015 

0.008 

0.012 

2.946 

2.884 

1.258 

0.322 

0.331 

0.307 

0.648 

0.635 

0.647 

2.578 

3.207 

2.747 

0.748 

0.541 

0.685 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

25 

25 

25 

7200 

7200 

7200 

12000000 12000000 12000000 

480 

600 

720 

0.2627 

0.2467 

0.2418 

0.0041 

0.0020 

0.0022 

0.0289 

0.0209 

0.0171 

0.0244 

0.0114 

0.0104 

0.0043 

0.0021 

0.0023 

0.750 

0.748 

0.750 

0.378 

0.378 

0.378 

0.931 

1.356 

1.398 

4898.9 

8797.5 

9221.7 

72593.1 

172996.0 218032.2 

208133.7 620002.2 937689.5 

52482.8 

49628.2 

49815.1 

0.983 

0.958 

0.941 

1.542 

3.965 

4.508 

36.486 

42.995 

50.295 

0.005 

0.006 

0.007 

2.246 

1.608 

1.142 

0.328 

0.338 

0.319 

0.650 

0.631 

0.629 

2.627 

4.968 

5.493 

0.640 

0.562 

0.543 
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Table 16. Characteristics of 180 rpm magnet machines 


nuiber of pole pairs 

i 

2 

3 

4 

5 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

1BOOOO 

180000 

180000 

180000 

180000 

voluie factor 

1286.1 

1286.1 

1286.1 

1286.1 

12B6.1 

aagnet factor 

25 

25 

25 

25 

25 

shaft rp* 

180 

180 

180 

180 

180 

stator current density 

7.80E+06 1.18E+07 

1.20E+07 

1.09E+07 

1.20E+07 

synchronous frequency 

3 

6 

9 

12 

15 

rotor radius 

0.4399 

0.1787 

0.9745 

0.7389 

0.8993 

gap diiension 

0.0413 

0.0036 

0.0201 

0.0166 

0.0214 

back iron depth 

0.1994 

0.0377 

0.1491 

0.0779 

0.0865 

stator slot depth 

0.0932 

0.0282 

0.0327 

0.0476 

0.0630 

aagnet radial dimension 

0.0433 

0.0038 

0.0212 

0.0174 

0.0225 

stator slot factor 

0.463 

0.531 

0.750 

0.750 

0.740 

rotor slot factor 

0.378 

0.378 

0.378 

0.378 

0.378 

envelope voluie 

29.375 

36.694 

34.709 

20.521 

22.314 

envelope weight 

211948.9 284972.1 

156312.7 

107071.7 

77934.3 

hysteresis loss 

20311.6 

53970.7 

47089.7 

41032.7 

36926.0 

eddy current loss 

364.0 

1934.3 

2531.5 

2941.1 

3308.5 

stator copper loss 

7.90E+05 2.54E+06 6.97E+05 

7.01E+05 

7.11E+05 

full load efficiency 

0.960 

0.881 

0.963 

0.963 

0.962 

active length 

12.100 

171.774 

3.195 

4.560 

1.866 

aagnet weight 

4072.4 

2079.6 

1941.0 

1048.8 

1087.7 

■agnet voluie 

0.550 

0.281 

0.262 

0.142 

0.147 

xs/turns-squared, p.u. 

2.996 

2.990 

2.981 

2.983 

2.870 

■agnet overhang 

0.003 

0.034 

0.332 

0.010 

0.339 

air gap flux density 

0.619 

0.619 

0.674 

0.618 

0.704 

overall length 

14.198 

172.518 

7.258 

7.651 

5.639 

overall diaieter 

1.548 

0.496 

2.353 

1.762 

2.140 
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Table 17. Characteristics of more 180 rpm magnet machines 


nutber of pole pairs 

6 

7 

8 

9 

10 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

volute factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

lagnet factor 

25 

25 

25 

25 

25 

shaft rpt 

180 

180 

180 

180 

180 

stator current density 

6.03E+06 1.20E+07 

1.20E+07 9.42E+06 

1.20E+07 

synchronous frequency 

18 

21 

24 

27 

30 

rotor radius 

0.6047 

1.2439 

0.9696 

0.9792 

0.9526 

gap diiension 

0.0139 

0.1161 

0.0089 

0.0031 

0.0079 

back iron depth 

0.0426 

0.1117 

0.0546 

0.0469 

0.0430 

stator slot depth 

0.0337 

0.0627 

0.0312 

0.0255 

0.0360 

■agnet radial ditension 

0.0146 

0.1219 

0.0093 

0.0033 

0.0083 

stator slot factor 

0.536 

0.750 

0.748 

0.575 

0.748 

rotor slot factor 

0.378 

0.378 

0.378 

0.378 

0.378 

envelope volute 

44.578 

53.524 

29.253 

43.668 

26.292 

envelope weight 

313139.6 

98601.8 

120336.5 234296.3 

102286.2 

hysteresis loss 

186214.5 

67518.3 

96060.7 214344.0 

101073.4 

eddy current loss 

20021.3 

8469.3 

13770.9 

34568.5 

18111.9 

stator copper loss 

3.91E+05 7.15E+05 

5.34E+05 3.91E+05 

5.31E+05 

full load efficiency 

0.970 

0.960 

0.968 

0.968 

0.967 

active length 

24.182 

0.651 

3.523 

7.417 

3.167 

tagnet weight 

3800.1 

3431.9 

932.5 

704.3 

714.1 

tagnet volute 

0.514 

0.464 

0.126 

0.095 

0.097 

xs/turns-squared, p.u. 

0.546 

0.525 

2.350 

2.990 

2.391 

tagnet overhang 

0.010 

0.329 

0.336 

0.340 

0.325 

air gap flux density 

0.619 

0.859 

0.669 

0.644 

0.671 

overall length 

26.715 

6.579 

7.474 

11.360 

7.042 

overall diateter 

1.390 

3.069 

2.128 

2.109 

2.079 
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$150,000, rendering the 180 rpm machines less economic to 
build. A machine such as these may cost between four and 
eight million dollars. Only those machines with magnet 
costs that are reasonable with respect to the other material 
costs should be candidates for design. 

In all the permanent magnet machines, stator current 
density went to the maximum. The overall length and overall 
diameter of the more reasonable machines would allow them to 
fit in machinery spaces aboard a ship. 

4.5. Verification 

No high-power permanent magnet machines were discovered 
during the search to find a benchmark. Because of the high 
material cost and the competition afforded by synchronous 
and induction machines, it seems none have been built. 
Several paper studies were found [13, 19, 20, 21, and 22], 
and the parameters resulting from this computer modeling 
seem to agree with them. The machine size is what was ex¬ 
pected, given the lower air gap flux density. The ef¬ 
ficiency was higher than the synchronous and induction 
machines. All-in-all, this modeling gave good machines. 
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Chapter Five. Induction machines 


The stator of an induction machine is the same as those 
of synchronous and permanent magnet machines. The rotor is 
significantly different. There is no independent mechanism 
to produce a rotor flux wave. The rotor winding is shorted, 
whether it is wound or cast, so that as the stator flux wave 
passes over the rotor, currents are induced in the winding. 
These currents produce only a small reaction flux, but it 
still tends to align with the stator flux wave. When at 
operating speeds, the rotor speed is a bit slower than the 
stator flux wave speed, and the difference in speeds is 
called slip. Typically, slip is a few percent of the stator 
frequency. The rotor currents are at slip frequency. If 
the rotor and stator speeds were the same, slip would be 
zero, there would be no tendency to align and torque would 
be zero. Then, the rotor would lag behind the stator until 
current was induced in the rotor winding by the passing 
stator flux wave and torque was again produced. 

If solid bars are used as the rotor winding, they are 
shorted at the ends of the rotor by end rings, to form what 
is called a "squirrel cage" rotor. If actual turns are 
used, the winding may be shorted through external resis¬ 
tances to affect the starting and torque-slip characteris¬ 
tics of the machine. Fitzgerald et al [23] and Alger [24] 
discuss induction machine characteristics in some detail. 

5.1. Assumptions 

A squirrel cage rotor was assumed for these machines. 
Copper was designated as the material for the rotor bars. 
However, these machines will be fed from a frequency 
changer, so only one layer of bars was used and the effects 
of magnetic diffusion ignored in the analysis (see 
Appendix D for a derivation of the components of an induc¬ 
tion machine equivalent circuit). The number of rotor bars 
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was arbitrarily set at 71. This quantity should not cause 
undesirable harmonics, as it will not be an integral mul¬ 
tiple of the number of poles or stator slots in any machine. 
The number and width of the rotor bars were inextricably 
entwined and could not be separated in the analysis. 

5.2. Machine description 

An attempt was made to design the same rpm and pole- 
pair machines as was performed for the synchronous and per¬ 
manent magnet models, but problems in limiting rotor current 
density allowed only a few of the machines to be designed. 
For example, no 7200 rpm machines were designed and 180 rpm 
machines could only be designed with up to twelve poles. 

All of the induction machines are listed in the tables 
starting on page . Only medium confidence should be placed 
in the induction machine designs, as there were some conver¬ 
gence difficulties in slip. (Slip is not listed in the 
tables for that reason.) 

5.2.1. Efficiency 

The higher rpm machines showed a slight increase in ef¬ 
ficiency as rpm increased. There was much movement around 
the average value of 97.5%. The movement decreased as rpm 
increased. With only six machines, it is hard to detect a 
trend in 180 rpm machine efficiency. Apparently, efficiency 
did increase with the number of pole-pairs, with all ef¬ 
ficiencies below 90%. Developmental studies for this thesis 
showed that off-design-point efficiencies for the 180 rpm 
machines were sometimes below 70% for the endurance speed 
condition. 
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EFFICIENCY 


Figure 14. High rpm induction machine efficiency 


INDUCTION MOTOR EFFICIENCY 


25775 HP 



NUMBER OF POLE-PAIRS 
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EFFICIENCY 


Figure 15. 180 rpm induction machine efficiency 


INDUCTION MOTOR EFFICIENCY 


DIRECT-DRIVE, 180 RPM, 25775 HP 



NUMBER OF POLE PAIRS 
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VOLUME (M~3) 


Figure 16. High rpm induction machine volume 


INDUCTION MOTOR VOLUME 


25775 HP 



NUMBER OF POLE-PAIRS 
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VOLUME (M~3) 


Figure 17. 180 rpm induction machine volume 


INDUCTION MOTOR VOLUME 


DIRECT—DRIVE, 180 RPM, 25775 HP 



NUMBER OF POLE PAIRS 
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WEIGHT (KG) 
(Thousands) 


Figure 18. High rpm induction machine weight 


INDUCTION MOTOR WEIGHT 


25775 HP 



NUMBER OF POLE-PAIRS 


75 
































































WEIGHT (KG) 
(Thousand s) 


Figure 19. 180 rpm induction machine weight 


INDUCTION MOTOR WEIGHT 


DIRECT—DRIVE, 180 RPM, 25775 HP 



NUMBER OF POLE PAIRS 
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5.2.2. Weight and volume 


The best of the high rpm machines rivaled the 
synchronous machines in weight and volume. The worst were 
very bad. Weight and volume decreased with an increase in 
the number of poles, but not necessarily with the increases 
in rpm. For the 180 rpm machines, both weight and volume 
decreased dramatically as the number of poles went from two 
to six, with much lower decreases after that. The 180 rpm 
machines were uncompetitive in the synthesis process. 

5.3. Machine characteristics 

The previously mentioned rotor current density dif¬ 
ficulty showed in the rotor slot factor, which was at the 
limit of 0.75 for almost every motor. The stator slot fac¬ 
tor gradually grew with the increase in poles, arriving at 
0.75. The length and diameter of both the 180 and higher 
rpm machines is such that they would fit in machinery 
spaces. 

5.4. Verification 

Induction machines were expected to be close to 
synchronous machines in volume, weight, and efficiency. 

They were, and this comparison served as the verification 
for the induction machine model. Because the confidence 
level in the designed machines is only medium, more work 
would be needed to verify that these machines would have the 
advertised properties if built. 
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Table 18. Characteristics of 1800 rpm induction machines 


nufiber of pole pairs 
power, hp 
efficiency factor 
volute factor 
shaft rp« 
priaary a«p-turns 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope volute 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
rotor copper loss 
taxitua torque 
terminal voits/turn 
air gap volts/turn 
Rl/turns-squared 
Xl/turns-squared 
Xi/turns-squared 
X2/turns-squared 
R2/turns-squared 
overall length 
overall diaaeter 


1 

2 

3 

25775 

25775 

25775 

1B0000 

180000 

180000 

1286.1 

1286.1 

1286.1 

1800 

1B00 

1800 

6.10E+05 

7.74E+05 

2.44E+05 

30 

60 

90 

0.9236 

1.0197 

0.4638 

0.0262 

0.0039 

0.0036 

0.1894 

0.0173 

0.1073 

2.1467 

0.5423 

0.3375 

0.2051 

0.3399 

0.1546 

0.285 

0.641 

0.355 

0.750 

0.750 

0.750 

148.863 

36.240 

7.558 

376603.7 

80669.2 

29447.0 

IB0951.2 

19122.1 

6883B.7 

32425.7 

6853.2 

37006.8 

1.03E+06 

1.45E+06 

498668.3 

0.939 

0.928 

0.969 

0.191 

0.090 

0.759 

8427.7 

8543.0 

19627.2 

185470 

185470 

185470 

234.31 

45.07 

266.45 

63.29 

32.92 

126.5B 

9.23E-07 

B.09E-07 

2.B0E-06 

3.70E-04 

3.89E-05 

9.60E-04 

2.51E-03 

4.41E-03 

1.23E-02 

2.03E-05 

3.66E-05 

3.76E-04 

4.50E-0B 

9.97E-09 

3.89E-07 

3.990 

4.184 

2.629 

6.572 

3.166 

1.824 


25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

1B00 

1800 

1800 

1.91E+05 

1.B6E+05 

1.79E+05 

120 

150 

180 

0.4352 

0.4060 

0.3945 

0.0029 

0.0049 

0.0051 

0.0093 

0.0405 

0.0460 

0.3490 

0.1902 

0.138B 

0.1451 

0.1353 

0.1176 

0.640 

0.699 

0.749 

0.750 

0.750 

0.748 

5.880 

3.760 

3.164 

21499.7 

15333.1 

13395.9 

39013.4 

42456.7 

48214.8 

27964.1 

38040.3 

51839.3 

165391.9 

264414.2 

320360.8 

0.987 

0.981 

0.977 

0.930 

1.000 

1.081 

22347.5 

30471.4 

37148.9 

185470 

185470 

185470 

233.56 

193.76 

189.53 

145.45 

145.98 

153.34 

1.50E-06 

2.54E-06 

3.32E-06 

9.53E-04 

6.B2E-04 

6.17E-04 

1.28E-02 

6.23E-03 

5.15E-03 

6.52E-04 

7.79E-04 

9.08E-04 

5.37E-07 

6.62E-07 

B.49E-07 

2.682 

2.644 

2.6B0 

1.593 

1.283 

1.169 
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Table 19. Characteristics of 2400 rpm induction machines 


nuaber of pole pairs 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

1B0000 

180000 

180000 

voluie factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

2400 

2400 

2400 

2400 

2400 

pritary aap-turns 

5.94E+05 

3.78E+05 

2.45E+05 

1.78E+05 

1.93E+05 

synchronous frequency 

80 

120 

160 

200 

240 

rotor radius 

0.7541 

0.7958 

0.4803 

0.4126 

0.3841 

gap ditension 

0.0069 

0.0020 

0.0094 

0.0087 

0.0081 

back iron depth 

0.0952 

0.0013 

0.0834 

0.0575 

0.0448 

stator slot depth 

0.6169 

0.5475 

0.2755 

0.2341 

0.1561 

rotor slot depth 

0.2514 

0.1979 

0.1601 

0.1375 

0.1197 

stator slot factor 

0.454 

0.477 

0.417 

0.413 

0.750 

rotor slot factor 

0.750 

0.750 

0.639 

0.568 

0.750 

envelope volume 

23.889 

21.180 

6.201 

4.554 

2.942 

envelope weight 

60526.1 

30880.0 

19931.3 

17243.1 

11459.0 

hysteresis loss 

78227.5 

26111.4 

80426.4 

86278.8 

54014.3 

eddy current loss 

37381.4 

18716.2 

76864.2 

103071.8 

77433.0 

stator copper loss 

1.10E+06 

334484.0 

376335.4 

306069.2 

282567.9 

full load efficiency 

0.940 

0.980 

0.972 

0.974 

0.978 

active length 

0.142 

0.189 

0.533 

0.909 

0.852 

rotor copper loss 

8819.4 

7426.0 

17135.5 

21589.8 

27954.3 

aaxiftui torque 

139102.5 

139102.5 

139102.5 

139102.5 

139102.5 

terminal volts/turn 

103.23 

113.69 

229.60 

346.92 

212.59 

air gap volts/turn 

51.27 

72.14 

122.69 

179.75 

156.78 

Rl/turns-squared 

1.04E-06 

7.82E-07 

2.08E-06 

3.23E-06 

2.52E-06 

Xl/turns-squared 

1.50E-04 

2.32E-04 

7.89E-04 

1.67E-03 

7.40E-04 

Xa/turns-squared 

3.88E-03 

1.25E-02 

3.37E-03 

4.31E-03 

3.35E-03 

X2/turns-squared 

6.22E-05 

1.56E-04 

4.63E-04 

1.07E-03 

9.36E-04 

R2/turns-squared 

2.84E-08 

4.44E-08 

2.97E-07 

7.69E-07 

6.73E-07 

overall length 

3.186 

3.380 

2.492 

2.594 

2.420 

overall diaaeter 

2.946 

2.693 

1.697 

1.426 

1.186 
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Table 20. Characteristics of 3000 rpm induction machines 


nuaber of pole pairs 

2 

3 

4 

5 

6 

power, hp 

25775 

25775 

25775 

25775 

25775 

efficiency factor 

180000 

180000 

180000 

180000 

180000 

voluae factor 

1286.1 

1286.1 

1286.1 

1286.1 

1286.1 

shaft rpa 

3000 

3000 

3000 

3000 

3000 

pnaary aap-turns 

3.10E+05 

3.62E+05 

2.80E+05 

1.76E+05 

1.37E+05 

synchronous frequency 

100 

150 

200 

250 

300 

rotor radius 

0.5198 

0.6366 

0.5580 

0.3660 

0.3325 

gap diiension 

0.0138 

0.0036 

0.0070 

0.0070 

0.0039 

back iron depth 

0.1739 

0.0134 

0.0178 

0.0503 

0.0388 

stator slot depth 

0.8926 

0.3785 

0.3668 

0.2274 

0.1288 

rotor slot depth 

0.1733 

0.2122 

0.1860 

0.1220 

0.1093 

stator slot factor 

0.340 

0.607 

0.545 

0.476 

0.721 

rotor slot factor 

0.750 

0.750 

0.750 

0.750 

0.750 

envelope voluae 

22.412 

10.122 

8.000 

3.379 

2.065 

envelope weight 

72906.0 

19541.1 

15961.5 

13000.9 

9016.0 

hysteresis loss 

185978.6 

25156.1 

44235.6 

74610.2 

52848.0 

eddy current loss 

111088.5 

22539.3 

52845.5 

111415.1 

94701.2 

stator copper loss 

322718.0 

359218.8 

221369.1 

269959.2 

226785.8 

full load efficiency 

0.968 

0.979 

0.983 

0.976 

0.980 

active length 

0.383 

0.189 

0.307 

0.818 

1.006 

rotor copper loss 

8574.1 

9089.1 

9006.9 

18292.9 

25568.1 

aaxiaua torque 

111282 

111282 

111282 

111282 

111282 

terminal volts/turn 

421.67 

101.44 

174.86 

361.45 

251.07 

air gap volts/turn 

119.20 

72.02 

102.80 

179.30 

200.52 

Rl/turns-squared 

1.12E-06 

9.14E-07 

9.39E-07 

2.91E-06 

4.01E-06 

Xl/turns-squared 

1.31E-03 

1.96E-04 

5.04E-04 

1.78E-03 

1.09E-03 

Xa/turns-squared 

4.47E-03 

6.93E-03 

3.82E-03 

5.33E-03 

9.00E-03 

X2/turns-squared 

1.82E-04 

1.72E-04 

3.25E-04 

9.97E-04 

1.54E-03 

R2/turns-squared 

1.60E-07 

5.15E-08 

1.08E-07 

6.66E-07 

1.01E-06 

overall length 

2.517 

2.750 

2.567 

2.310 

2.352 

overall diaaeter 

3.210 

2.064 

1.899 

1.301 

1.008 
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Table 21. Characteristics of 3600 rpm induction machines 


nusber of pole pairs 
power, hp 
efficiency factor 
volume factor 
shaft rpa 
priiary aap-turns 
synchronous frequency 
rotor radius 
gap diiension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
rotor copper loss 
saxiauft torque 
terainal volts/turn 
air gap volts/turn 
Rl/turns-squared 
Xl/turns-squared 
Xa/turns-squared 
X2/turns-squared 
R2/turns-squared 
overall length 
overall diaeeter 


l 

2 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

3600 

3600 

3600 

2.72E+05 

1.70E+05 

2.99E+05 

60 

120 

180 

0.4964 

0.4350 

0.5150 

0.0061 

0.0027 

0.0033 

0.3428 

0.0133 

0.0238 

1.3891 

1.0052 

0.3606 

0.1654 

0.1308 

0.1717 

0.290 

0.338 

0.585 

0.750 

0.750 

0.750 

41.461 

17.960 

6.532 

164461.5 

54127.3 

14922.3 

261300.3 

145818.0 

34381.6 

93647.7 

104519.7 

36966.2 

341655.1 

107764.3 

288857.8 

0.965 

0.981 

0.981 

0.393 

0.700 

0.247 

7562.3 

7220.5 

10787.1 

92735 

92735 

92735 

475.12 

702.61 

143.82 

140.25 

218.89 

91.39 

1.54E-06 

1.25E-06 

1.07E-06 

1.67E-03 

3.93E-03 

3.70E-04 

2.39E-02 

4.19E-02 

9.66E-03 

1.20E-04 

4.65E-04 

2.54E-04 

1.98E-07 

4.57E-07 

1.03E-07 

2.403 

2.451 

2.320 

4.469 

2.913 

1.805 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1236.1 

1286.1 

3600 

3600 

3600 

2.95E+05 

2.05E+05 

1.23E+05 

240 

300 

360 

0.5305 

0.4338 

0.3144 

0.0088 

0.0079 

0.0030 

0.0238 

0.0258 

0.0367 

0.3236 

0.2792 

0.1181 

0.1768 

0.1446 

0.1045 

0.688 

0.553 

0.750 

0.750 

0.750 

0.749 

6.536 

4.331 

1.722 

14074.1 

11549.4 

7593.8 

42814.0 

62801.7 

52473.1 

61376.6 

112537.8 

112835.4 

210564.6 

173450.1 

192320.6 

0.983 

0.982 

0.981 

0.249 

0.480 

0.964 

9954.8 

11162.4 

23436.6 

92735 

92735 

92735 

142.32 

275.51 

261.98 

94.86 

149.83 

218.01 

8.0BE-07 

1.37E-06 

4.25E-06 

3.59E-04 

1.13E-03 

1.18E-03 

2.80E-03 

3.94E-03 

1.26E-02 

3.05E-04 

7.24E-04 

1.84E-03 

9.68E-0B 

2.79E-07 

1.07E-06 

2.406 

2.247 

2.234 

1.773 

1.494 

0.944 
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Table 22. Characteristics of 180 rpm induction machines 


rubber of pole pairs 
power, hp 
efficiency factor 
voluae factor 
shaft rpa 
priaary aaps-turns 
synchronous frequency 
rotor radius 
gap diaension 
back iron depth 
stator slot depth 
rotor slot depth 
stator slot factor 
rotor slot factor 
envelope voluae 
envelope weight 
hysteresis loss 
eddy current loss 
stator copper loss 
full load efficiency 
active length 
rotor copper loss 
aaxiaua torque 
terainal volts/turn 
air gap voits/turn 
Rl/turns-squared 
Xl/turns-squared 
Xa/turns-squared 
X2/turns-squared 
R2/turns-squared 
overall length 
overall diaaeter 


1 

n 

L 

3 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1286.1 

1286.1 

180 

180 

180 

1054077 

564571.5 

569333.6 

3 

6 

9 

0.9956 

0.7236 

0.6896 

0.0025 

0.0060 

0.0195 

0.0227 

0.0240 

0.1568 

0.6714 

0.7116 

0.2536 

0.3319 

0.2412 

0.2299 

0.590 

0.563 

0.734 

0.750 

0.750 

0.750 

47.668 

36.B37 

20.316 

217524.1 

168344.2 

98690.0 

4551.5 

16194.6 

19356.0 

81.6 

580.4 

1040.6 

4.44E+06 

1.96E+06 

3.05E+06 

0.808 

0.901 

0.854 

0.825 

2.047 

1.854 

1.08E+05 

1.29E+05 

2.08E+05 

1854700 

1854700 

1854700 

61.849 

104.401 

58.731 

29.510 

53.243 

45.967 

1.33E-06 

2.05E-06 

3.14E-06 

5.08E-05 

1.58E-04 

6.01E-05 

1.21E-02 

4.61E-03 

8.17E-04 

1.89E-05 

6.91E-05 

7.91E-05 

1.16E-07 

4.44E-07 

4.33E-07 

4.817 

4.965 

4.691 

3.384 

2.930 

2.239 


4 

5 

6 

25775 

25775 

25775 

180000 

180000 

180000 

1286.1 

1236.1 

1286.1 

180 

180 

180 

510008.9 

512534.4 

324264.5 

12 

15 

18 

0.6616 

0.7411 

0.6962 

0.0334 

0.0235 

0.0054 

0.1158 

0.1038 

0.0812 

0.2231 

0.2200 

0.1466 

0.2039 

0.1845 

0.1501 

0.750 

0.694 

0.718 

0.683 

0.696 

0.701 

19.450 

21.357 

17.600 

94554.9 

93084.9 

87558.9 

24294.6 

31609.4 

36548.6 

1741.4 

2832.1 

3929.6 

3.15E+06 

2.59E+06 

2.09E+06 

0.850 

0.872 

0.891 

2.486 

2.159 

3.090 

2.23E+05 

1.94E+05 

2.14E+05 

1854700 

1854700 

1854700 

76.530 

76.097 

87.418 

59.124 

57.510 

77.345 

4.04E-06 

3.28E-06 

6.61E-06 

9.01E-05 

9.33E-05 

1.12E-04 

*r 

<=> 

i 

LU 

<=> 

-o 

5.08E-04 

2.49E-03 

1.37E-04 

1.27E-04 

2.40E-04 

7.40E-07 

6.20E-07 

1.15E-06 

5.266 

5.217 

5.897 

2.068 

2.177 

1.859 
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Chapter Six. Nominal ship design 


6.1. Technology sensitivity analysis 

Technology sensitivity analyses, such as this thesis, 
must be able to quantitatively compare similar technologies. 
The pertinent differences must be made apparent through ap¬ 
propriate analysis. Inherent in the analysis must be the 
consideration of the global system complexity. Naval ships 
are extremely complex and the effects of various tech¬ 
nologies can be lost in the complexity. One methodology for 
technology characterization in naval ships has been proposed 
by Goddard [25]. This method has been followed to show the 
benefits of electric drive. 

6.2. ASSET 

The Advanced Surface Ship Evaluation Tool (ASSET) was 
developed over several years to be the U. S. Navy’s premier 
ship design computer program. It has its roots in HANDE, a 
hydrofoil design program developed in the 1970’s. The ship 
design spiral is traversed in an iterative fashion until 
convergence on a number of parameters is achieved. Boeing 
Computer Services is the contractor for ASSET, under the su¬ 
pervision of the David Taylor Naval Ship Research and 
Development Center, at Carderock, Maryland. i6 it is divided 
into large modules by ship type. These modules include 
Monohulls, Hydrofoils, and SWATH ships. 

The geometry of a particular ship is input to ASSET. 

The following were used as characteristic ship traits: full 
load displacement, certain Ship Work Breakdown Structure 
(SWBS) weight groups, 1? endurance fuel load at 20 knots, 

16. Greenwood, R.W., and Fuller, A.L., "Development of a 
Common Tool for Ship Design and Technology Evaluation," 

Proc. SNAME New England Section Marine Computers 1986. 

17. See Appendix E for more information on SWBS. 
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draft, maximum and sustained speed for 51,550 installed 
horsepower, and transverse GM. From the changes in these 
parameters during the various computer runs, the effects of 
transmission choices were noted. 

The inputs to ASSET describe the ship that is being 
designed. The outputs cover the range of calculations pos¬ 
sible in structures, volume, space, machinery, propeller 
characteristics, resistance, powering, and weight. There 
are ASSET performance modules on cost, stability, hydros¬ 
tatics, seakeeping, manning and space but the usual syn¬ 
thesis output is of more use during a technology sensitivity 
analysis. 

The descriptions of several ships are contained in an 
ASSET data base. For a particular ship, a Current Model is 
maintained that holds all of the parameters to describe that 
particular ship. In ASSET Version 2.0, over 380 parameters 
are used for each ship description. User control over most 
of these parameters is possible, or control may be given to 
the executive program which will then "design" a ship sub¬ 
ject to whatever constraints the user desires. 

Some intricacies of ship design are not handled well by 
ASSET. For example, the program is not able to handle 
equipment re-arrangements easily, and almost all equipment- 
level volumes are approximated from studies of past ships.i 8 
For this reason, some equipment-level weights and volumes 
need to be calculated off-line and input to the program 
through its weight adjustment facility, especially if ac¬ 
curacy in these areas is important to the study being per¬ 
formed. 

The baseline ship used in these studies has a full load 
displacement of 5485 LT, carries 272 crew members, is 425 
feet long, and has a primary mission of anti-submarine war- 

18. The Enhanced Machinery Module [27], in the process of 
being made available, will improve the situation dramati¬ 
cally. Some of the relationships from that module were used 
in calculating electric propulsion weights. 
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fare (ASW). It is armed for that purpose and has equipment 
in keeping with its size. The baseline ship is described in 
more detail below. 

6.2.1. Margins 

A naval ship design has margins in weight, vertical 
center of gravity (KG), space, ship service electrical gen¬ 
eration, propulsion power, accommodations, and structural 
strength which allow for equipment, mission, and system 
growth over its projected thirty year life. Without these 
margins, the ship would be difficult to modify, because 
whatever might be added in these areas would have to be paid 
for by a removal. For example, if a 50 ton radar system 
were added, the original 40 ton radar system and 10 tons of 
fuel might be removed to leave the ship at its original 
weight. With margins, the 10 tons of fuel might not be 
removed. 

Margins are typically split into Acquisition and Serv¬ 
ice Life allowances. Acquisition margins recognize the fact 
that ship specifications change over the design cycle and 
during construction. For example, the fourth ship built may 
have a different weapons system than the first, with a dif¬ 
ferent electrical requirement. If the electrical generation 
plant had to be changed during construction to accommodate 
the new weapons system, the total cost might be prohibitive. 
If an Acquisition margin is built into the original design, 
this may not occur. A Service Life margin makes allowance 
for configuration changes over the life of the ship. 

The ASSET program uses margins when synthesizing a 
ship. The margins are under operator control. The margins 
suggested by Goddard and used in this analysis are listed in 
the table.1® 

19. Goddard, C.H., "A Methodology for Technology Charac¬ 
terization and Evaluation for Naval Ships," S.M. NAME AND 
O.E. thesis, MIT, 1985, p. 31. 
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Table 23. Recommended technology assessment design margins 

for a monohull surface combatant 


Acquisition Service Life 


Weight 8 

KG 

Space 

Electrical^ 
Propulsion power 0 

Accommodations 

Strength 


12.5% of SWBS 1-7 10% 

12.5% of KG 1-7 1.0 ft. 

0 (no excess volume) 0 

20 % 20 % 


10% total EHP prior to prelim body plan 
8% prior to self-propelled model tests 
Accom = 1.1 x ship manning at delivery 
2.24 KSI of marginal stress at delivery 
(Max primary stress for hull material) 


Notes: 

a. The service life weight margin applies only to naval 
architectural limits of the ship (reserve buoyancy, 
stability, structures), not to the final design weight. 

b. In sizing the electric plant, the calculated maximum 
electric load plus these design margins shall be met with 
one generator out of service. The remaining generators 
shall not be loaded in excess of 90%. Note that the service 
life margin is not applied to SWBS group 200 which would be 
expected to remain stable over the life of the ship. 

c. Performance requirements (sustained speed, endurance 
range) are met at delivery full load displacement. 


6.3. Philosophy of effort 

The nature of this technology characterization required 
that certain limits be imposed on the total effort. (If the 
Naval Sea Systems Command were to do this study, many people 
would simultaneously be employed to investigate every 
detail.) Some items were fixed, some were allowed to float 
with the design. 

The hardest item to handle is volume. There are very 
few ways of adjusting volume as easily as weight is ad¬ 
justed. One way is through the use of Marginal Volume Fac¬ 
tors, which equate a weight penalty with every increase in 
volume. (See Howell [26].) 

The differences among transmission systems appear 
primarily in the machinery spaces and fuel tanks of a ship. 
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ASSET can handle the tankage changes without aid. However, 
it does no equipment arrangement or space analysis inside 
the machinery spaces, leaving that for arrangement experts 
to do off-line. It was decided to keep the machinery space 
volume the same in all ships of this study, no matter how 
that volume was divided into spaces. Changes in equipment 
volume will be noted in the analysis and left for the advan¬ 
tage of others in machinery space rearrangement. 

Moderate to high technical risk has been accepted in 
specifying equipment cooling. Current densities at the 
limits of cooling technology are used, assuming that liquid 
cooling of both the stator and rotor can be performed. 

Lower current densities would be required if such cooling 
were not possible, resulting in slightly larger, heavier, 
less efficient machines. This last statement was proven 
during the course of the thesis research, as the first 
(chronologically) current densities used were two-thirds 
larger than those listed herein. 

Other risk areas include the use of advanced vacuum 
switchgear and the assumed efficiencies of reduction gears 
and power converters. These are low risk items; the 
technology is well understood and commercially available. 


Table 24. Ship design items held constant during analysis 


Endurance speed 
Endurance range* 

Machinery box volume 
Installed horsepower 

Payload weight, volume, and electrical 

Length 

Beam 

Ship electrical load 
Ship molded lines 
Manning 

Deckhouse and superstructure geometry 


20 knots 

5500 nautical miles 
109,670 ft3 
51,550 hp 
requirements 
425 ft. 

55 ft. 

2030 kW (24 hr. avg) 
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Table 25. Ship design items allowed to float during analysis 


Maximum and sustained speeds 
Endurance fuel load 
Ship full load displacement 
Ship draft 

Ship resistance and powering 

Ship arrangement outside machinery spaces 


Note: 

a. Endurance range was allowed to float as a comparison 
between two electrical transmission ships and the mechanical 
transmission baseline ship. It was held constant for the 
rearranged ship. 
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Most of the ship synthesis has been left to ASSET 
though parameters from Goddard were used where possible. 
ASSET designs a reasonable, generic ship with good seakeep¬ 
ing characteristics. As shown above, some ship characteris¬ 
tics were frozen to ensure transmission comparisons were not 
performed with different ships. Leaving the ship synthesis 
to ASSET allowed concentration on the specifics of the 
propulsion plant. 

6.4. Baseline ship 

The baseline ship has a mechanical transmission, i.e., 
two power trains consisting of a gas turbine, clutch and 
coupling, reduction gear, shafting, and propeller. There 
are two machinery rooms, each containing one gas turbine. 

The gas turbine used as the model is the General Electric 
LM-2500, rated at 25,775 brake horsepower. This is a very 
common marine gas turbine. It and its predecessors are 
powering the latest classes of naval combatant, such as the 
DD963, FFG7, and DDG51. 

The locked-train double-reduction gears are reversible, 
allowing the use of fixed-pitch propellers. There is no 
mechanical cross-connect allowed between the shafts. Except 
for the power level, this is the gear system being employed 
on the DDG51 class. Gears of this sort are about one per¬ 
cent inefficientzo per reduction stage. Since these are 
double-reduction gears, an efficiency of 98% was used. 

An endurance speed of 20 knots has been specified. 

This is in keeping with standard fleet practices. An en¬ 
durance range of 5500 nautical miles permits ocean crossings 
without refueling. The lack of a cross-connect capability 
between the two shafts means at least two gas turbines will 


20. Inefficiency = 1 - efficiency. Information on stage in¬ 
efficiency is from a conversation with Mr. Samuel Shank, the 
author of the ASSET Enhanced Machinery Module [27]. 
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be on-line during endurance cruising. This is inherently 
inefficient; some operators alleviate the inefficiency by 
declutching one shaft and free-wheeling that propeller, 
reaching desired speeds by loading up the operating gas tur¬ 
bine . 

A measure of initial static stability is the ratio 
GM/Beam. GM is the vertical distance between the center of 
gravity and the metacenter of the ship. Typical values for 
this ratio are 8-10%. A lower value (6.5%) has been ac¬ 
cepted for purposes of comparison with the variants. A 
large ship redesign effort would have been necessary to 
bring GM/B into a better range. 

The electrical generation plant consists of three gas 
turbines, each driving a 2000 kW generator. The data used 
for the gas turbines was taken from the Detroit Allison 501K 
turbine-generator set used aboard the DD963 ships. 

Both the deckhouse and main hull are constructed of 
High Tensile Strength (yield stress = 50,000 psi) steel. 
Active stabilizing fins and a sonar dome are included in the 
design. The payload is listed in Table 26. A coarse layout 
of the machinery spaces is shown in Figure 20. 


Table 26. Payload for baseline and variant ships 

FFG7 Command and control suite 
Satellite, UHF, and HF communications 
SLQ32V3 electronic countermeasures 
NIXIE acoustic countermeasure 

SPS-49 two-dimensional air search and tracking radar 

SPS-55 surface search and tracking radar 

SQR-19 towed array surveillance system 

MK92 missile and gun fire control system 

Harpoon fire control system 

LAMPS III helicopters and support system 

JP-5 aviation fuel 

MK32 over-the-side torpedo system 

MK13 guided missile launching system 

MK75 76mm gun 

Close-In-Weapon-System 

Small arms 

Appropriate ammunition and reloads 
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Figure 20. Baseline ship subdivision and machinery arrange¬ 


ment 
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6.5. Backfit ship 


If the mechanical transmission of the baseline ship 
were replaced with an electrical transmission without chang¬ 
ing the subdivision of the ship, it would be as if an older 
ship had been updated, or "backfitted,” with new technology. 
This is the idea behind the first variant ship. 

Propulsion motors and generators were added into the 
propulsion plant, in the original spaces. Shafting still 
runs from the machinery spaces to the propellers. A reduc¬ 
tion gear is still necessary for higher rpm propulsion 
motors, but electrical cross-connect may improve the en¬ 
durance fuel efficiency. 

Some rearrangement within the machinery space is neces¬ 
sitated by the backfit. The machinery spaces in the 
baseline ship are not long enough to contain the stack-up 
length of gas turbine, generator, motor, and reduction gear, 
without "folding" the power train. This may be accomplished 
by changing the design of the reduction gear or by placing 
the gas turbine and generator (which require a mechanical 
connection) side-by-side (or transversely) with the motor 
and reduction gear (another mechanical connection), using 
transmission line to electrically connect them. Since there 
are a variety of ways to rearrange the machinery box, and 
since the chosen method has no effect on the analysis, the 
rearrangement was not specified. 

The propulsion generators and motors may operate at 
different rotational speeds, and therefore different elec¬ 
trical frequencies, so power converters must be used between 
them. Power converters change the frequency of the power 
being transferred between the generators and motors, through 
the use of cycloconverters or thyristors. They add another 
inefficiency to the transmission. A reasonable estimate of 
the efficiency of an 18 MW power converter is 97%. 21 

21. Professor John Kassakian, MIT, private communication. 
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6.6. Rearranged ship 


The second variant ship takes advantage of the benefits 
of machinery space rearrangement. The same machinery space 
total volume is preserved but split into five spaces. The 
propulsion motors are placed very near the propellers, 
resulting in much shorter shafting runs and an increased 
GM/B ratio. The decrease in shaft length means a decrease 
in shaft weight and more space available outside the 
machinery box. The shafts previously ran through shaft al¬ 
leys that may be returned to other uses. This rearrangement 
is almost certainly not the optimum one and can be improved 
in the sense of space efficiency. It is an arrangement, 
however, that can demonstrate the benefits to be expected of 
a ship designed for electric drive. Power converters again 
connect the generators and motors. Transmission line forms 
the connections, at a much lower weight than shafting. 

Many other choices in large components are possible for 
this rearranged ship. For example, three propulsion gener¬ 
ators and gas turbines driving two or four propulsion motors 
might have been chosen. The number of prime movers was 
retained from the mechanical baseline ship, however, to make 
the comparison of transmissions realistic. Too many changes 
might have obscured the fundamental differences in 
efficiency, weight and volume. 

Both geared and direct-drive propulsion motors were 
used in this variant. When geared motors are used, the 
reduction gears are also placed near the propellers. A 
coarse layout of the rearranged ship is shown in Figure 21. 

6.7. Weight and volume algorithms 

Few components in this thesis are exact commercial 
models. The weight, volume, and other characteristics are 
taken from those for which data was available. The equa¬ 
tions for shafting and reduction gears were taken from the 
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ASSET theory manuals [28], while the switchgear and braking 
resistor equations came from the ASSET Enhanced Machinery 
Module [27], which is not yet generally available. These 
equations represent much study by ship and equipment desig¬ 
ners, incorporating equipment which .is commercially avail¬ 
able. Where possible, the ASSET equations were verified 
against other studies and actual equipment [29, 30, 31]. 

For example, the machinery in the FFG7 was used as a model 
and verification for reduction gears and shafting. 

Appendix E gives more explanation, as well as presenting a 
computer program used to generate weight and volume figures. 
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Figure 21. Layout of machinery spaces on rearranged ship 
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Chapter Seven. Machine design and system synthesis 


7.1. Machine matrix 

During this study, synchronous, permanent magnet, and 
induction machines were designed at shaft speeds of 180, 
1800, 2400, 3000, 3600, and 7200 rpm. The number of pole- 
pairs for the 180 rpm machines ran from one to twenty-five. 
For the higher rpm machines, from one to six pole-pairs were 
used. If every machine could be a generator or a motor, a 
165x165 machine matrix results. From these 27,225 possible 
combinations, two were chosen and input to the two variant 
ships for synthesis in ASSET. 

These particular rpms were chosen partially because of 
the choice of the gas turbine. The LM-2500 operates at a 
full-load speed of 3600 rpm, making multiples and "nice" 
fractions of that speed desirable. A 3600 rpm, two-pole 
machine has a synchronous frequency of 60 hertz, the stand¬ 
ard in the United States. A four-pole machine at 1800 rpm 
is also a 60 hertz machine. If "nice" frequencies result 
from rpm choices, results may not be obscured. The low rpm 
machine is tied to propeller rpm. For the baseline ship, 
maximum speed propeller rpm is 170 rpm. A ten rpm allowance 
for "battle override" gives a requirement for 180 rpm. 

The reduction in the number of machine combinations is 
a bit more difficult to explain. First, it was observed 
that reduction gears add greatly to the weight and volume of 
the transmission and detract from its efficiency. Second, 
induction generators are notoriously difficult to control. 

It was then decided that generators would not be induction 
machines and any generators used would operate at the same 
shaft speed as the gas turbine, eliminating a possible 
reduction gear. The matrix then measured 12x165 and had 
1980 combinations. 

From this point on, the decision theory espoused in 
Schweppe and Merrill [13] was used, specifically using "knee 
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sets" to eliminate uncompetitive designs. The considera¬ 
tion, or Figure-of-Merit, was always minimum weight, volume, 
and inefficiency. These external characteristics are those 
"seen" by a ship design (since propulsion current and volt¬ 
age were specified in terms of turns -- fertile ground for 
another tradeoff study). Initially, individual motors of 
all of the higher rpms were considered together to select 
the best geared propulsion motors. Since the weight and 
volume of a reduction gear varies with shaft rpm, the weight 
and volume of the reduction gear was included with that of 
the machine to select the best machine-gear system. Direct- 
drive machines were also selected. The 3600 rpm synchronous 
and permanent magnet machines were also considered 
separately as propulsion generators. 

Once the initial selection of machines was made, the 
matrix measured 5x11, or 55 combinations. All of the com¬ 
binations were plotted in knee curves that showed the 
volume, weight, and inefficiency of each transmission. The 
volume and weight of shafting, braking resistors, cooling 
systems, power converters, and the inefficiency of the power 
converters were common to all combinations and were not in¬ 
cluded at this level. The inefficiency of any reduction 
gears was included where appropriate. Three of the gener¬ 
ators were synchronous machines and two were permanent mag¬ 
net machines. Of the motors, two were 1800 rpm synchronous 
machines, two were 1800 rpm permanent magnet machines, three 
were 1800 rpm induction machines, and two were 180 rpm 
synchronous machines. The 1800 rpm motors clearly dominated 
the higher rpm machines, largely because of the differences 
in reduction gear weight and volume. 

Since there was no single dominating combination, ten 
of the 55 combinations were selected. These ten were among 
the best at least twice on the knee curves. These ten com¬ 
binations were composed of only synchronous machines. 

Programs to calculate off-design-point efficiency were 
written. Each motor and generator was evaluated at the 
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power level and rpm appropriate for the sustained and en¬ 
durance speed conditions of the ship, using the delivered 
horsepower (DHP) and propeller rpm taken from the ASSET syn¬ 
thesis run on the baseline mechanical ship. Some of the 
combinations had very low endurance efficiencies. There was 
a sharp division evident between the geared (lower weight, 
much higher volume) and the direct-drive systems. 

The ten combinations, with their maximum, sustained, 
and endurance speed transmission inefficiencies, were again 
made the subject of knee curves. A simple scoring scheme 
was devised to rank the combinations according to their 
grouping on these last knee curves. If a combination was in 
the best group on a particular knee curve, it was given two 
points for that curve. If it was in the second best group, 
it received one point. If it was in neither the best or 
second best group, it received no points. When the scores 
were totaled, two combinations stood out. One was a geared 
drive system and one a direct-drive system. These two com¬ 
binations were used in both of the variant ships and are the 
subject of the next chapter. 

7.2. Knee curves 

Figures 22 through 24 show the knee curves for the 
propulsion generators. The first letter of the generator 
ID indicates whether it is a synchronous machine or a per¬ 
manent magnet machine. All of the generators are 3600 rpm 
machines. The generators selected were SB (four poles), SC 
(six poles), SD (eight poles), PB (four poles), and PC (six 
poles). Generator SA was not selected because of its poor 
showing on the volume-weight curve, even though it was com¬ 
petitive on the volume-efficiency curve. 

The 180 rpm, direct-drive propulsion motor curves are 
in Figures 25 through 27. They were not combined with 
geared motors because one of the points being investigated 
was whether or not geared motors were "better" than direct- 
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drive motors. The clumping of the machines necessitated 
other graphs on different scales to distinguish between the 
machines. Machines 1-25 are synchronous, 26-50 are per¬ 
manent magnet, and 51-56 are induction machines. The number 
within a group indicates the number of poles in the machine, 
e.g., machine 32 has (32-25)x2=14 poles. Valid designs with 
over twelve poles were not achieved for the induction 
machines. Machines 3, 8, 11, and 16 were selected for fur¬ 
ther work. These are all synchronous motors. 

The higher rpm motors were combined with their reduc¬ 
tion gears to form system knee curves. In all of these knee 
curves, machines 1-6 are 1800 rpm, 7-12 are 2400 rpm, 13-18 
are 3000 rpm, 19-24 are 3600 rpm, and 25-30 are 7200 rpm. 
Figure 28 is the volume-efficiency curve for synchronous 
machines, showing the distinct grouping of the machine-gear 
systems due to the high volume of the reduction gears. Note 
the high values of the permanent magnet machines in the 
volume-weight curve Figure 29. The clumping of induction 
motors around the low inefficiencies is shown in the weight- 
efficiency curve of Figure 30. From these curves, the 
motors on page 97 were selected. 

The initial motor and generator combinations were made 
and plotted on more knee curves (Figures 31 to 33). On 
those curves, the high-volume or high-weight nature of the 
combinations can be seen. Since the multiple-attribute 
decision theory embodied in knee curves does not say how to 
select between high-volume or high-weight, the best of each 
were selected. Combinations 1, 2, 8, 9, 12, 13, 19, 20, 30, 
and 31 were chosen. The off-design-point efficiencies were 
calculated and all of the information was plotted. 

Figure 34 is a bar-graph of the maximum, sustained, and en¬ 
durance speed transmission efficiencies of the various com¬ 
binations, including the reduction gears, if any, and power 
converters. The final combination knee curves are sum¬ 
marized in the scoring scheme of Table 27, which was ex¬ 
plained on page 98 . Combinations 12 and 20 were chosen to 
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use in the ships of the study. 

This is a good method to choose among the possible 
machines. During the course of this thesis, the above path 
was followed through several complete iterations and a few 
partial ones. As stated in Schweppe and Merrill, knee 
curves serve very well to eliminate uncompetitive options, 
allowing concentration on the better ones. 

Table 27. Final combination knee curve scores 


Combo ID Firsts Seconds Total 


1 

2 

8 

9 

12 

13 

19 

20 
31 


3 

3 
1 
0 
6 

4 
1 
3 
0 


1 

0 

0 

1 

1 

0 

1 

4 

2 


7 
6 
2 
1 

13 

8 

3 

10 

2 


Conclusion: test combinations 12 and 20 
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Figure 22. Curve of volume-efficiency for 3600 rpm gener¬ 
ators 


PROPULSION GENERATORS 


25775 HP, 3600 RPM 
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Figure 23. Curve of volume-weight for 3600 rpm generators 


PROPULSION GENERATORS 


25775 HP, 5600 RPM 
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Figure 24. Curve of weight-efficiency for 3600 rpm gener¬ 
ators 


PROPULSION GENERATORS 


2S775 HP. 3600 RPM 
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Figure 25. Curve of volume-efficiency for 180 rpm motors 


PROPULSION MOTORS 

DIRECT-DRIVE, 180 RPM, 25775 HP 
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Figure 26. Curve of volume-weight for 180 rpm motors 


PROPULSION MOTORS 


DIRECT-DRIVE, 180 RPM, 25775 HP 
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Figure 27. Curve of weight-efficiency for 180 rpm motors 


PROPULSION MOTORS 


DIRECT-DRIVE, 180 RPM, 25775 HP 
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Figure 28. Curve of volume-efficiency for geared motors 


SYSTEM KNEE CURVE 


SYNCHRONOUS MACHINES. 2S775 HP 
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Figure 29. Curve of volume-weight for geared motors 


SYSTEM KNEE CURVE 
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EFFICIENCY) 


Figure 30. Curve of weight-efficiency for geared motors 


SYSTEM KNEE CURVE 


INDUCTION MACHINES, 25775 HP 



109 



































































TOTAL EFFICIENCY) 


Figure 31. Volume-efficiency curve for initial PM and PG 
combinations 
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Figure 32. Volume-weight curve for initial PM and PG com¬ 
binations 
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Figure 33. Weight-efficiency curve for initial PM and PG 
combinations 
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EFFICIENCY 


Figure 34. Final combination transmission efficiencies 
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Chapter Eight. Analysis 


This chapter presents an analysis of the ships with 
electric transmissions. Standard naval architectural 
methods have been used to observe and comment on the variant 
ships, comparing them to the mechanical baseline ship. Con¬ 
clusions and recommendations follow the analysis. 

The names used to describe the various ships imply 
their internal arrangement, equipment, and ASSET Design Mode 
Indicator (DMI). The two DMIs used were ENDURANCE, when en¬ 
durance range was held constant at 5500 NM, and FUEL WT, 
when the usable fuel weight was held constant at 996.6 LT.22 
The ship names are as follows: 


MECH 23 BASELINE: Baseline, mechanical transmission 
ELEC 23 BASELINE 12: Backfit ship, geared motors 
ELEC 23 BASELINE 20: Backfit ship, direct-drive motors 
NEW MR ELEC 12: Rearranged ship, geared motors 
NEW MR ELEC 20: Rearranged ship, direct-drive motors 
CONSTANT FUEL ELEC 23 BASE 12: Backfit ship, geared 

motors 

CONSTANT FUEL ELEC 23 BASE 20: Backfit ship, direct- 

drive motors 


8.1. Direct effects 

The direct effects of an electric transmission are the 
changes in weight and volume of the propulsion system, as 
well as the transmission efficiency. Included are the 
weight and volume of the propulsion motors and generators, 
transmission lines, cooling systems, switchgear, power con¬ 
verter, exciter, braking resistor, any reduction gears and 
their associated gear oil, and shafting. These items are 
listed in Table 31. A positive difference from the baseline 
ship means a heavier and/or larger ship. 


21. Professor John Kassakian, MIT, private communication. 

22. Not all fuel in a ship is usable. There are nooks and 
corners of fuel tanks that are inaccessible to the fuel sys- 
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Note that the only variant that has a lower direct 
weight effect than the baseline is NEW MR ELEC 12. The ac¬ 
cumulation of weight increases in the others makes them 
heavier, while NEW MR ELEC 12 has lower motor and shafting 
weight than the rest. Geared drive is always lighter than 
direct-drive, largely due to the high weight of the direct- 
drive propulsion motors. With respect to volume, direct- 
drive is always smaller than geared drive, because of the 
reduction gears. All electric transmissions are larger than 
the mechanical baseline, but the smallest variants, within 
motor type groups, are the rearranged ships. 

Table 28 contains the maximum, sustained, and endurance 
speed transmission efficiencies of the two generator-motor 
combinations. Note that the off-design-point efficiency of 
the direct-drive combination is significantly lower than the 
geared combination, even though it does not have the added 
inefficiency of reduction gears. This is in large part due 
to the poor efficiency of the slowly rotating direct-drive 
motor at the endurance speed. 

Table 28. Transmission efficiencies 

Combo Maximum Sustained Endurance 

12 0.9307 0.9266 0.8817 geared combination 

20 0.9209 0.9093 0.7754 direct-drive combination 

Endurance efficiency with one generator driving two motors. 

Table 29. Propulsion generator efficiencies 

PG ID Maximum Sustained Endurance 

SC 0.9891 0.9870 0.9737 geared combination 

SC 0.9891 0.9872 0.9768 direct-drive combination 

Table 30. Propulsion motor efficiencies 
PM ID Maximum Sustained Endurance 

S4 0.9898 0.9876 0.9526 geared combination 

SL8 0.9598 0.9496 0.8184 direct-drive combination 
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The above tables show the efficiencies of the motor and 
generator used in two particular combinations. The 
efficiency of the motor, reduction gears, and power con¬ 
verter have a direct effect on the efficiency of the gener¬ 
ator, as they change the loading point of the generator. 
Generally, motors and generators are more efficient when 
they are loaded more closely to their design point. (The 
same is true of gas turbines.) The inefficiencies of the 
power converters and reduction gears, if any, are included 
in the transmission efficiencies. 
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Table 31. Direct volume and weight effects 


CONSTANT CONSTANT 








FUEL 

FUEL 


MECH 

ELEC 

ELEC 

NEW 

NEW 

ELEC 

ELEC 


23 

23 

23 

MR 

MR 

23 

23 

BASELINE BASELINE 

BASELINE 

ELEC 

ELEC 

BASE 

BASE 



12 

20 

12 

20 

12 

20 

Electric Propulsion weights: 

in LT 






PMs 

0 

18.08 

102.94 

18.08 

102.94 

18.08 

102.94 

P6s 

0 

10.27 

10.27 

10.27 

10.27 

10.27 

10.27 

Trans lines 

0 

0.18 

0.18 

0.77 

0.77 

0.18 

0.18 

Cooling sys 

0 

5.98 

5.98 

5.98 

5.98 

5.98 

5.98 

Switchgear 

0 

1.56 

1.56 

1.56 

1.56 

1.56 

1.56 

Power converter 

0 

7.16 

7.16 

7.16 

7.16 

7.16 

7.16 

Exciters 

0 

3.50 

3.50 

3.50 

3.50 

3.50 

3.50 

Braking resistor 

0 

10.00 

10.00 

10.00 

10.00 

10.00 

10.00 

Reduction gears 

78.90 

41.19 

0.00 

41.19 

0.00 

41.19 

0.00 

Shafting 

69.00 

66.84 

66.84 

41.72 

41.72 

66.B4 

66.84 

W298 (op fluid) 

14.7 

14.7 

0 

14.7 

0 

14.7 

0 

Direct effect 

162.60 

179.46 

208.43 

154.93 

183.90 

179.46 

208.43 

Diff fa Baseline 

0 

16.86 

45.83 

-7.67 

21.30 

16.86 

45.83 

Electric Propulsion volumes: 

in cubic 

feet 





PMs 

0 

108.86 

909.13 

108.86 

909.13 

108.86 

909.13 

PGs 

0 

61.49 

61.49 

61.49 

61.49 

61.49 

61.49 

Trans lines 

0 

20.9 

20.9 

6.39 

6.39 

0.33 

0.33 

Cooling sys 

0 

200 

200 

200 

200 

200 

200 

Switchgear 

0 

70.2 

70.2 

70.2 

70.2 

70.2 

70.2 

Power converter 

0 

1089.73 

1089.73 

1089.73 

1089.73 

1089.73 

1089.73 

Exciters 

0 

245.02 

245.02 

245.02 

245.02 

245.02 

245.02 

Braking resistor 

0 

1422.88 

1422.88 

1422.88 

1422.88 

1422.88 

1422.88 

Reduction gears 

2731.04 

1425.75 

0.00 

1425.75 

0.00 

1425.75 

0.00 

Shafting 

517.71 

501.5 

501.5 

306.06 

306.06 

501.5 

501.5 

Direct effect 

3248.75 

5146.33 

4520.85 

4936.38 

4310.90 

5125.76 

4500.28 

Diff fa Baseline 

0.00 

1897.58 

1272.10 

1687.63 

1062.15 

1877.01 

1251.53 


8.2. Indirect effects 

Indirect effects are again composed of weights and 
volumes, but these are the ripple effects of the propulsion 
system through the ship. For example, if a transmission is 
more efficient at endurance speed, it should be expected 
that less onboard fuel would be needed to achieve the same 
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endurance range as a less efficient transmission. This is 
indeed the case. Another important indirect effect is the 
change in full load displacement, which is tied to the 
powering characteristic of the ship. 23 The following tables 
list the characteristics and indirect effects of the various 
ship configurations. 

Every electric transmission had lower maximum and sus¬ 
tained speeds than the baseline ship (by about 0.43 knots), 
but also lower EHP requirements at those speeds. The lower 
EHPs are a reflection of lower drafts (less resistance). 

The lower speeds show that the transmission efficiencies of 
the variant ships are lower than the mechanical baseline 
ship. There are more components in the electrical power 
trains, hence the lower efficiencies. The speed difference 
of 0.43 knots may be regarded by some as significant; it is 
about the speed reduction to be expected by a fouled bottom. 

The endurance range of all ships except those with con¬ 
stant fuel load is 5500 NM. The fuel load in the others 
varies greatly, showing the benefit of electrical cross- 
connection. In the two constant fuel ships, there was an 
increase in the endurance range of 1350 and 1400 NM, respec¬ 
tively, for the geared and direct-drive transmissions. This 
is an indication of fuel cost savings from the electric 
transmission. If a ship refuels every three steaming days 
(receiving a third of its tank capacity), steams 100 days 
each year, and fuel is priced at about $18 per barrel, this 
represents about a $600,000 savings per ship per year. 

The initial static transverse stability of the variants 
was degraded by the change in propulsion equipment. As pre¬ 
viously stated, the 6.5% GM/Beam ratio of the mechanical 
baseline ship is not as large as desired for an actual ship, 
but provided a benchmark to measure relative changes. 

tem. Typically, 95% of the onboard fuel is usable. 

23. For the same molded lines, ships with higher displace¬ 
ments will have greater wetted surface areas and higher 
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Table 32. General ship characteristics 



NECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEW 

NR 

ELEC 

12 

NEW 

NR 

ELEC 

20 

CONSTANT 

FUEL 

ELEC 

23 

BASE 

12 

CONSTANT 

FUEL 

ELEC 

23 

BASE 

20 

LBP 

♦25 

425 

425 

425 

425 

425 

425 

Ship disp FL 

5485 

5297.5 

5417.1 

5234.7 

5353.4 

5484.1 

5530.2 

Diff fa Baseline 

0 

-187.5 

-67.9 

-250.3 

-131.6 

-0.9 

45.2 

Nav beat 

55 

55.04 

55.02 

55.04 

55.03 

55 

55 

Nav draft 

16.44 

16.07 

16.31 

15.95 

16.18 

16.44 

16.53 

Depth sta 10 

38 

38 

38 

38 

38 

38 

38 

Accoaaodations 

272 

272 

272 

272 

272 

272 

272 

GNt/B 

0.065 

0.059 

0.058 

0.061 

0.064 

0.063 

0.061 

ghi 

950.01 

973.41 

959.1 

979.76 

967.39 

950.07 

943.41 

Full load KG 

23.23 

23.87 

23.75 

23.84 

23.53 

23.34 

23.44 

LCB/LBP 

0.516 

0.514 

0.515 

0.513 

0.514 

0.516 

0.516 

LCF/LBP 

0.565 

0.565 

0.565 

0.564 

0.565 

0.565 

0.565 

Hetted surface 

23347.7 

23010 

23230.5 

22886.3 

23115.7 

23346.6 

23425.8 

Cp 

0.6 

0.596 

0.598 

0.594 

0.597 

0.6 

0.601 

Waterplane area 

17641 

17575.9 

17624.8 

17539.8 

17602 

17640.9 

17653.2 

Hull volute 

618504 

618504 

618504 

618504 

618504 

618504 

618504 

Dkhs volute 

121528 

12152B 

121528 

121528 

121528 

121528 

121528 

Total volute 

740032 

740032 

740032 

740032 

740032 

740032 

740032 

Endur fuel vol 

44365 

36592 

39647 

36367 

39396 

44364 

44364 

rtachy box vol 

109671 

109671 

109671 

109624 

109624 

109671 

109671 
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Table 33. Powering 


CONSTANT CONSTANT 



HECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEW 

HR 

ELEC 

12 

NEM 

HR 

ELEC 

20 

FUEL 

ELEC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

Powering: 

V*ax 

29.13 

29.02 

28.83 

29.09 

28.9 

28.82 

28.7 

EHP 9 Vsax 

34821 

33131 

32741 

33154 

32763 

33066 

32700 

Vsus 

27.94 

27.82 

27.6 

27.88 

27.66 

27.63 

27.48 

EHP 9 Vsus 

28163 

26684 

26165 

26700 

26181 

26636 

26136 

Endurance 

5500 

5500 

5500 

5500 

5500 

6550 

6086 

EHP 9 Vend 

6851 

6632 

6772 

6557 

6698 

6851 

6905 

HPi 

51550 

51550 

51550 

51550 

51550 

51550 

51550 

kWi 

6000 

6000 

6000 

6000 

6000 

6000 

6000 

Avg 24 hr load 

2030 

2030 

2030 

2030 

2030 

2030 

2030 


120 
































































In the constant fuel ships, the high weight of the 
propulsion motors and propulsion motors and propulsion gen¬ 
erators, combined with a smaller amount of vertically-lower 
fuel reduced GM/B by over one-half percent, a not incon¬ 
siderable amount. For the rearranged variants, GM/B 
decreased less than the two previous ships. They also have 
less low fuel, but the propulsion generators are lower than 
in the backfit ships and the propulsion motors are very much 
lower. Only a very small decrease in GM/B was seen in the 
constant fuel ships because the constant fuel load compen¬ 
sated for the increased high weight of the electric trans¬ 
mission. The longitudinal metacentric height, GM1, in¬ 
creased for all variants, though it seems it should have 
decreased with the decrease in waterplane area and draft. 

The machinery space volume was the same for all ships. 

There were no big surprises in the area of weight. The 
structural weight (W100) encloses the same volume in every 
ship, so it was about constant. The propulsion plant weight 
(W200) varied with the type of transmission. Weight groups 
W300, W400, W500, W600, and W700 were virtually identical in 
every ship, and the variable loads were dominated by the 
change in fuel weight. The Design and Builders Margin is a 
function of the light ship weight (summation of W100 through 
W700), so the margin weight moved with the light ship 
weight. The miles-per-gallon figure of NM/LT of fuel shows 
the endurance efficiency of electrical cross-connection. 

Only a few comments need be made regarding Table 35. 

The structural weight fraction shows the changes in full 
load displacement, remembering that the W100 weights were 
all about the same. The same may be said for the weight 
fraction of the W300 through W700 groups and payload weight. 
Higher propulsion plant weights in the variant ships drove 
up the W200 fraction, except for NEW12. The fuel weight 
fraction shows the same behavior as the miles-per-gallon 
figure. 
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Table 34. Ship weights 


CONSTANT CONSTANT 



MECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEH 

HR 

ELEC 

12 

NEH 

HR 

ELEC 

20 

FUEL 

ELEC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

Height: 
DFH onboard 

1049.1 

865.3 

937.5 

860 

931.6 

1049.1 

1049.1 

Usable fuel nt 

996.6 

822 

890.6 

817 

885 

996.6 

996.6 

Diff ft Baseline 

0 

-183.8 

-111.6 

-189.1 

-117.5 

0 

0 

NH per LT fuel 

5.52 

6.69 

6.18 

6.73 

6.21 

6.57 

6.11 

Payload Height 

571.2 

571.2 

571.2 

571.2 

571.2 

571.2 

571.2 

H100 

1684.8 

1686 

1692 

1660.6 

1666.6 

1686.2 

1692.1 

H200 

343.6 

351.8 

387.1 

326.1 

361.3 

351.9 

387.1 

H300 

236.5 

236.5 

236.5 

236.5 

236.5 

236.5 

236.5 

H400 

302.2 

302.2 

302.2 

302.2 

302.2 

302.2 

302.2 

H500 

615.4 

613.5 

614.3 

613.5 

614.2 

615.4 

615.4 

H600 

426.6 

426.6 

426.6 

426.6 

426.6 

426.6 

426.6 

H700 

95.9 

95.9 

95.9 

95.9 

95.9 

95.9 

95.9 

Loads 

1317.1 

1120.9 

1193.2 

1115.6 

1187.2 

1305.1 

1304.9 

D&B targin 

463.1 

464.1 

469.3 

457.7 

462.9 

464.3 

469.5 

Disp FL 

5485.2 

5297.5 

5417.1 

5234.7 

5353.4 

5484.1 

5530.2 
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Table 35. Naval architectural analysis indices 


CONSTANT CONSTANT 


HECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

NEH 

MR 

ELEC 

12 

NEW 

HR 

ELEC 

20 

FUEL 

ELEC 

23 

BASE 

12 

FUEL 

ELEC 

23 

BASE 

20 

L/D 

11.18 

11.18 

11.18 

11.18 

11.18 

11.18 

11.18 

L/B 

7.73 

7.72 

7.72 

7.72 

7.72 

7.73 

7.73 

B/T 

3.345 

3.425 

3.373 

3.451 

3.401 

3.345 

3.327 

GH1/L8P 

2.235 

2.290 

2.257 

2.305 

2.276 

2.235 

2.220 

NlOO/Dfl 

0.307 

0.318 

0.312 

0.317 

0.311 

0.307 

0.306 

N200/DT1 

0.063 

0.066 

0.071 

0.062 

0.067 

0.064 

0.070 

N300/DT1 

0.043 

0.045 

0.044 

0.045 

0.044 

0.043 

0.043 

W400/DT1 

0.055 

0.057 

0.056 

0.058 

0.056 

0.055 

0.055 

N500/Df1 

0.112 

0.116 

0.113 

0.117 

0.115 

0.112 

0.111 

N600/DT1 

0.078 

0.081 

0.079 

0.081 

0.080 

0.078 

0.077 

moim 

0.017 

0.018 

0.018 

0.018 

0.018 

0.017 

0.017 

Hfuel/Ofl 

0.182 

0.155 

0.164 

0.156 

0.165 

0.182 

0.180 

Wpayload/Bfl 

0.104 

0.108 

0.105 

0.109 

0.107 

0.104 

0.103 

Hld/Dfl 

0.240 

0.212 

0.220 

0.213 

0.222 

0.238 

0.236 

V«b/Vtot 

0.148 

0.148 

0.148 

0.148 

0.148 

0.148 

0.148 

HPi/Dfl 

9.398 

9.731 

9.516 

9.848 

9.629 

9.400 

9.322 

HPitVaax/DTl 

273.76 

282.39 

274.35 

286.47 

278.29 

270.91 

267.53 

kHi/Dfl 

1.094 

1.133 

1.108 

1.146 

1.121 

1.094 

1.085 

Nt /DTI 

0.050 

0.051 

0.050 

0.052 

0.051 

0.050 

0.049 

Vtot/Dfl 

134.914 

139.695 

136.610 

141.370 

138.236 

134.941 

133.816 

moo/vtot 

5.100 

5.103 

5.122 

5.026 

5.045 

5.104 

5.122 

N200/HP1 

14.930 

15.287 

16.821 

14.170 

15.700 

15.291 

16.821 

V»b/HPi 

2.127 

2.127 

2.127 

2.127 

2.127 

2.127 

2.127 

W300/kNi 

0.039 

0.039 

0.039 

0.039 

0.039 

0.039 

0.039 

Vib/(HPi+kWi) 

1.906 

1.906 

1.906 

1.905 

1.905 

1.906 

1.906 

H500/Vtot 

1.863 

1.857 

1.859 

1.857 

1.859 

1.863 

1.863 

H600/Vtot 

1.291 

1.291 

1.291 

1.291 

1.291 

1.291 

1.291 

Dfl/Vtot 

16.603 

16.035 

16.397 

15.845 

16.204 

16.600 

16.739 



































































For this technology characterization, everything 
devolves to total values. What is the total effect on the 
ship, once the individual pieces are put together? Table 36 
gives the bottom line. The electric propulsion plants are 
larger and heavier than their mechanical drive cousin; 
however, the extra weight and volume are more than compen¬ 
sated by the savings in fuel weight and volume. If a ship 
is designed from the beginning to be an "optimized" electric 
drive ship, over 6300 cubic feet of volume and 250 LT may be 
saved. The savings might be used for other systems, to 
reduce the overall size and cost of the ship (maybe allowing 
a larger buy, since 30 ships times 250 LT is a 7500 LT 
ship), or to extend the naval architectural limits of the 
ship design. 

If a ship is backfitted with this technology, it is un¬ 
likely that tank volume can be recovered. However, the dual 
benefits of increased time-on-station and better fuel 
economy are realized. In this case, the choice between 
geared or direct-drive systems can be made by selecting the 
system with the most leverage, i.e., if the ship is volume- 
limited, use the lower volume direct-drive system (since the 
shafts are already in place). 

To put the volume and weight savings in perspective, 
note that 6300 cubic feet and 250 LT translates to twenty 
Tomahawk missile cells. The ship would be volume limited, 
with about 200 LT of weight savings still unused. This is a 
significant addition to the firepower of any ship, and the 
unused weight allows for ship growth. 
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Table 36. Total differences 



HECH 

23 

BASELINE 

ELEC 

23 

BASELINE 

12 

ELEC 

23 

BASELINE 

20 

Total voluies: 
Fuel volume 

44365 

36592 

39647 

Fuel diff 

0 

-7773 

-4718 

Prpln voluae 

3240.75 

5146.33 

4520.05 

Prpln diff 

0.00 

1897.58 

1272.10 

Total 

47613.75 

41738.33 

44167.05 

Total diff 

0.00 

-5875.42 

-3445.90 


Total weights: 

Dfl difference 

0.00 

-187.50 

-67.90 

prpln diff 

0.00 

16.86 

45.83 

fuel diff 

0.00 

-183.80 

-111.60 

other 

0.00 

-20.56 

-2.13 




CONSTANT 

CONSTANT 



FUEL 

FUEL 

NEW 

NEH 

ELEC 

ELEC 

NR 

NR 

23 

23 

ELEC 

ELEC 

BASE 

BASE 

12 

20 

12 

20 


36367 

39396 

44364 

44364 

-7998 

-4969 

-1 

-l 

4936.38 

4310.90 

5125.76 

4500.28 

1687.63 

1062.15 

1877.01 

1251.53 

41303.38 

43706.90 

49489.76 

48864.28 

-6310.37 

-3906.85 

1876.01 

1250.53 


-250.30 

-131.60 

-0.90 

45.20 

-7.67 

21.30 

16.86 

45.83 

-189.10 

-117.50 

0.00 

0.00 

-53.53 

-35.40 

-17.76 

-0.63 
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8.3. Closure 


8.3.1. Conclusions 

This thesis has demonstrated the usefulness of electric 
drive transmissions in reducing ship weight and volume. 
Electric drive transmissions are better than mechanical 
drive transmissions on a ship basis. They provide, besides 
the weight and volume advantages, substantial arrangement 
flexibility and the opportunity to use new technologies in 
the ship design arena. The technical risk associated with 
these different technologies is minimal, as there is much 
industrial experience with electric machines, advanced 
switchgear, and the like. If the weight and volume reduc¬ 
tions are reinvested in the ship design through more optimum 
arrangements and subdivision, a substantially more efficient 
ship may be realized. Such a ship could successfully com¬ 
pete with the best of current ships. 

Small, light, high-power motors can be designed to a 
fair degree of detail with a computer optimization scheme if 
a meaningful objective function can be devised. For a ship 
system, the objective function should contain measures for 
volume, weight, efficiency, and relative cost (if a par¬ 
ticular material is significantly more expensive than other 
used). A steepest-descent scheme can be combined with a 
Monte Carlo scheme to quickly converge on the objective 
function. 

The use of electric drive, and its consequent electri¬ 
cal cross-connect, can reduce the endurance fuel load by as 
much as 17.5%. When used in combination with an improved 
machinery arrangement and subdivision, that percentage can 
rise to 18%. If the fuel load stays constant, the endurance 
range may increase as much as 25%. 

On both an equipment weight basis and a ship weight 
basis, systems composed of a direct-drive propulsion gener¬ 
ator (with the same shaft rpm as the prime mover) and a 
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geared propulsion motor are better than those systems using 
no gears. Regarding volume, a non-geared system has lower 
equipment volume but higher ship volume due to the lower en¬ 
durance fuel efficiency. 

Geared propulsion motors have better off-design-point 
efficiencies than those in direct-drive systems, primarily 
due to their higher rpm. A reduction in output power (in a 
motor) of 75% means only a few percent reduction in 
efficiency for a geared motor, while the same power reduc¬ 
tion means a 20% or more efficiency reduction for a direct- 
drive motor. 

Permanent magnet machines do not appear attractive for 
ship propulsion systems. They are both heavier and larger 
than candidate systems using synchronous, and, to a smaller 
extent, induction machines. Their low air gap flux density 
is the main detractor. Current permanent magnet materials 
cannot develop the energy product to compete with other al¬ 
ternatives, even though the NdFeB magnets are now in the 
marketplace. Induction machines may be useful as propulsion 
motors, but in this thesis they did not appear so. There¬ 
fore, ship propulsion generators and motors should only be 
synchronous machines. 

8.3.2. Recommendations 

The same modeling approach used in this thesis should 
be taken with variable reluctance machines (VRM). Although 
no VRMs have been built at ship propulsion power levels, it 
is not inconceivable that they could serve in such a 
capacity. 

The induction machine model used here needs refinement, 
especially in the area of limiting maximum rotor current 
density. All of the machines need an analysis of their 
transient and dynamic characteristics. 

The recent advent of liquid hydrogen temperature super¬ 
conducting materials may signal an era where conventional 
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machines are overshadowed by the smaller, higher-flux 
machines possible with superconducting technology. However, 
if these new materials fail to provide the required current 
density, design of satisfactory machines may not be pos¬ 
sible . 

Integrated electric ship service propulsion plants may 
be beneficial additions to electric drive technology. Their 
influence on the systems suggested here may be an area of 
interest for future ship designs. 
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Appendix A. Definitions of machine variables and constants. 


bd 

BETA 

br 

BR 

BR1 

BREM 

BRNGS 

BSAT 

CP 

CRHO 

cw 

D 

dcore 

DCU 

DMAG 

doa 

dr 

ds 

eaf 

ef f cy 


ew 

f req 
g 


magnet operating point flux density, Tesla. Used 
in permanent magnet machines. 

hysteresis loss factor of M19 magnetic steel. The 
figure 2.5 was used. 

air gap flux density, Tesla. Used only in per¬ 
manent magnet machine model. 
air gap flux density, 1.05 rms Tesla, 
residual induction, 7.89 kilogauss, of M19 mag¬ 
netic steel,at BSAT and T1. 
remanence flux of NdFeB, 1.21 Tesla, 
weight percentage of rotor shaft bearings, 1.03, 
or three percent of rotor weight, 
max flux density anywhere, 1.5 rms Tesla, 
stator coil pitch. The figure 0.8 was used, 
copper electrical resistivity, 1.724E-8 ohm- 
meters . 

copper weight, kg. 

density, 7.65e3 kg/m3, of M19 magnetic steel. 

back iron depth, meters. 

copper density, 8968.0 kg/meter 3 . 

density of NdFeB, 7.4e3 kg/meter 3 . 

overall machine diameter, meters. 

depth of rotor slots, meters. 

depth of stator slots, meters. 

p.u. internal voltage, used in syn only. 

efficiency, defined as 

output power 

effcy = - 

output power + ph + pe + i2r + i2rr 

effective weight of machine, a combination of 
weight, volume, and efficiency. Used as the ob¬ 
jective criteria for the optimization scheme, 
machine synchronous frequency, Hz. 
air gap dimension, meters. 
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gams 

stator geometric factor, non-dimensional. Used to 

find convergence on active length. 

gamr 

rotor geometric factor, non-dimensional. Used to 

find convergence on active length. 

GMIN 

minimum machine air gap, 0.002 meters. 

HC1 

coercive force, 0.48 oersteds, of M19 magnetic 

steel, at BSAT and Tl. 

il 

induction machine primary current, amps. 

i2 

induction machine load current, amps. 

i2r 

copper loss, watts, on stator. 

i2rr 

copper loss, watts, on rotor. 

jr 

full load rotor current density, amps/m 2 . 

jrnl 

no load rotor current density, used in syn only. 

js 

full load stator current density, amps/m 2 . 

JSMAX 

maximum stator current density, 12e6 amps/meter 2 

ke 

efficiency weighting factor. 

KM 

magnet material cost weighting factor, defined as 

ks[] 

$ per pound of magnet material 

KM - 

$ per pound of magnetic iron 

harmonic winding factors. 

kv 

volume weighting factor. 

1 

machine active length. 

le 

rotor winding space factor, used in induction 

machine. 

ler 

combined length of rotor end windings. 

lm 

radial dimension of permanent magnet. 

loa 

overall machine length, meters. 

lr 

ratio of rotor slot width to slot pitch. In per¬ 
manent magnet machines, defined as ratio of magnet 

width to rotor "slot'' pitch. 

lrat 

length ratio, used in permanent magnet program to 

calculate the effect of magnet overhang. 

Is 

ratio of stator slot width to slot pitch. 
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MAX_TIP_ 

minpwr 

MU 

MUR 

nr 

NU 

overhang 

P 

PC 

pe 

PF 

PI 
ph 
PS I 

r 

rlpl 

r2pl 

rev 

relpl 

RHO 

riv 

RMAX 

rpm 

RSF 

sev 

siv 

slip 


SPEED maximum allowable rotor tangential velocity, 
200 meters per second. 

minimum required power of the machine, horsepower. 
MU of air, 4rre-7 Henries/meter. 

NdFeB relative reversible permeability. The 
figure 1.05 was used. 

number of rotor bars in induction machine, 
anomalous loss factor of M19 magnetic steel. The 
figure of 2 was used, 
percentage of magnet overhang, 
number of machine pole-pairs. 

permanence coefficient of NdFeB. The figure 1.1 
was used. 

eddy current loss, watts. 

power factor of all machines. The figure 0.8 was 
used. 

3.141592654. 
hysteresis loss, watts. 

Trated/Tpullout for induction motor. The figure 

0.55 was used. 

rotor radius, meters. 

per length primary resistance, ohms/meter. 

per length secondary resistance, ohms/meter. 

rotor copper volume, meters 3 . 

per length Thevenin equivalent resistance. 

electrical resistivity, 52 micro-ohm-cm, of M19 

magnetic steel. 

rotor iron volume, meters 3 . 

maximum rotor radius, 2.0 meters. 

machine shaft revolutions per minute, referred to 
as rotor speed. 

rotor slot space factor. The figure 0.35 was 
used. 

stator copper volume, meters 3 . 

stator iron volume, meters 3 . 

guessed slip of the induction machine. 
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slipl 

larger derived slip. 

slip2 

smaller derived slip. 

smax 

maximum machine slip. 

SSF 

stator slot space factor. The figure 0.35 was 

used. 

T1 

thickness, 0.014 inch, of M19 magnetic steel. 

tmaxpl 

per length maximum torque. 

vlapl 

per length Thevenin equivalent voltage. 

va 

iP+jQ!, VA rating, volt-amps. 

vagpl 

vol 

per length air gap voltage. 

machine envelope volume, meters 3 , with a margin. 

VOLALL 

volume allowance for frame and foundation. A ten 

percent allowance was used. 

volmag 

volume of permanent magnet material, meters 3 . 

vtpl 

per length terminal voltage. 

w 

electrical frequency, radians per second. 

wr 

width of rotor slots. 

wm 

mechanical angular velocity of rotor, radians per 

second. 

wt 

weight of copper and iron in a machine, plus a 

margin. 

WTALL 

weight margin for frame and foundation. A ten 

percent margin was used. 

wt_iron 

iron weight, kg. 

wtmag 

weight of permanent magnet material, kg. 

xlpl 

per length primary impedance. 

x2pl 

per length secondary impedance. 

xbeltpl 

per length belt impedance. 

xelpl 

per length Thevevin equivalent impedance. 

xmpl 

per length air gap magnetizing impedance. 

xrdpl 

per length rotor differential leakage impedance. 

xrspl 

per length rotor slot leakage impedance. 

xs 

p.u. synchronous impedance. 

xsepl 

per length stator end turn impedance. 

xslot 

slot impedance. 

xsdpl 

per length stator differential leakage impedance. 
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xsspl 

xzzpl 


per length stator slot leakage impedance, 
per length zig-zag impedance. 
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Table 37. Listing of various functions used in computer 
programs 


Ittt a pseudo randoa nuaber generator tttl 

Idefine 

MULTIPLIER 

25173 

tdefine 

MODULUS 

3276B 

Idefine 

INCREMENT 

13849 

Idefine 

H0DFLT 

32768.0 

double randoaO 



{ 


extern long int seed; 

seed=(MULTIPLIERtseed+INCREMENT) Z MODULUS; 
return(seed/MODFLT); 

} 


double swf(n) It stator winding factor tl 

int n; It haraonic order tl 

{ 

double kp, kb; 

extern double cost), sin(); 

kp=cos(0.3142ln); It pitch factor, assuaes 0.8 coil pitch tl 

kb=(sin(0.5236Jn))/(0.5236tn); It breadth factor, frot 

Kirtley’s "Basic Formulas ... 8 and assuaes an 
electrical winding angle of 60* tl 

return(kptkb); 

} 


double abs(q) 

double q; 

{ 

if (q < 0.0) 

q= (-1.0)tq; 
return(q); 

} 


double sinh(u) 

double u; 

{ 

double exp(), ans; 

ans = 0.51(exp(u) - exp(-u)); 
return(ans); 

} 


double cosh(u) 

double u; 
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{ 

double exp(), ans; 

ans = 0.51(exp(u) + exp(-u)); 
return(ans); 

} 

double rwf(n) It rotor winding factor, saae conents as swfi) tl 

int n; 

{ 

double kp, kb; 
extern double sin(); 
kp=l; 

kb=(sin(0.5236tn))/(0.5236ln); 

return(kptkb); 


Idefine toi l.Oe-6 
Idefine pi 3.141592654 

double besi (p,x) It bessei I(order, arg) tl 

double x; 
int p; 

{ 

double bi; 
int bia(); 

double exp(),sqrt(),bi5(); 
if (p<0) abort ( B besi: negative index'); 
if (x<0) abort (' besi: negative argument*); 
if (x>60) bi=exp(x)/sqrt<2lpitx); 
if ((x<12> i ibia(p,x,ttbi)) bi=bis(p,x); 
return (bi); 


double bis(p,x) 
double x; 
int p; 

{ 

double tabs!); 
double xx; 
int i,fk,k; 
double bi=0; 
double t=l; 
xx=x/2.; 

for (i=l;i<p+l;i++) t=tlxx/i; 

if (t<l.0e-36> bi=0; 

else 

{ 

bi=t; 

xx=xx$xx; 

for (k=lj(k<100mt(ffab5(t)-fab5lbiltol))>0);k++) 
{ 

fk=kt(p+k>; 


lo 5 




t=tlxx/fk; 

bi=bi+t; 

} 

} 

return (bi); 

> 

int bia(p,x,pbi) 
double x,tpbi; 
int p; 

{ 

int ist,fn,k,fk; 

double xx,bi,t; 

double fabs(),sqrt(),exp(); 

fn=4ip*p; 

t=1.0; 

bi=l.0; 

xx=0.125/x; 

{or (k=l; (k<30)U((fabs(t)-tabslbittol)) >0) ;k++) 

{ 

fk=(2lk-l)l(2lk-l); 

t=tlxxt(fk-fn)/k; 

bi=bi+t; 

} 

if lk==31) ist=l; 
else 

{ 

ist=0; 

bi=bilexp(x)/sqrt(2.0lpilx); 

> 

Jpbi=bi; 
return (ist); 

> 

double besk(p,x) It aodified bessel function Kptx) 1/ 

double x; 
int p; 

{ 

double bk; 

double exp(),sqrt(),kO(),kl(); 

if (p<01 abort (’ besk: negative index *); 

if (x(0) abort (* besk: negative argument *); 

if (x>60) bk=exp(-x)/sqrt(2.0lx/pi); 

else if (p==0) bk=k0(x); 

else if (p==l) bk=ki(x); 

else f 

double g0,gl,gj; 
int j; 
g0=k0(x); 
gl=kl(x); 

for (j=2;j<p+I;j++) 

{ 

gj=2l(j-l)lgl/x+g0; 

g0=gl; 
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9t=gj; 

) 

bk=gj; 

\ 

i 

return (bk); 

} 

double kO(x) 
double x; 

{ 

double bk; 

double log(),sqrt(l,exp(); 
if (x<l) 

{ 

double a,b,:,c,g0,x2j,f,hj,rj; 

int j; 

b=0.5lx; 

a=.57721566+1og(b); 

c=blb; 

gO=-a; 

x2j=l; 

f=l; 

hj=0; 

for (j=l;j<7;j++) 

{ 

rj=1.0/j; 
x2j=x2jtc; 
f=flrjlrj; 
hj=hj+rj; 

g0-g0+x2j If t (hj-a); 

) 

bk=gO; 

) 

else 

{ 

double t[12],a,b,c,pa; 

int 1; 

a=exp(-x); 

b=1.0/x; 

c=sqrt(b); 

tC03=b; 

for (1=1;1<12;1++) tUMU-Ulbt 

pa = 1.2533141-. 1566642IU01; 

pa += .0B811128tt[11-.09139095* t£2J; 

pa += .1344596ttl33-.2299850lt[4); 

pa += .379241tt[53-.5247277$tt63; 

pa += .557536BttI73-.4262633ltt03; 

pa += .218451B*t[9]-.066B0977*tI103+.009l893B3tt£113 

bk-atctpa; 

} 

return (bk); 
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double kllxl 
double x; 

{ 

double bk; 

double log(),exp(),sqrt(); 
if (x<l) 

( 

double a,b,c,ql,x2j,f,hj,rj; 

int j; 

b=x/2; 

57721566+log(b); 
c=blb; 
x2j=b; 
fM; 
hj=l; 

qI= 1.0/x+x 2j<(.5+a-hj); 
for (j=2;j<9;j++) 

{ 

x2j=x2jlc; 
rj=1.0/j; 
f=fIrjlrj; 
hj=hj+rj; 

gl=gl+x2jlfl(.5+(a-hj>lj); 

} 

bk=gl; 

} 

else 

{ 

double a,b,c,ttl2],pa; 

int 1; 

a=exp(-x); 

b=1.0/x; 

c=5qrt(b); 

tfO]=b; 

for (1=1;1<12;1++) til]=tll-l]*bj 

pa = 1.2533l4l+.4A99?7tt[0]-.1468583lttl]; 
pa += .1280427ltt2]-.1736432lt[3]+.2847618lt[41j 
pa ♦* -.4594342ltt5]+.62B338llt[6]-.4632295Ut7]; 
pa += .5050239lt[8]-.258l304ltt91; 
pa += .0788000ltttl0] 0I082418ltt113; 

bk=atctpa; 

} 

return (bk); 

1 

double besip (p,arg) /I deriv of bes 1(order, arg) 1/ 

int p; 
double arg; 

{ 

double x,y,z; 
x = besi (p-1,arg); 
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y = p I besi (p,arg)/arg; 
z = x - y; 
return (z); 

} 

double beskp(p,arg) It deriv of bes K(order, arg) 1/ 

int p| 
double arg; 

{ 

double x,y,z; 

x = -besk (p-1,arg); 

y = - p t besk (p,arg)/arg; 

z = x + y; 

return (z); 

> 
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Appendix B. Synchronous Machines and General Relations 
An equivalent circuit for a synchronous machine is 


Ra Xs 



Figure 35. Synchronous machine equivalent circuit 

where Ra is the stator resistance and Xs is the synchronous 
reactance. The internal voltage of the machine, Eaf, is 
developed between the stator and rotor across the air gap. 
The current direction is shown as if the machine is a motor. 

Vta = Rala + jXsIa + Eaf 

Eaf represents a mutual coupling between the stator and 
rotor, and 


w M If 4 Po ki kf 1 r Ns Nf 

Eaf =- where M =- 

42 tt g p* 


If is field current, Nf is the number of series field turns 
and kf is the fundamental rotor winding factor. Since it is 
never desired to specify the number of turns on either the 
stator or the rotor, a scheme has been devised so that 
derivations are conducted in volts-per-turn, ampere-turns, 
and ohms-per-turn-squared, which results in power in watts. 
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Figure 36. Phase belt conductor area 


For the stator, the transverse area can be divided into 
phase belts. The area labelled *A’ is the area occupied by 
phase 'a’, in one direction. Then, Ns la = Aa Ja, where Ja 
is the current density in phase 'a' conductors. The area 
has some effective conductor area, subject to the need for 
conductor cooling passages and insulation area. Therefore, 
for stator and rotor currents, analysis yields 


la = 


Aa Ja 
Ns 


and 


If = 


Af Jf 
Nf 


It is also desired never to specify the number of slots on 
either the stator and rotor. Accordingly, slot space fac¬ 
tors are defined as 

U slot J 

r width i 


k 


slot pitch 


Figure 37. Slot space factors 

stator slot width 

ls = - 

stator slot pitch 
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rotor slot width 
rotor slot pitch 































































For a typical turbogenerator conductor bar, the copper area 
is about thirty-five percent of the conductor envelope area. 
Variables titled SSF (stator slot space factor) and RSF 
(rotor slot space factor) embody this thirty-five percent. 
The conductor area of a single stator phase in one direction 

(r+g) 2tt Is ds SSF 

Aa =- and, for the rotor, 

6 

r 2tt lr dr RSF 

Af =- There is only one phase on the 

2 rotor. 

The armature resistance is 

rho It Ns 

Ra =- 

Aa/Ns 

where rho is the electrical resistivity of copper and It is 
the turn length. If a stator turn can be modeled as 



Figure 38. Stator turn 


and the circumfrence is 2rr(r+g), then 


It 


21 + 2 { 


2tt _CP_-£ ( 4/3 > 

P 


> 


where -1(4/3) is (l/sin(60 0 )) and CP is the coil pitch. Then 

Ra __ - _2_£ho_C_l_+__C 2 w i( r +«>CP/p)*£C4/’3) >_ 

Ns 2 it (r+g) Is ds SSF 
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The stator and rotor copper losses are 


Rala 2 - rho It Aa Ja 2 Rflf 2 = rho ltf Af Jf 2 

The rotor is full pitched. 

The vector diagram for the synchronous machine equiv¬ 
alent circuit is (assuming Ra ==> 0) 



Figure 39. Synchronous machine vector diagram 

By use of the law of cosines, 

Eaf 2 = Vt 2 + (Xsla) 2 - 2VtXsIa cos(r + 90°) 

where r is the power factor angle. If this equation is put 
into per unit form, with eaf = Eaf/Vt and xs = Xsla/Vt, then 

eaf 2 = 1 + xs 2 + 2 xs sin(r) 

xs can be calculated, allowing the calculation of eaf. 

There is simple linear relationship between eaf, the no-load 
rotor current density (jrnl), and the full-load rotor cur¬ 
rent density (jr). It is jr = (eaf)(jrnl), because Eaf is 
directly proportional to If and If is directly proportional 
to Jf . Ampere’s Law states 



If one chooses an integration path around half of the rotor, 

4 Nf If kf 

Po H 2g - Po-- 

TT P 
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If the no-load condition is chosen, and it is recognized 
that Br is essentially constant with increases in Vt 
(because when Vt increases, so must Eaf), then 


Br 


2 Mo r lr dr RSF kf jrnl 
g P 


and jrnl 


_Br_g_p_ 

2 Mo r lr dr RSF kf 


The synchronous reactance is a measure of flux linkage 
within the machine. It is Xs = w Ls = w (Lai + 1.5 Laa), 
where Laa is the armature single phase inductance and Lai is 
the slot leakage inductance. The factor of 1.5 is derived 
from the 120° separation between the three phases of these 
machines, as 


Xsla = w (la Lai +IaLaa + IbLab + IcLac) 


Due to the symmetry of the machine, 


Lab = Lba = Lac = Lea = Lbc = Lcb = -0.5 Laa 


Then, Xsla = w la (Lai + 1.5 Laa). 

Self-inductance is 


4 Mo ks 2 Ns 2 1 r 

Laa =- 

tt g P 2 


from Ampere’s Law and L = (flux/current). If a single con¬ 
ductor per slot is postulated and the effects of slot teeth 
are ignored, then the leakage inductance is 


slots 

Lai =- (Pself + Pmutual) Ns 2 

(pole)(phase) 
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where P is reluctance. For this conductor configuration, 


Mo 1 

Pself = ds ws 

3 

Mo 1 

and Pmutual = ds ws, 

2 

with ws equal to slot width. 

Since ws is not known, we use 

the stator slot space factor, 

Is, multiplied by the number 

of slots to yield 


(slots) ws = Is tt (r+g) 

Therefore, 


5 Ns 2 Mo 1 ds Is tt (r+g) 
Lai =--- 


36 p 


The real power developed 

by the machine is 

Pwr = 3 Vt la pf, where pf is 

the power factor. Through the 

use of Lens’ Law and Ampere’s 

Law, terminal voltage may be 

expressed as 


2 r 1 w Br ks 

Vt =- 

P 



Using our previous relation for la, the expression for Pwr 


is 


2tt r 1 w Br ks Ja 

Pwr =- 

P 

SSF (r+g) Is ds pf 

Finally, winding factors 

need to be derived. The wind- 

ing breadth factor, kb, is 


sin(m n P/2) 

kbn = 

m sin(n r/2) 



where P is the electrical angle between adjacent slots, n is 
the harmonic order, and m is the number of slots per pole 
per phase. The winding pitch factor, kp, is 
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kpn = sin (n a/2) where a is the electrical angle 


between sides of the coils (pitch angle). For a three phase 
winding, a = 2 tt p. If a 0.8 coil pitch is assumed (which 
will rid the machine of certain harmonics during balanced 
operation) then a = 180° (1 - 0.8) = 36° and 
kpn = sin(0.3142 n). Assumptions are needed to calculate 
kbn without specifying the number of stator slots or turns. 
If the winding angle is specified as did Kirtley [32], then 


kb n 


sin(n Ow/2) 
(n ©w/2) 


A reasonable electrical winding angle is 60°, since most of 
the stator periphery will contain turns. The breadth factor 
devolves to kbn = sin(0.5236 n)/(0.5236 n), for which the 
fundamental harmonic factor equals 0.955. The winding fac¬ 
tor is the product of the breadth and pitch factors. For 
the rotor, a pitch factor of one is assumed. 
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Table 38. Listing of synchronous machine design program 


•include "stdio.h" 
•include "def.h" 


ft prograa naae: syn.c for synchronous, round rotor aachmes tf 

long int seed; ft start point for randoa nuaber generator tl 

double bC263£261, hCHJC2&] 9 ks[81, kr[81; 

ft b is “best" array, h is "hold" array, ks/kr are winding factors tf 

•ain() 

{ 

double design^point0, rnd_walkl), swf(), rwf(), ke, kv, ainpwr, 
stepsize, randoaO, abs(), freq, rpa; 
int p, iteration, i, j, best, printout0, loops; 

FILE tfopenO, tfp; 

printf("\nReading input data froa SYN.DAT . . ."); 


fp=fopen("syn.dat","r"); ft input seed for randoa nuabers tl 
fscanf(fp,"Xd",iseed); 
fscanf(fp, "Xd", ip); 
fscanf (fp, "Xlf", iainpwr); 

■inpwr1=746.0; 
fscanf(fp, "Ilf", ike); 
fscanf (fp, "Xlf", ikv); 
fscanf(fp, # Xlf i , irp»); 
fdose(fp); 


ft input nuaber of pole pairs tf 
ft input aachine power, derived fa ASSET tf 
ft convert to watts tf 
ft CERs for Effective Weight tf 

ft aachine aax shaft rpa tf 


printf( B \nHow aany loops do you want? "); 
scanf("Id", iloaps); 

printf( a \n\nDoing prograa calculations . . .\n"); 


for (i=t; i < 8; i+=2) ft haraonic winding factors tf 

( 

ks[ i]=5wf(i); 
krC i]=rwf(i); 

} 


freq=rpa$p/60.0; ft aax electrical frequency tl 

ft MAIN BODY OF THE PROGRAM tf 


for (i=1; i <= loops; ++i) 

( 

stepsize=0.1; 

iteration^; 

design^point(ainpwr, p, ke, kv, freq); 

ft put stuff in the hold array tf 

while (iteration <= 10) 

{ 
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rnd_xalklainpxr, p, stepsize, ke, kv, freq); 

/I stagger around 1/ 

best=0; It index to best EH of the lot 1/ 

for (j=l; j<=10; ++j) 

if (h[j][181 < hEbestj[18]) 

best=j; It find the best «achine 1/ 

if labsllh[0][181 - h[best3[18J)/h[0][18]) < 0.005) 

It saall iaproveaent in EH I/ 

{ 

stepsize/=2.0; 

♦♦iteration; 

) 

else /» transfers best to 0 position 1/ 

{ 

for (j=l; j <= 25; ++j) 

h[0)[j) = hibestlljl; 

} 

} 

for (j=l; j <= 25; ++j) 

bti][j]=h[best][j]; It keep the best aachine t/ 


best=l; 

for (i=1; i <= loops; ++i) 

if (b[i][181 < b[best 1[183) 

best=i; It find and keep the best of the best 1/ 

ainp»r/=746.0; It turn back into hp 1/ 

print_out(best, p, ainpxr, ke, kv, rpa); It output to disk file I/ 


fp=fopen("syn.dat*,V); It output seed 1/ 

fprintf(fp,"Zd*, seed); 

fprintf(fp,“Sn2d", p); 

fprintf(fp,’\nJlf", ainpxr); 

fprintf (fp,’\nXlF, ke); 

fprintf(fp, B \n21f', kv); 

fprintf(fp,"\nXlf“, rpa); 

fdoselfp); 

} 

It END OF MAIN PROGRAM; ALL THAT FOLLOW ARE FUNCTIONS 1/ 


double desiqn_point(ainp»r, p, ke, kv, freq) 

It deteraines a randoa design point >/ 

double ainpur, ke, kv, freq; 
int p; 

{ 

double r, jrnl, jr, js, Is, lr, dcore, ds, dr, g, », 1, xs, eaf, i2rr, 
va, ph, pe, i2r, vol, wt, effcy, ex, xsl, xs5, xs7, xsal, 
cx, siv, scv, riv, rev, loa, doa, xt_iron, find_sivl); 
extern double sqrtO, randoaO; 
int c=0, d=0; 

xhile Id != 1) { 
xhile Ic != 1) 
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{ 

r-randoiOIRPIfiJ; 

if ((2IPIIrlfreq/p) < HAX_T1P_SPEED) /I check rotor tip speed I/ 
break; 

> 


w=2tPl*freq; 


lr=randoa()!0.5 + 0.25; 

It rotor slot factor 1/ 

dr=randoa()tr/5.0; 

It slot no deeper than 20 l of rotor radius tl 

dcore=(BR*r)/(BSATtp); 

It back* iron depth 1/ 

ds'randoa()t0.9tdcore; 

It slot depth < 901 of body depth 1/ 

while (c != 1) 

It gap dimension 1/ 


{ 


g=randoa()l(0.1»r - GH1N) + SHIN; 
if (g > 0) 

break; 

> 

ls=randoa()!0.5 + 0.25; IX stator slot factor XI 

js=randoa()JJSHAX; IX full load stator current density XI 

jrnl = (BRtgtp)/(8tHUIRSFtrldrtlrtkrU3); IX noload rotor current density XI 

xsl=ks[H*ksU]; 

xs5=(ks[51lks[5]/25); 

ss7=(ksC73tkst73/49); 

xsal=(5tlslPllds/l8); 

xs=(HU»jstSSFt(r+g)*(r+g)ldslPllls I (12l(xsl + xs5 + xs7)/(Pl*glp) + xsal)> 

/(12*rIBRIksl13); 

IX p.u. synch iapedance XI 

if (xs > 2.0) 

continue; IX don't want xs too big XI 

else 

++d; 

eaf=sqrt(l + xstxs + 2txst0.fc); IX 0.6 is sin(r), pur factor angle, 
eaf is p.u. internal voltage at full load XI 

jr=eafJjrnl; IX jr full load, linear with eaf XI 

l=(ainpwrtp)/(2»BRtw»ksU]trlSSFtdstP[t(r+g)ll5tjsSPF); IX active length XI 

va=2tPl*rtltwIBRlksCl3<js<SSF*(r+g>*dslls/p; IX va rating XI 

siv = find_siv(l, r, g, ds, dcore, Is); IX stator iron voluae 1/ 

riv = UPltrMr - 2tdrllr); IX rotor iron voluae XI 

scv = 2*PIt(r+g)»ds»ls»(l + 2.3094tPll(r+g)tCP/p); 

IX stator copper voluae XI 
rev = 2IPltrtdrtlrl(l + 2.3094tPltr/p); IX rotor copper voluae XI 

cm = (rev + scvllDCU; IX total copper weight XI 

loa = 1 + 4t(r+g); IX length-over-all XI 

doa = 21(r+g+ds+dcore); IX over-all-diaaeter XI 

vol = VQLALL4(1oatP1fdoatdoa/4); IX aachine envelope voluae XI 

wt = HTALL* (cm + DKBRNGStriv + siv)); IX aachine weight in kg XI 
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wtjron = DMriv ♦ siv); 


IX iron weight only 1/ 


ph = 31.86225IBETAIfreqlBRllHCllwt Jron/D; 

It hysteresis loss in watts, uses iron weight of aachine 1/ 

pe = (106236.9INUIBSATIBSATIT1IT 1Ifreqlfreqtwt_iron)/(RHOID)? 

It eddy current loss in watts, uses iron weight of aachine X/ 

i2r = 2.01SSFIdsICRHOlPII (1 ♦ 2.3094IPI$(r+g)ICP/p) Ijsljsl(r+q>lls; 

IX stator copper loss in watts 1/ 
IX revised 1-12-87 XI 

i2rr = 2.OIRSFtdrtrIIrtjrIjrICRHOtP11 (1 + 2.30941PI*r/p); 

IX rotor excitation losses, 1-12-87 XI 
i2r+=i2rr; /I total copper losses XI 


effcy=(ainpwr)/(ainpwr + ph + pe + i 2r); 


ew=wt + kel(l-effcy) + kvtvol; IX Effective weight 1/ 


hEOHIHs; 

hE03 E23=freq; 

hE03 E33 =h; 

ht03E4]=r; 

h E03 E51=9j 

ht03E63=dcore; 

hE03 E73=ds; 

hE03[8]=dr; 

hE03E91=ls; 

hEOlE103=lri 

hE03Ell]=vol; 

h[03 E123=wt; 

hE03 E133=ph; 

hE03E143=pe; 

ht03tl53=i2r 

hE03E163=va; 

hE03 E173=effey; hE03E183=ew; 

hE03 E193=1; 

hE03E203=jr; 

hE0][21]=jrnl; 

hE03E22]=xs; 

hE03E233=eat; 

hE03E243=loa; 

ht03t253=doa 


IX this section just changed all the variables in the “hold 8 array XI 


return; 

» 


double rnd_walk(ainpwr, p, stepsize, ke, kv, freq) 

It walks about design_point 10 tiaes 1/ 
double stepsize, ainpwr, ke, kv, freq; 
int p; 

{ 

double r, jrnl, jr, js, Is, 1r, dcore, ds, dr, g, w, 1, xs, eaf, i2rr, 
va, ph, pe, i2r, vol, wt, effcy, ew, xsl, xs5, xs7, xsal, 
cw, siv, scv, riv, rev, loa, doa, wtjron, find_siv(); 
extern double sqrtO, randoaO; 
int i=l; 


while (i <= 10) 

{ 

It read in the walk around the design point 1/ 

js=ht03[l]t(l ♦ stepsizeMrandoaO - 0.5)); 
if (js > JSrtAX) 

js = JSflAX; It reset to liait 1/ 

w=2IPItfreq; 

r=h[0][43l(l ♦ stepsizeKrandoaO - 0.5)); 
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if <(air/p) > HAKJIP.SPEED) 

continue; It go to next try if violated 1/ 

g=hCOlE53H1 + stepsizeKrandoaO - 0.5)1; 
if (g < 6H1N) 

g=GH]N; It reset to the li«it 1/ 

dcore=(BRIr)/(BSATtp); It aost efficient use of iron tl 

ds=h[0][7]t(l + stepsizeKrandoaO - 0.5)); 

if (ds > dcore) It can’t have too-deep slots tl 

ds=dcore; It reset to the liait 1/ 

dr=h[0][81l(l + stepsizeKrandoaO - 0.5)); 
ls=hC0]C9]l(1 + stepsizeKrandoaO - 0.5)); 
if (Is > 0.75) 

ls=0.75; It reset to the liait 1/ 

if (Is < 0.25) 

1 s=0.25; 

1 r=hC03C10]t(1 + stepsizeKrandoaO - 0.5)); 
if (lr > 0.75) 

1r=0.75; It reset to the liait 1/ 

if (lr < 0.25) 

lr=0.25; 

It coaputation section of the walk 1/ 

jrnl=(BR*gtp)/(BIHUtRSFtrtdrllrlkrt11); It noload rotor current density tl 

xsl=k$[ims[l}; 
xs5=(ksC53*ksC5]/25); 
xs7=(ks[71tkst7]/49); 
xsal = (5llslPIIds/18); 

xsKHUIjstSSFKr+glKr+gltdsIPIlls I (12Kxsl + xs5 + xs7)/(PItglp) + xsal)) 

/(12*rIBRtksC 13); 

It p.u. synch iapedance tl 

if (xs > 2.0) 

continue; It can’t have xs too big tl 

eaf=sqrt(l + xstxs + 2lxst0.6); It O.i is sin( t), par factor angle, 

eaf is p.u. internal voltage at full load tl 

jr=eaftjrnl; It jr full load, linear aith eaf 1/ 

lKainpartp)/(2JBRtaJk5[l]lr»SSFidstPII(r+g)tlsljstPF); It active length 1/ 

va=2*PI*r*llwtBRIksC13*jstSSF*(r+g)tdstls/p; It va rating 1/ 


siv = find_siv(l, r, g, ds, dcore, Is); 
riv = HPIlrKr - 2ldrllr); 
scv = 2IPII(r+g)tdstlst(1 + 2.3094*PIt(r+g)*CP/p); 


rev = 2IPI»rldrllrKl + 2.3094»PItr/p); 

ca = (rev + scvlIDCU; 

loa = 1 + 4Kr+g); 

doa = 2Kr+g+ds+dcore); 

vol - VOLALLKloaiPIldoaldoa/4); 

at = HTALLKca + DKBRNGSIriv + siv)); 


It stator iron voluae tl 

It rotor iron voluae tl 

It stator copper voluae tl 
It rotor copper voluae tl 
It total copper aeight tl 
It length-over-all tl 
It over-all-diaaeter tl 
It aachine envelope voluae tl 

It aachine aeight in kg tl 
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wtjron = DMriv + si v); 


It iron weight only tl 


ph = 31.86225tBETAIfreqlBRliHCllwtjron/D; 

It hysteresis loss in watts, uses iron weight of tachine 1/ 

pe = (106236.9tNUtB5ATI8SATtTltilIfreqlfreqtwt_iron)/(RHOID); 

It eddy current loss in watts, uses iron weight of «achine 1/ 

i2r = 2.0ISSFtdstCRH0tPII (1 + 2.3094IPII(r+g)tCP/p) Ijsljst(r+g)tls; 

ft stator copper loss in watts tl 
It revised 1-12-87 tl 

i2rr = 2.0IRSFIdrlrtlrtjrljrtCRH0tPll (1 + 2.3094IPIIr/p); 

It rotor excitation losses, 1-12-87 tl 
i2r+=i2rr; It total copper losses tl 

effcy=(«inpwr)/(ftinpwr + ph + pe + i2r); 

ew=wt + ket(l-effcy) + kvlvol; It Effective weight tl 


hCi][l]=js; 

hEi H23=f req; 

hEiDE31 =h; 

hEi][4]=r; 

hti3E53=g; 

hti][6]=dcore; 

h[i3[7]=ds; 

hEi][8D=dr; 

hti3 E9]=ls; 

hti3E103=lr; 

h EiIt 11D=vol; 

hii3 E123=wt; 

hEi3E13]=ph; 

hEi31143-pe; 

hti 3E153=i2r 

hCiIE16}=va; 

hCi]C17]-effcy; 

hEi3E18]=eM; 

hti3E193=1; 

hEi3E203=jr; 

hEi3E21 ]=jrnl; 

hEi3[223-xs; 

hEi H233=eaf; 

hEi3E243=laa; 

h£i3£253=doa 


It this section just changed all the variables in the "hold 8 array tl 


++i; It go to the next hti][] tl 

} 


return; 

> 


printout (best, p, tinpwr, ke, kv, rp«) 
int best, p; 

double linpwr, ke, kv, rpi; 

{ 

char outfiletMl; 

FILE Ifpo, tfopenO; 
int i; 

printf("\ntfhat is the nase of the file where you want the output? "); 
scanf("Is", outfile); 

fpo=fopen(outfile, V); 


fprintf(fpo,-%d", p); 
fprintf(fpo,“\nXlf", «inpwr); 
fprintf(fpo,"\nllf", ke); 
fprintf(fpo,"\n21f", kv); 
fprintf(fpo,"\nllf", rp«); 
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for (i=i; i <= 25; ++i) 

fprintf (fpo, "NnXHbtbest3ti ]); 
fprintflfpo, *\n‘); 
fclose(fpo); 

) 


double find siv( 1 , r, 9 , ds, dcore, Is) 
double 1 , r, g, ds, dcore, Is; 

{ 

double one, two, three, four; 

one = lr+g+ds+dcore)tlr+g+ds+dcore) - lr+g)tlr+g); 
t*o = 2tPIt(r+g)tdstls; 

three = (r+g+ds+dcore>*(r+g+ds+dcore) - lr+g+ds)f(r+g+ds) 

four = IMPitone - tNo) + PII4ilr+g)tthree; 
returnlfour); 

} 
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Table 39. Listing of synchronous efficiency program 


tinciude “stdio.h" 
■include H def.h B 


IX progra* naae: seft.c to find efficiency of synchronous aachines XI 
It works with only a single aachine 1/ 

aainO 

{ 

FILE IfopenO, ifp; 

double r, jrnl, jr, js, Is, lr, dcore, ds, dr, g, w, 1, xs, eaf, i2rr, 

va, ph, pe, i2r, vol, wt, effcy, ew, xsl, xs5, xs7, xsal, ks[81, 
siv, riv, wt_iron, find_siv(), paeff, pajs, parpa, dhp, rpa, ainpwr, 
ke, kv, freq; 

extern double swf(), sqrtO; 
int e=0, f, p, i; 
char infiletl43; 

for (i=1; i < 8 ; i+=2) It haraonic winding factors XI 

ksti]=swf(i); 

} 

X5l=kstl3ikstl3; 

xs5=lks[53iks[53/25); 

xs7=iksC73tks[73/49); 

xsal=(5llslPl*ds/lB); 


printf(“\nCalculates efficiency of a single aotor.\n - ); 
while (e != 1) 

{ 

f = 0; 

printf PHhat is the naie of the input file? “); 
scanfrZs", infile); 
fp = f open (inf ile, V); 

fscanflfp, "2d", (p); It input number of pole pairs tl 

fscanflfp, “ZIP, (iinpwr); 

•inpwr 1= 746.0; /1 now in watts tl 

fscanflfp, “ZIP, (ke); 

fscanflfp, 'Ilf, (kv); 

fscanflfp, “ZIP, (rp«); 

fscanflfp, "HP, (js); 

fscanflfp, "Xlf", (freq); 

fscanflfp, "ZIP, (w); 

fscanflfp, 'ZIP, (r); 

fscanflfp, “ZIP, (g); 

fscanflfp, “ZIP, (dcore); 

fscanflfp, “ZIP, (ds); 

fscanflfp, “ZIP, (dr); 

fscanflfp, “ZIP, (Is); 

fscanflfp, “ZIP, (lr); 

fscanflfp, “ZIP, (vol); 

fscanflfp, “ZIP, (wt); 


154 























































fscanfIfp, 'Xlf', iph); 
fscanflfp, "Xlf, &pe); 
fscanf(fp, ■Ilf , fci2r); 
fscanf(fp, "XU", iva); 
fscanf (fp, "Xlf", ieffcy); 
fscanf(fp, 'Xlf', lew); 
fscanf (fp, 'Xlf, il); 
fscanf(fp, 'Xlf', lor); 
fscanf(fp, 'Xlf', fcjrnl); 
fdose(fp); 

siv = findsivd, r, g, ds, dcore, Is); It stator iron voluae 1 / 

riv = HPUrKr - 2Jdrllr); /I rotor iron voluae 1 / 

wt_iron = DKriv + siv); /I iron weight only 1/ 

while If != 1) 

{ 

printf(’\nWhat is the sustained speed aachine horsepower? ’); 
scanffXlf', idhpi; 
dhp 1= 746.0; 

printf(“What is the sustained speed aachine rpa? "); 
scanffXlf, iparpa); 

freq = parpalp/60.0; It aax electrical frequency II 

pajs = jsldhp/ainpwr; It PH stator current 1/ 

xs=(HUIjslSSFIlr+g)t(r+g)IdsIPlIls I (12*(xsl + xs5 + xs7)/(Pllglp) + xsal)) 

/<12*r*BR*ks£13); /I p.u. synch iapedance 1/ 

eaf=sqrt11 + xslxs + 2lxsl0.6); It 0.6 is sin(r), pwr factor angle, 
eaf is p.u. internal voltage at full load tl 
jr=eafIjrnl; It jr full load, linear with eaf tl 

ph = 31.86225IBETAIfreqtBRllHCllwt_iron/D; 

It hysteresis loss in watts, uses iron weight of aachine 1/ 
pe = (106236.9INUIBSfiTIBSfiTITllTltfreqtfreqlwt_iron)/(RH0iD); 

It eddy current loss in watts, uses iron weight of aachine tl 
i 2 r = 2.0*SSF*ds*CRH0*PE* (1 + 2.3094IP1*(r+g)*CP/p) Ijsljsllr+gltls; 

It stator copper loss in watts 1/ 
i2rr = 2.OIRSFIdrlrllrljr*jrlCRHOlPI* (l + 2.3094IPIIr/p>; 

It rotor excitation losses, 1-12-87 tl 
paeff = dhp/ldhp + ph + pe + i 2 r + i 2 rr); 
printf(”\n Sustained speed efficiency is Xlf, paeff); 


printfC\nHhat is the endurance speed aachine horsepower? '); 
scanffXlf, fcdhp); 
dhp t= 746.0; 

printf I'Hhat is the endurance speed aachine rpa? '); 
scanffXlf, iparpa); 

pajs = jsldhp/ainpwr; It PM stator current tl 

freq = parpalp/60.0; It aax electrical frequency tl 

xs=(HUtjstSSFIIr+g)l(r+g)lds!PIIls I (12*(xsl + xs5 + xs7)/(PIIgtp) + xsal)) 

/(12lrlBRU:s[ll); It p.u. synch iapedance tl 

eaf=sqrt(l + xslxs + 2lxsl0.6); It 0.6 is sin(t), pwr factor angle, 
eaf is p.u. internal voltage at full load tl 
jr=eafIjrnl; It jr full load, linear with eaf tl 

ph = 31.86225* BETAIfreqlBRlIHC1twt_iron/D; 
pe = (106236.91NU*BSATIBSAT*THTllfreqtfreq*wt iron)/(RH01D); 
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i2r = 2.0tSSFtdstCRH0tPlt (1 + 2.3094tPIt(r+g>»CP/p) tjstjst(r*q)lls; 
i2rr = 2-OtRSFIdrtrfir IjrljrICRHOtPI$ (1 + 2.3094IPIlr/p>; 
paeff = dhp/(dhp + ph * pe + i 2 r + i 2 rr); 
pnntf("\n Endurance efficiency is XIf*, paeff); 

printf("\nSaue aachine 7 *); 
scanf{"Xd\ if); 
if (f == 0) 

continue; 
else if (f == 2 ) 

{ 

e = 1 ; 
break; 

} 

} It end of f-loop 1/ 

} ft end of e-loop Xi 

} it end of tain proqraa ii 

double find_siv(l, r, q, ds, dcore, Is) 
double 1 , r, q, ds, dcore, Is; 

{ 

double one, two, three, four; 

one = (r+g+ds+dcore)l(r+g+ds+dcore) - (r+g)t(r+g); 
two = 2IPII<r+g)lds$ls; 

three = (r+g+ds+dcore)l(r+g+ds+dcore) - (r+g+ds)l(r+g+ds); 

four = iKPItone - two) + PIMJ(r+g)tthree; 
return(four >; 

} 
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Appendix C. Permanent magnet machine 


The equivalent circuit for a permanent magnet machine 
is almost identical to that of a synchronous machine. The 
only difference is the source of the internal voltage, which 
develops the field flux wave that interacts with the arma¬ 
ture flux wave. The field flux wave is a result of per¬ 
manent magnets built into the rotor to develop magnetic 
poles. 

A typical magnetic circuit, combined with Ampere’s Law, 

shows 



Hmlm + Hg g — Ni 


Figure 40. Typical magnetic circuit 


We use a constitutive law of B = PH and assume that any 
steel has P = ®. If flux leakage is about zero, BmAm = BgAg 
since flux is solenoidal. Then, 

Ni - Hgg Po Hg Ag 

lm -- and Am = - 

Hm Bm 

The magnet volume is 

„ , A (Ni - Hgg) Po Hg Ag 

Hm Bm 

Minimum magnet volume occurs when the magnet’s maximum 
energy product (MEP), HmBm, is a maximum. If current is 
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zero and only magnitudes are used, 


Vm 


g_Bg 2 __Ag 
P o Hm Bm 


The load line of the magnet is developed as 


1m 

g 

Bm 


Am 

Po 

Hm Ag 


Bm 

Po 

lm 

Ag 

Hm 

g 

Am 


The permanance coefficient, Pc, is the slope of the load 
line. On a magnet diagram 



Figure 41. Magnet operating point diagram 

A good algorithm for machine design is to adjust the 
geometric dimensions to operate at the MEP, on the load 
line. If operation at MEP is assumed, the needed slope is 
determined, the dimensions are randomly generated, PoHg is 
calculated, and the design is maximized for PoHg and mini¬ 
mized for magnet volume, then a search technique has been 
delineated. 

Magnets may "overhang" the active length at either end 
to account for manufacturing imprecision and to permit a 
smaller armature diameter. This overhang affects the 
developed flux. 
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Figure 42. Permanent magnet machine diagram 

The flux per pole, 0, with overhang and a given Ja and Jm, 
is the same as would exist in a configuration in which 
Ja’ = Jm’ = Ja + ND and there were no overhang. N is a non¬ 
linear parameter promulgated in graph form by Ireland [33]. 

A good power fit for N is 

(S’ Jm 

N = 0.38558 (- )o. 7061 3 

D 

6 Jm 

where - = overhang and 0 i overhang i 0.34. Then, 

D 

Ja SJm 

-+- 

D D 

Y = - and 0wi t h = 0w/o Y 



D 


This flux-with-overhang is applied to the problem as would 
be the usual flux. What is the usual flux? A permanent 
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magnet hysteresis diagram shows the residual flux density to 



The magnet does not operate at that point, but rather 
on the load line. Also, there is not magnet material at 
every point along the circumference of the rotor. A Fourier 
series is a good way to find the flux-without-overhang. 

Using a developed rotor, 



Figure 44. Magnet material on a developed rotor 


PM = Bd 

wr p wr 
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) dx for n = 1 , 2, 3 ... 


2 P T rot x 


An =- 

T 

Jo B(X) 

cos ( - 

) dx for 

T 

B(x) = 

Bd 

for 

0 i x i (tt r 


0 

for 

(tt r lr/p) i 


Bd 

for 

(ttt/p) < x < 


0 

for 

(Trr/p) (1+lr) 


lr/p) 

x < (TTr/p) 
(TTr/p) (1 + lr) 
i x i ( 2 tt r/p) 


After integration and manipulation, 


2Bd tt lr it lr 

An =-[ sin(-) - cos(-) + 1] 

tt 2 2 


The equation represents only the fundamental term of the 
Fourier series. The flux-without-overhang is 


2Bd 

Br =- 

TT 


sin( 


tt lr 



tt lr xp 

cos (-) + 1] cos (-) 

2 2r 


As is usual in steady-state analysis, the magnitude is used. 

The next quantity to find is Bd. If the magnet 
operates on the load line, the operating point flux is 


Brem 

Bd =- where Pr is the relative 

1 + p r /p c 

reversible permeability of the magnet and Brem and Pc are as 
previously defined. With this relative permeability, the 
magnet length is lm = g Pr . The magnetic machine can now be 
specified. An end view with dimensions is given in Chapter 
Four. 
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Table 40. Listing of permanent magnet machine design program 


linclude "stdio.h* 
♦include "def.h" 


ft prograa naae: paa.c for permanent aagnet aachines tl 

long int seed; ft start point for randoa nuaber generator If 

double bE 1013C271, hEinC273, ksCBI; 

ft b is "best" array, h is ■hold" array, ks are winding factors tf 

aainO 

{ 

double designpointO, rnd_walk(), swf(), ke, kv, ka, ainpwr, 
stepsize, randoaO, abs(), freq, rpa; 
int p, iteration, i, j, best, print_out(), loops; 

FILE tfopenO, tfp; 


printf("\nReading input data froa PHH.DAT . . 


fp=fopen( > paa.datV r>) 5 ft input seed for randoa nuabers tf 
fscanf(fp, B Xd\iseed); 
fscanfUp, ■«", ip); 
fscanf(fp, "Xlf", iainpwr); 
ainpwrt=746.0; 
fscanf(fp, "Ilf", ike); 
fscanf(fp, "Xlf", ikv); 
fscanf{fp, "Xlf", iks); 
fscanf(fp, "Xlf B , irpa); 
fclose(fp); 


It input nuaber of pole pairs tf 
ft input aachine power, derived fa ASSET 1/ 
ft convert to watts t / 
ft CERs for Effective Height tf 


ft aachine aax shaft rpa tf 


printf( B \nHow aany loops do you want? "); 
scanf("Zd", iloops); 

printf("\n\nDoing prograa calculations . . .\n"); 

for (i=l; i < 8; i+=2) ft haraonic winding factors tf 

{ 

ksEi3=swf(l); 

} 


freq=rpatp/60.0; ft aax electrical frequency tf 

ft MAIN BODY OF THE PROGRAM tf 


for (i=1; i <= loops; ++i) 

{ 

printf( , \ni=Xd a , i); 
stepsize=0.1; 
iteration^; 

design_point(ainpwr, p, ke, kv, ka, freq); 

ft put stuff in the hold array tf 

while (iteration <= 10) 

{ 

rnd_walk(ainpwr, p, stepsize, ke, kv, ka, freq); 
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It stagger around tl 

It index to best EH of the lot I/ 


best=0; 

for (j=l; j<=10; ++j) 

if <hfj][18] < hCbestlC181) 

best=j; It find the best aachine 1/ 

if (abst<hE01E 181 - hEbestlClBU/(htOlUBl+.OOl) 1 < 0.005) 

It saall iiproveaent in EH tl 
It 0.001 takes care of div by zero tl 

{ 

stepsize/=2.0; 

++iteration; 

) 

else It transfers best to 0 position tl 

{ 

for (j=l; j <= 26; ++j) 

htOlfjl = htbestlCjl; 

} 

> 

for (j=l; j <= 26; ++j) 

bti]tj]=htbest][j]; It keep the best aachine tl 

} 

best=l; 

for (i=1; i <= loops; ++i) 

if (bti]C18i ( bCbestlC18]) 

best=i; It find and keep the best of the best tl 

ainpar/=746.0; It turn back into hp tl 

printout (best, p, ainpar, ke, kv, ka, rpa); It output to disk file tl 

fp=fopen( , paa.dat’, VI; It output seed tl 

fprintf(fp,*Id', seed); 

fprintf(fpj'tnXd 1 , p); 

fprintf(fp,"\nllf*, ainpwr); 

fprintf(fp,'\nllf’, ke); 

fprintf(fp,"\nXlf", kv); 

fprintf(fp,"\nXlf", ka); 

fprintf(fp,“\nXlf*, rpa); 

fdose(fp); 

} 

It END OF MAIN PROGRAM; ALL THAT FOLLOH ARE FUNCTIONS tl 

double design_point(ainpwr, p, ke, kv, ka, freql 

It deteraines a randoa design point tl 

double ainpar, ke, kv, ka, freq; 
int p; 

{ 

double r, js, Is, lr, dcore, ds, g, w, 1, xs, va, ph, ataag, 

pe, i2r, vol, at, effcy, ea, xsl, xs5, xs7, xsal, la, volaag, 
bd, find_lr(), ca, siv, scv, riv, loa, doa, find_siv(), atjron; 
extern double sqrtO, randoal); 
int c=0, d=0; 

while (d != 1) ( 
ahile (c != 1) 
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{ 

r=randoa()JRMAX; 

if <<2IPI*r»freq/p) < HAX.TIP.SPEED) It check rotor tip speed »/ 
break; 

} 


M=2*Pl»freq; 


bd=BREH/(1 + HUR/PC); 

It 

lr=find_lr(bd); 

It 

nhile (c != 1) 

It 


{ 

g=randoa()l(0.1tr - GH1N) 
if <g > 0) 

break; 


} 


+ 


air gap flux, operating point, froi aagnet 
characteristics tl 
find aagnet slot factor 1/ 

gap diaension tl 

6H1N; 


la=gtMUR; 


It radial length of aagnet, see PH-B, 12/15/86 1/ 


ls=randoa 010.5 + 0.25; 
js=randoaOIJSHAX; 
dcore=lbdt(r+la))/IBSATlp); 
ds=randoaO10.91dcore; 


It stator slot factor tl 

It full load stator current density tl 

It back iron depth 1/ 

It slot depth < 90Z of body depth tl 


xsl-kstlltksU]; 

xs5=(kst5ms[5J/25); 

xs7=(ksE71lksE71/49); 

xsal=l5tlstPllds/18); 

xs=(HUIj stSSFt(r+g+la)l(r+g+la)tdstPl*ls 

* (12*(xsl + xs5 + xs7)/(Pltgtp) + xsal))/(12t(r+la)»bdlksm); 

It p.u. synch iapedance tl 


if Ixs ) 3.0) 

continue; It don’t nant xs too big tl 


else 


++d; 


It ay escape hatch — aission coaplete tl 


l=(ai npwrtp)/(2tbdtwtksU]t Ir+la) ISSFSds tPl t (r+g+1 a) t Is t js tPF); 

It active length tl 


va=2tPlt(r+la>tll«tbdtksm»jstSSFt(r+q+la)tdstls/p; It VA rating tl 


siv = find_siv(1, r, g, ds, dcore, Is, la); It stator iron voluae 1/ 
riv = HPltrtr; It rotor iron voluae tl 
scv = 2*Pl»(r+g+la)*ds*ls*(l + 2.3094tPl»(r+g+la)tCP/p); 


cm = scvtDCU; 

loa = 1 + 41fr+g+la); 

doa = 2t(r+g+la+ds+dcore); 

vol = VOLALLI(loaJPlldoaldoa/4); 

Mt = HTALLIIcm + DKBRNESIriv + siv)); 
at_iron = DKriv + siv); 


It stator copper voluae I/ 
It total copper Height tl 
It length-over-all 1/ 

It over-all-diaaeter tl 
It aachine envelope voluae 1/ 

It aachine Height in kg tl 
It iron Height only 1/ 


ph = 31.86225IBETAIfreq»BRltHCltHt_iron/D; 
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It hysteresis loss in Matts, uses iron weight of aachine 1/ 


pe = (106236.9INUIBSATtBSATITlJTltfreqtfreqlMt_iron)/(RH0*D); 

It eddy current loss in Matts, uses iron Height of Machine tl 

vol«ag=2tPllr*ltl«*lr; It aagnet voluae 1/ 

Nt«ag=voloagtDHAG; It aagnet Height I/ 

i2r = i.OtSSFIdstCRHOtPIt (1 + 2.3094IPIS(r+g+la)iCP/p) $jsljst(r+g+laHIs; 

It stator copper loss in natts t/ 
It revised 1-12-87 tl 


effcy=(ainpMr)/(ainpHr + ph + pe + i2r); 


en=Ht + katntaag + kel(l-effcy) + kvt(vol + kalvolaag); 

It Effective Height I/ 


hfOlllHs; 

hl0][6]=dcore; 

ti [OH 113=vol; 

ht01[l61=va; 

h[0][201=Htaag; 

h[0H23]=0.0; 

ht0][26]=doa; 

It this 


h[0][2]=freq; 

h[0][71=ds; 

h[0][121=Ht; 

ht0][17]=effcy; 

hl0][2l]=volaag; 


hlO][3]=H; 

hlOK8]=la; 

hl0][13]=ph; 

htOH18]=eH; 


hlOH4]=r; 

h[0im=ls; 

h[0H14]=pe; 

hL0H19]=l; 


h[011221=xs; 


It this one Bakes overhang=0 1/ h[0H241=bd; 


h[0][51=g; 
hC031103=lr; 
hl0][15]=i2r; 


hlOH25]=loa; 


section just changed all the variables in the "hold" array tl 


return; 

)) 


double rnd_walk(ainpar, p, stepsize, ke, kv, ka, freq) 

/I nalks about design_point 10 tiaes 1/ 
double stepsize, ainpnr, ke, kv, ka, freq; 
int p; 

{ 

double r, js, Is, lr, dcore, ds, g, h, 1, xs, va, ph, pe, i2r, 

vol, Ht, effcy, en, xsl, xs5, xs7, xsal, la, volaag, ntaag, bd, 
find_lr(), find.brO, br, overhang, lrat, cm, siv, scv, riv, loa, 
doa, find_siv(), nt.iron; 
extern double sqrtO, randoaO; 
int i=l, j, check=0, ccj; 

while (i <= 10) 

{ 

It read in the Halk around the design point tl 

js=h[0imt(l + stepsizetfrandoaO - 0.5)); 
if (js > JSHAX) 

js = JSHAX; It reset to liait tl 

H=2tPltfreq; 

r=h[0][4U(l + stepsizet(randoa() - 0.5)); 
if ((H»r/p) > HAX_TIP_SPEED) 
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continue; /I go to next try if violated 1/ 

g=h 10]C5]t11 + stepsizeltrandoaO - 0.51); 
if (g < 6HIN) 

g=6MN; It reset to the liait 1/ 


overhang=randoa()10.34; 


l=h[0][19]; 
la=gtHUR; 

bd=BREH/(l + HUR/PC); 
lr=find_lr(bd)j 
1s=h[0][9]I(1 + stepsizellrandoaO - 0.51); 
if (Is > 0.75) 

15=0.75; 
if (Is < 0.25) 

ls=0.25; 


It aagic number "0.34* fro* the book 
on permanent aagnets by Jaaes Ireland. 

The effect of overhang is to concentrate air 
gap flux, or reduce leakage. 1/ 

It nea la, based on the nea g 1/ 

It aagnet operating point flux density tl 
It find nea lr tl 


It reset to the liait 1/ 


for (j=l; j <= 10; ++j) It aagnet char convergence loop 1/ 

{ 

lrat=l/(2llr+g+la>); 
br=bd$find_br(overhang, lrat); 

It find effect of overhang on air gap flux tl 
dcore=(brl(r+la))/(BSATIp); It aost efficient use of iron tl 
ds=h10117]Idcore/hC0116]; It keep thea in the saae ratio 

as upon exit fa design pt 1/ 


xsl=ks[l]*kst 11; 
xs5=(ks[5]»ks[5]/25); 
xs7=(ks[7]*ksC7]/49); 
xsal=(5llstPItds/18>; 

xs=(fHJtjslSSFI(r+g+lB)<(r+g+laUdslPItls I (12<(xsl ♦ xs5 ♦ xs7)/(PIIgtp> 

+ xsal))/(12l(r+la)tbrtkst1]>; 

It p.u. synch iapedance 1/ 

l=(ainparlp)/(2lbrlalks[l]t(r+la)»SSFIds*PII(r+g+la)llslj5<PF); 

It active length tl 

} It end of aagnet char convergence loop tl 

if (xs > 3.0) 

{ 

++check; 

if (check ) 25) It can’t close on decent xs tl 

{ 

for (ccj=i; ccj <= 10; ++ccj) 
htccj][181=10000000.0; 

It aake this randoa aalk undesirable tl 
printfl" burp*); 

break; It go to next design point tl 

} 

continue; It can’t have xs too big 1/ 

} 
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va=2tPIt(r+U>mwtbr»!csm«jstSSFt<r+g+l*)tdsU5/p; 


It VA rating »/ 


siv = find_siv(l, r, g, ds, dcore, Is, la); It stator iron voluae 1/ 
riv = ltPItrlr; It rotor iron voluae 1/ 
scv = 2IPIt(r+9+1a)Jdsllst <1 + 2.3094tPI*(r+g+ia)tCP/p); 

It stator copper voluae 1/ 


ch = scvtDCU; 

loa = 1 + 4tlr+g+l»); 

doa = 2*(r+g+la+ds+dcore); 

vol = VOLALLttloatPItdoatdoa/4); 


It total copper Height t / 
It length-over-all 1/ 

It over-all-dia«eter 1/ 

It aachine envelope voluae 1/ 


nt = HTALLKch + Dt(BRN6Slriv + siv)); 
nt_iron = DMriv + siv); 


It aachine Height in kg 1/ 
It iron Height only 1/ 


ph = 31.86225tBETA* freqlBRltHCltnt_iron/D; 

It hysteresis loss in watts, uses iron Height of aachine I/ 


pe = (106236.9tMU»BSAT tBSATtTlIT1llreqtfreqt»t_iron)/(RHOtB); 

It eddy current loss in natts, uses iron Height of aachine 1/ 

volaag=2*PltrtlI(1.0+(overhang*(2t(r+g+la))))tlatlr; It aagnet voluae t/ 

Htaag=volaagtDHA6; It aagnet Height t/ 

i2r = 2.0tSSFtdslCRH0IPII (l + 2.3094IPIt<r+g+la)*CP/p) *jstjs*(r+g+la)*ls; 

It stator copper loss in natts 1/ 

It revised 1-12-87 »/ 


effcy=(ainpHr)/(ainpwr + ph + pe + i2r); 

en=nt + katntaag + ket(l-effcy) + kvMvol + katvolaag); 

It Effective Height 1/ 


hli Hl]=js; 
hli l[61=dcore; 
h[i][11]=vol; 
htiH16J=va; 


htiJ[21=freq; 
hfi 1E7 J=ds; 
htiJ tI21=nt; 
htiJtl71=effcy; 


h[ill20]=Htaag; hfi H21)=volaag; 
hti3 t231=overhang; 


hti ][3 J=h; 
hti ][81=la; 
hfiIt131=ph; 
hfi1f18]=en; 


hfilf4]=r; 
hfill9)=ls; 
hfi lfl41=pe; 
hti)t191=1; 


hfi][22 )=ks; 


hfi]f51=g; 
hti If10)=lr; 
hfi If151=i2r 


hfi)[24)=br; hfi1t25)=loa; hfi)(26)=doa 

It this section just changed all the variables in the 'hold" array tl 


++i; It go to the next hti HI I/ 

check=0; 

} 


return; 

) 


printout (best, p, ainpnr, ke, kv, ka, rpa) 
int best, p; 

double ainpnr, ke, kv, ka, rpa; 

{ 

char outfiletUl; 
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FILE l-fpo, Ifopen(>; 
int i; 

printf(“\nWhat is the na#e of the file where you want the output? "); 
scanfCXs", outfile); 

fpo=fopen(outfile, V); 

fprintf(fpo,'%d", p); 
fprintf(fpo,“\nXlfainpwr); 
fprintf(fpo,"\nXl-f", ke); 
fprintf(fpo,'\nXlf', kv); 
fprintf(fpo,’\nXlf", ki); 
fprintf(fpo,’\nXlf', rp«); 
for (i=t; i <= 26; ++i) 

fprintf (fpo, ’\nXlf', blbestHi)); 
fprintf(fpo, ’\n"); 
fcloselfpo); 

} 


double find_lr(br) It bracketing il 

double br; 

{ 

extern double cos(), sin(); 
int c=0; 

double 1 r=0.5, topM.O, tbr=0.0, abs(), bt»=0.0; 

while (c!=l) 

{ 

tbr=(2»BREHI (sin(O.SIPIllr) - cos(0.6»PI»lr) + 1.01I/P1; 
if (abs((br-tbr)/br) <= 0.001) /I check for convergence 1/ 
break; 
if (br > tbr) 

( 

bt»=lr; 

lr-bti + (top- bta)/2.0; 

} 

else It br < tbr »/ 

( 

top=lr; 

lr=bt« + (too-bt»)/2.0; 

} 

if (abs((top-1r)/top) <= 0.0005) /I no convergence M 
{ 

printf(’\nXlf, lr); 
abort(‘no solution for lr’); 

} 

} 

return(lr); 

) 
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double find_br(overhang, lrat) 
double overhang, lrat; 

{ 

extern double powO; 

double expon=0.706133501, factor=0.385576838, eon, yyy; 

/I enn is Ireland's N, yyy is Ireland's Y \l 

enn = factor I powloverhang, expan); 
yyy = (lrat + overhang)/(1rat + enn); 
return(yyy); 

) 

double find_siv(l, r, g, ds, dcore, Is, li) 

double 1, r, g, ds, dcore, ds, Is, li; 

{ 

double one, two, three, four; 

one - (r+g+li+ds+dcore)t(r+g+l«+ds+dcore) - (r+g+l«)I(r+g+l«); 
two = 2tPlt(r+g+l«)tdstls; 

three = (r+g+li+ds+dcore)t(r*g+U+ds*dcore) - (r+g+l«+ds)l(r+g+l«+ds); 

four = II(PI lone - two) + PII4l(r+g+l*)lthree; 
return(four); 

) 
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Table 41. Listing of permanent magnet efficiency program 


iinclude "stdio.h 8 
tindude "def.h 8 


It program name: peff.c to find efficiency of permanent magnet machines tl 
it works with a single machine tl 


iain() 

{ 

FILE Ifopen(), Ifp; 

double r, js, Is, Ir, dcore, ds, dr, g, w, 1, xs, li, kc, sgeff=1.0, 

va, ph, pe, i2r, vol, wt, effcy, ew, xsl, xs5, xs7, xsal, ks[83, 
siv, nv, wtjron, find_siv(), pmeff, pgeff, teff, pgjs, pmjs, 


pmrp®, geff=0.98, pceff=0.99, gr, dhp, rpa 
wtfflag, volmag, overhang, bd; 
extern double swf(), sqrtO; 
int e=0, f, p, i; 
char infiieU43; 


iinpwr, ke, kv, freq, 


siv = find_siv(l, r, g, ds, dcore, Is, 1 b); 
riv = llPIlrtr; 
wt_iron = DXriv + siv); 


It stator iron volume tl 
It rotor iron volume tl 
It iron weight only tl 


printf(“\nCalculates efficiency for a single machine.\n H ); 
while (e != 1) 

{ 

f = 0; 

printf("What is the name of the input file? "); 
scanfl^s 8 , infile); 
fp = f open (infile, V); 

fscanf(fp, H, £d B , &p); It input number of pole pairs tl 

fscanf(fp, "XIf", fcminpwr); 
minpwr t= 746.0; It now in watts tl 

fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 


fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 


fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 
fscanf(fp, 


■xif, 

ike); 

■nr, 

(kv); 

■xif, 

Ike); 

"Xlf", 

(rpe); 

"Xlf", 

(js>; 

■Ilf*, 

req); 

■Xlf, 

(w); 

“Xlf*, 

(r); 

■Xlf, 

lg); 

■Xlf, 

(dcore) 

■Xlf, 

Ids); 

■Xlf, 

Me); 

■Xlf, 

(Is); 

"Xlf", 

Ur); 

"Xlf", 

(vol); 

"Xlf", 

(wt); 

"Xlf", 

(ph); 

■Xlf, 

(pe); 

"Xlf", 

(i2r); 

"Xlf", 

(va); 
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fscanf(fp, ’Xlf, ieffcy); 
fscanf(fp, “Xlf*, iew); 
fscanf (fp, ’Xlf’, tl); 
fscanf (fp, “Xlf, iwtsag); 
fscanf(fp, "Xlf*, ivolaag); 
fscanf (fp, “Xlf, Ixs); 
fscanf(fp, “Xlf", ioverhang); 
fscanf (fp, “Xlf, &bd); 
fdoselfp); 

while (f != 1) 

{ 

printf C\nHhat is the sustained speed aachine horsepower? '); 
scanffXlf, idhp); 
dhp »= 746.0; 

printffNhat is the sustained speed »achine rpa? '); 
scanffXlf*, iparpa); 

freq = parpalp/60.0; /I aax electrical frequency 1/ 

p«js = jsldhp/ainpwr; It PM stator current 1/ 

ph = 31.86225IBETA*freq»BRUHCUwt.iron/D; 

It hysteresis loss in watts, uses iron weight of aachine tl 
pe = (106236. ?INUIBSAT*BSATITltTUfreqlfreq»wt.iron)/(RHOID); 

It eddy current loss in watts, uses iron weight of aachine t/ 
i2r = 2.0*SSF»ds»CRH0*Pll (1 + 2.3094IPI*(r+g+laUCP/pl ljs»js»(r+g+la)»ls; 

It stator copper loss in watts I/ 

paeff - dhp/(dhp + ph + pe + i2r); 

printf(“\n Sustained efficiency is Xlf, paeff); 

printf C\n\nHhat is the endurance speed aachine horsepower? •); 
scanffXlf*, idhp); 
dhp »= 746.0; 

printf(“What is the endurance speed aachine rpa? *); 
scanffXlf’, iperpa); 

pajs = jsldhp/ainpwr; It PM stator current 1/ 

freq = parpalp/60.0; It aax electrical frequency 1/ 

ph = 31.862251BETAtfreqlBRUHCllwt_iron/D; 

pe = <106236.91NUlBSATtBSATlTUTllfreqlfreqlwt_iron)/(RH01D); 

i2r = 2.01SSFldslCRH01Pll (1 + 2.3094IPI*(r+g+r«>*CP/p) ljsljsl(r+g+la)tls; 

paeff = dhp/(dhp + ph + pe + i2r); 

printf(*\n Endurance efficiency is Xlf, paeff); 

printff\nSaae aachine? *); 

5 canffXd\ if); 
if (f == 0) 

continue; 
else if (f == 2) 

C 

e = l; 
break; 

} 

} It end of f'loop 1/ 

> It end of e-loop 1/ 

} It end of aain prograa tl 
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double find_siv(l, r, g, ds, dcore, Is, la) 
double 1, r, g, ds, dcore, Is, la; 

{ 

double one, two, three, four; 

one - (r+g+la+ds+dcore) Mr+gHa+ds+dcore) - (r+g+la)*(r+g*la); 
two = 2tPU(r+g+la)Sdstls; 

three = (r+g+li+ds+dcore)t(r+g+li+ds+dcore) - (r+g+la+dsllir+g+la+ds); 

four = lKPIJone - two) + PU4l(r+g+l«)lthree; 
return(four); 

} 
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Appendix D. Induction machines 


The equivalent circuit for an induction machine is 
shown in the figure. 


Ri X, 



Figure 45. Induction machine equivalent circuit 

The stator leakage reactance, XI, is the sum of the 
stator differential leakage reactance (Xsd), the stator slot 
leakage reactance (Xss), and the stator end turn leakage 
reactance (Xse). In turn, Xsd is the sum of the belt and 
zig-zag reactances. 

The rotor leakage reactance, X2, is the sum of the 
rotor differential leakage reactance (Xrd), the rotor slot 
leakage reactance (Xrs), and the skew leakage (Xskew). Xm 
is the magnetizing reactance and Rc is the core resistance. 
Rl and R2 are the stator and rotor resistances, respec¬ 
tively. 

A derivation of the properties of electric machines, 
using Ampere’s Law and the constitutive relation earlier 
postulated, yields the result that 

® Po 6 Ns 2 1 r 

flux = E I --ks 2 cos(wt) 

n = 6k±l g tt n 2 p 2 
k=0 


173 









































From this important stator inductances are taken. 

Belt leakage reactance (Xbelt) is the sum of the reac¬ 
tances due to phase-belt harmonics of an "infinite" slot 
winding. In most machines, the most important harmonics 
present are the fifth and seventh, as the third is canceled 
in balanced operation. Then the belt reactance is 


6PowNs 2 lr ks5 ks7 

Xbelt =- [(-) 2 + (-) 2 ] 

TT g p 2 5 7 


This is the harmonic form of the fundamental mutual reac¬ 
tance of Appendix B. The winding factor is again the 
product of the pitch and breadth winding factors. 

Zig-zag reactance is leakage due to all the air gap 
harmonics that would be produced if the winding had one slot 
per pole per phase. For a phase belt of one slot, with each 
slot carrying the same current and equally separated in time 
and space phase, the zig-zag reactance alone would be 
present. Belt leakage occurs because phase belts are ac¬ 
tually several slots wide. Zig-zag reactance has harmonic 
orders higher than seven, with the same form as Xbelt. No 
even or triplen harmonics will be present. 

The fundamental harmonic of the flux yields the mag¬ 
netizing reactance, Xs, which can be viewed as that required 
to "energize" the air gap. 

6 Po w Ns Irks 2 

Xm =- 

g P 2 

Figure 46 shows a typical stator slot. The stator slot 
leakage reactance, summing the self and mutual reactances, 
is 


Xss = 


18 w Po 1 Ns 2 d2 
ns 2 ( w2 


ds ns v , v 

+ --+ PN P ) 

2 ws 6 
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where ns is the number of stator slots and Np is the coil 
throw, in slots. Np = (CP ns)/(6 p) , a result available by- 
manipulation. w2 



ws 


Figure 46. Stator slot diagram 


Then, 


3 w Mo 1 Ns 2 d2 ds 

Xss =- + ---)(1 + CP) 

ns w2 2 ws 

and 

(1 - Is) 2 t r (r+g) 

ns =- 

ws 

where we have traded the need for the knowledge of the num¬ 
ber of stator slots for the need to know the width of an in¬ 
dividual slot. A reasonable relationship between slot 
dimensions is d2 = ds and w2 = ws. Then, 


Is 2tt (r+g) 

The stator end turn leakage reactance may be estimated 
by treating the two end regions as a single helically shaped 
winding. If the active region of the machine is ignored and 
the helix given air core properties, the inductance can be 
found from standard sources. 


Xse 


-16 w tt Mo r 2 sin 2 (9w/2) No 2 
e« 2 le p 2 


Ip’ ( 


pn r 
le 


)Kp’ ( 


prr r 
le 


) 


where le is the combined length of both end windings. 


175 































Presume the winding radius to be r, with helix pitch 
6w = tt/3 and 

le = 2Trr tan(Tr/3) Np /ns 

Ip’ and Kp’ are the first derivatives with respect to their 
arguments of the hyperbolic Bessel functions Ip and Kp. 

When the three phases are summed, a multiplier of 1.5 will 
be realized. Finally, using the previous results for Np, 


Xse 


81 wPoNs 2 r pirr pirr 

- Ip > ( - )Kp .( - 

pH 2 CP tan(TT/3) le le 


Stator resistance is R1 = (rho It Ns)/(As/Ns), as in 
Appendix B, leading to 


R1 - cp/p) 

SSF ds it (r+g) Is 

When the actual calculation is performed, 1 is not known. A 
guess-and-iterate scheme is used. Iteration continues until 
convergence on 1 is achieved. 

Rotor resistance uses a similar scheme, but the 
presence of rotor bars and end rings instead of turns 
changes it somewhat. A model of a rotor bar is below. 


dr 

v 

wr 



Figure 47. Typical rotor bar configuration 


Induction motor transformer models provide a way to 
find rotor resistances and inductances. The flux density 
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produced by the stator and rotor is 


Mo 

B =-(Fs + Fr ) 

g 

The stator mmf is Fs = -j 


3 Ns ks ws 
tt p 


Ti , where Ti is stator 


current and ws is slip frequency. Rotor mmf is 
ws r 2 Br 

Fr = -j- 

p 2 Zs 

where Zs is the rotor surface impedance. The air gap volt¬ 
age, or voltage across Xm, is defined in terms of the flux 
density and rotor mmf as 


Vag = -j 


w 

Ws 


2 p 1 Ns ks 


Fr 


If rotor mmf is now identified with rotor current referred 
to the stator winding, 


3 Ns ks 


Fr = -j - Iz , where 12 is rotor current. Then, 

tt p 

Zs w 6 1 Ns 2 ks 2 

Vag = ---I 2 

ws tt r 


Separating Zs into its real and reactive parts and using a 
rotor surface model to describe the relation between rotor 
electric field amplitude and rotor surface current yields 


12 1 Ns 2 ks 2 

R2 = -rs lot 

nr 

12 1 Ns 2 ks 2 

Xrs =-xs lot 

nr 

6MowlNs 2 rks 2 1 1 _ , 

Xrd =-I {(-) 2 + (---- -) 2 } 

p+k*nr p-k*nr 


tt g 
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If magnetic diffusion is ignored, end ring resistance 
can be calculated by comparing losses in the rings and slot. 
The ratio of current densities is found by the ratio of the 
areas. This is then squared and multiplied by the ratio of 
volumes. When summed, 

12 1 Ns 2 ks 2 nr r wr 

R2 =-rslot [1 +-] 

nr tt 1 ler p 2 

rslot = CRHO/(dr wr), nr is the number of rotor bars 
and ler is the end ring length, approximated as 
ler = 2 tt (r - wr/4 - ds/2). The rotor bar width, wr, is 
found by specifying nr and observing that (nr wr) = ( 2 tt lr) . 

Rotor skew leakage arises when the rotor slots are 
skewed angularly along the axial length to prevent rotor 
cogging. Then, flux does not fully link the bars. When the 
effect is integrated over the rotor, it is seen that 

Xskew = Xm [1 - (2 sin(skew/2)/skew) 2 ] with the amount 
of skew measured in radians. When typical values of skew 
are input, we see that (Xskew/Xm) ~ 0.5%. This is a negli¬ 
gible effect and will be ignored. 

From the previous rotor bar model, it is seen that 
xslot = w Po (d2/w2 + dr/2wr). Assume that d2 = wr/4 
and w2 = wr/4. Then xslot =wPo (1+ dr/2wr). 

Fitzgerald et al [23] state that only small errors 
result if Rc is omitted. Therefore, the core branch may be 
omitted. 


Once the components of the equivalent circuit have been 
calculated, the designer must turn to power and torque con¬ 
siderations. The internal mechanical power of the machine 
is 


P 


(1 - flip) 

slip 


3 12 2 R2 


The air gap voltage has been previously defined. The ter¬ 
minal voltage, Vt, may be found from Vag by means of a volt- 
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age divider and a Thevenin equivalent circuit developed. 



Figure 48. Induction machine Thevenin equivalent circuit 


R1 + j(Xl + Xm) 

Rel + jXel = (R1 + jXl) in parallel with jXm. Then operat¬ 
ing point torque is 

1 3 Via 2 (R2/slip) 

T =-- 

w (Rel + R2/slip) 2 + (Xel + X2) 2 

Torque is a maximum when the power delivered to (R2/slip) is 
a maximum. By matching load and Thevenin impedances, the 
power is a maximum and a slip-at-maximum-torque is found 


_R2_ 

Sma * ” (Rel 2 + (Xel + X2) 2 )0.5 

and the corresponding torque is 

1 1.5 Via 2 

Tin ax ________ 

w Rel + (Rel 2 + (Xel + X2) 2 )0-5 

Typical induction motors have the ratio between maximum and 
operating point torque as 

T 

- ~ 0.55 

Tma x 


179 












































which is used to find the machine active length to use in 
the circuit component calculations. The equation for rated 
torque can be manipulated to yield a quadratic expression 
for operating point slip, or a Newton’s method convergence 
can be used to find operating point, slip. Convergence on 
slip and active length through Newton’s method is used to 
generate the machines of this thesis. 

Finally, rotor copper losses are (1 - slip)P; stator 
copper losses are found using Ii and R1. 
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Table 42. Listing of induction machine design program 


(include “stdio.h' 
(include “def.h* 


/I prograa naae: ind.c for induction aachines, 4/19/87 tl 

long int seed; It start point for randoa nuaber generator I/ 

double bl261[331, htlll[33], xs[42); 

It b is "best" array, h is “hold" array, xs are winding factors 1/ 

aaint) 

{ 

double designjointt), rnd_walk(), swfl), ke, kv, ainpwr, 
stepsize, randonl), absl), freq, rpa, ksC431; 
int p, iteration, i, j, best, print,out(), loops, flag; 

FILE IfopenO, (fp; 

pnntf C\nReading input data froa IND.DAT . . .*); 


fp-fopenl’ind.dat“,“r“); 
f scant (fp.'ZdVseed); 
fscanftfp, ’ld‘, ip); 
fscant(fp, “Ilf", iainpwr); 
ainpwrl=746.0; 
fscanftfp, “Ilf, fcke); 
fscanftfp, •Ilf", &kv>; 
fscanftfp, ■Jlf, Lrpa); 
fclosetfp); 


It input seed for randoa nuabers tl 

It input nuaber of pole pairs tl 

It input aachine power, derived fa ASSET tl 

It convert to watts tl 

It CERs for Effective Height tl 

It aachine aax shaft rpa (/ 


printf("\nHow aany loops do you want? “); 
scantt-Zd-, Lloops); 

printf(”\n\nDoing prograa calculations . . .\n“); 


for ti=l; i < 42; i+=2) It haraonic winding factors tl 

{ 

kstihswf (i); 

xsfil = Ikstiltkstill/titi); 

} 

freq=rpatp/60.0; It aax electrical frequency tl 


It MAIN BODY OF THE PROGRAM tl 


for (i=1; i <= loops; ++i) 

{ 

printft’\nld“, i); 
stepsi 2 e= 0 .1; 
iteration=0; 
flag = 0; 

designpointtainpwr, p, ke, kv, freq, fcflag); 

It put stuff in the hold array tl 

if (flag == 1) 

{ 

hlOHlSl = 10000000.0; 
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best = 0; 
} 


hhile ((iteration <= 10) tt (flag != 11) 

I 


rnd_walk(ainpnr, p, stepsize, ke, kv, freq); 

It stagger around tl 

best=0; It index to best EH of the lot 1/ 

for (j=l; j<=10; ) 

( 

if (h(j]U8] < 1.0) 
continue; 

if (htj]E181 < hEbest1C 181) 

best=j; It find the best lachine tl 

} 

if (absE(hC01[181 - hEbestl£181)/h£01[181) < 0.005) 

It saall iaproveaent in EH tl 

{ 

stepsize/=2.0; 

++iteration; 

) 

else It transfers best to 0 position tl 

{ 

for (j=l; j <= 32; ++j) 

hEOlCj] = hEbest1E j1; 

} 


} 


for (j=l; j<= 32; ++j) 

bEi1Ejl=h[best 1Ej1; It keep the best aachine 1/ 

} 


best=l; 

for (i=1; i <= loops; *+il 
{ 

if (bEi1E1B1 < 1.0) 


continue; 

if (bEiIE 181 < bEbestlElBl) 


} 


best=i; 


■i npwr/=746.0; 

print_out(best, p, ainptir, 


It find and keep the best of the best »/ 
It turn back into hp tl 

ke, kv, rpa); It output to disk file tl 


fp=fopen(’ind.dat‘,V>i /< output seed tl 

fprintf(fp,*Zd*, seed); 

fprintf(fp,"\nld“, pi; 

fprintf(fp,*\nllf", ainpnr); 

fprintf(fp,“\nllf", ke); 

fprintf(fp,*\nXlf", kv); 

fprintf(fp,*\nXlf*, rpa); 

fdose(fp); 

} 

It END OF MAIN PR06RAH; ALL THAT F0LL0H ARE FUNCTIONS tl 


double design_point(ainpHr, p, ke, kv, freq, fla) 


182 

























































































/I deternnes a randoi design point I/ 


double tfiinpwr, ke, kv, freq; 
int p, tf1 a; 

{ 


double r, Is, lr, dcore, ds, dr, wr, g, w, ha, 1, i2rr, 
rl, r2, xbelt, xzz, xss, xse, xi, xrd, xrs, 
xa, vag, vt, via, rei, xei, taax, radical, sb, ila, lib, 
x2, va, ph, pe, i2r, voi, wt, effcy, ew, le, tieinpwr, 
sa, sc, ch, siv, scv, riv, rev, loa, doa, r2siax, ttaax; 
double ler, fl, f2, f3, slip, saax, jr, js, 
il, i2, xslot, wt_iron; 

extern double sqrtl), randoal), cosl), sin(), cosh(), sinhO, vratl), 
floor!), besipl), beskpO, tanl), find_siv(), smH); 
int cslip=0, nr=71, ccr=l, ccg, ccs=l, k; 

w=2tPIIfreq; It synchronous frequency in rad per sec tl 

w«=w/p; It aech angular velocity in rad per sec 1/ 

1 = 10.0; It initial guess on length tl 

taax = ainpwr/(walPSI); It pull-out torque tl 


while (cslip != 1) 
{ 


r=randoa(>IRMAX; 

if ((2iPIIrlfreq/p) < HAX_T1P_SPEED) It check rotor tip speed tl 
break; 

} 


lr=r andoa 010.5 + 0.25; 
dr=randoa()tr/3.0; 
dcore=(BRIr>/(BSATIp); 
ds=randoa()10.9ldcore; 
g=randot01(0. Hr - GRIN) 
if (g < SHIN) 

g = SHIN; 

ls=randoa()10.5 + 0.25; 


It rotor slot factor tl 

It slot no deeper than 33X of rotor radius tl 

It back iron depth tl 

It slot depth < 90Z of body depth to start tl 
+ SHIN; 


It stator slot factor tl 


f 1=0. 
f 2=0- 
f 3=0. 


for (k=i; k < 100; ++k) 
{ 


fi = (i./(p + ktnr)) t (i./(p + ktnr)); 
f2 = (i./(p - klnr)) t (l./(p - klnr)); 
f3 += fl + f2; 

} 

while (ccr <= 5) It start convergence loop on r tl 

{ 

if (ccs > 5) It we've had trouble with jr tl 

{ 


r t= jr/JSMAX; 

if (r > (MAXJlP_SPEEDtp/w)) 
r = KAX_TIP_SPEEDSp/w; 
if (r > RHAX) 


r = RHAX; 
lr 1= jr/JSHAX; 
if (lr > 0.75) 

lr = 0.75; 
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} 


wr 2 (BtPItlr*(r - dr/2.))/(4lnr + Plllr); It rotor slot width tl 
ler = 2tPll(r - wr/4.0 - dr/2.0); It end ring length 1/ 

ccs = 1; 

while (ccs <= 5) It start convergence loop on ds and dr tl 

{ 

ccg 2 1; 

while (ccg <= 10) It start convergence loop on length tl 

{ 

rl 2 (6ICRH0I(1 + 2.3094tPlllr+g)»CP/p)>/(SSFIdslPI»lr+g)»Is); 
r2 2 (6txst11ICRH0II*(1 ♦ (2trtrlIr)/(Iertptp)ll/(Pllrllrldr); 
xbelt 2 (ilwtllt(Ulrt(xst5] + xst7])/(Pllgtplp)); 
xzz = (itwtllMUIrKxstll] + xstUl + xsC17] + xstl9] + xst23] + 
xsl291 + xs[31] + xs1373 + xs[41]l)/(Pllglptp); 
xss 2 (3tHUtwll»dsl(l+CP))/(Isl2tPII(r+g)); 

Ie = (PllrICPttan(PI/3))/(3lp); 

xse = (27lwlIII1Ulr| r «(-beskp(p, plPIIr/Ie)tbesip(p, ptPIlr/le)))/ 
(Pltletptp); 
xl = xbelt + xzz + xss + xse; 
xrd = (6IHUIwlIlxsmtrlf3)/(PI«g); 

xrs = (itxsUltwtltNUI(wr+0.5ldr) 1/(IrlPIl(r - (dr/2.0) - (wr/8.0))l 
x2 = xrd + xrs; 

xa = (6tHUtwlltrtxsU))/(Pltgtptp); 

vag 2 2trfBRtwtlfsw-f(1)/p; It air gap voltage tl 

rel = (x>txa$rl)/(rltrl + (xl+xa)l(xl+xa)); 

It thev equiv resistance tl 

xel = (xaltrltrl + xltxl + xllxa)l/(rllrl + (xl+xa)t(xl+x«)); 

It thev equiv inductance tl 

r2saax 2 sqrtlreltrel + (xel+x2lt(xel+x2)l; It r2 at saax tl 
vt 2 vag/vrat(r2saax, x2, xl, xa, rl); 

It tersinal voltage at saax tl 
via 2 (vtlxa)/sqrt((xl + xa)l(xl + xa) + (rltrl)); 

It thevenin equivalent voltage tl 
ttaax = (1.5lvlatvla)/(wal(rel+sqrt(reUrel+(xel+x2)l(xel+x2)))); 

It test aaxiaua torque tl 
if (abs((teax-tteax)/taax) <= 0.005) 

break; It we have convergence tl 

1 I 2 taax/ttaax; It reset I, reiterate 1/ 

++ccg; 

if (ccg > 10) 

{ 

tfla 2 1; 

printf(*\n flag set on length’); 
return; 

) 

) It end convergence loop on length tl 

saax = r2/sqrt(rellrel + (xel+x2)l(xel+x2)i; 

slip 2 saax/3.0; It starting point for slip converge tl 

cslip 2 1; 

while (cslip <= 20) 

{ It start convergence loop on slip tl 
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tainpwr = <3tvlalvla*r2/slip)/( (rel+r2/slip)t (rel+r2/slip) + 
(xel+x2)t(xel+x2)); 

if (abs((ainpwr-tainpwrl/ainpwr) <= 0.005) 

break; It we have convergence tl 

slip 1= tainpwr/ainpwr; 

++cslip; 

) It end convergence loop on slip tl 


i2 = vagtsqrt((r2/slip)t(r2/slip) + x2tx2)/< <r2/slip)t(r2/slip) + *21x21; 
jr = <3ti2)/(PItlrtdrtrIRSF); It rotor current density tl 

if (jr > JSHAX) 

{ 

dr jr/JSHAX; 
if (dr > (r/31) 

dr = r/3; It reset to liait tl 

} 

ila = (r2/slip)t(xa+x2) - x2lr2/slip; 
ilb = *2t(xa+x2) + (r2/slip)l(r2/slip); 

il = vagtsqrt(ila»ila + i lbti Ib)/(x«l(x2lx2 + (r2/slip)Kr2/slip))); 
js = (3.0lil)/lPll(r+g)ldsllslSSF); 
if (js > JSHAX) 

ds 1= js/JSHAX; 

if ((js <= JSHAX) W (jr <= JSMflXl) 
break; 

++ccs; 

} It end convergence loop on ds and dr tl 

if (ccs <= 5) 
break; 

++ccr; 

if (ccr > 5) 

/ 

l 

ifla = 1; 

printfCkn flag set on r’l; 
return; 

) 

} It end convergence loop on r tl 

It calculations tl 


saa* = r2/sqrt(rellrel + (xel+x2)l(xel+x2)); 


siv = find_siv(1, r, g, ds, dcore, Is); It 
riv = MPUrMr - 2ldrtlr); It 
scv = 2XPlt(r+g)»ds»ls»(l + 2.3094*Plt(r+g)ICP/p1; 

It 

rev = wrtdrMnrll + 211 er 1; It 
cw = (rev + scvlIDCU; It 
loa = 1 + 4t(r+g); It 
doa = 2l(r+g+ds+dcore); It 
vol = V0LALLI(loatPItdoatdoa/4); It 
wt = HTALLKcw + DKBRNGSIriv + sivl); It 
wt_iron = DKriv + siv); It 
ph~= 31.8A225tBETAtf reqlBRlIHCltwt_iron/D; It 
aachine tl 


stator iron voluae tl 
rotor iron voluae 1/ 

stator copper voluae tl 
rotor copper voluae tl 
total copper weight tl 
length-over-all tl 
over-all-diaaeter 1/ 
aachine envelope voluae 1/ 
aachine weight in kg tl 
iron weight only tl 

hysteresis loss in watts, uses iron weight of 


pe = (106236.9«NU*8SATtBSATtTlIT Ilfreqtfreqtwt_iron)/(RH0*D); 
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It eddy current loss in Matts, uses iron weight of aachine tl 
i2r=3.0tiHiUrl; It stator copper loss in watts tl 

i2rr = sliplainpwr; It rotor copper loss in watts tl 

vt = sqrt((vag+i 1*r 1)*(vag+i 1 >r 1) + (iltxDMiltil)); It rated voltage tl 
va = 3*vt>i 1; It VA rating at terainals tl 

effcy=(einpwr)/(ainpwr + ph + pe + i2r ♦ i2rr); ew=wt ♦ ket(l-effcy) + kvtvol; 
fecfive weight tl 


hfO]Ul=il; hC0H21=freq; ht0][3]=w; 
hC0]I6]=dcore; hC0][7]=ds; h[0KB]=dr; 
htO] 111 ]=vol; htOH121=wti h[0][131=ph; 
htOH16]=va; ht0]tl71=effcy; h[Q][18]=ew; 
h[0H211=saax; htOH221=taax; hI03C231=vt; 
h[0K26]=rl; hC0][271=xl; h[0K28]=xa; 

htOH31]=loa; hC0K321=doa; 

It this section just changed all the 


h[0H4]=r; 

ht0][91=ls; 

htOH 141=pe; 

ht01tl91=l; 

htOH24]=slipi 

htOH291=x2; 


htOK51=g; 
h[031101=1r; 
ht0][151=i2r; 
htOlC201=i2rr; 
htOH25]=vag; 
h C01t301=r2; 


variables in the 'hold' array tl 


It Ef- 


return; 

> 


double rnd_walk(*inpwr, p, stepsize, ke, kv, freq) 

It walks about design_point 10 tiaes tl 
double stepsize, sinpwr, ke, kv, freq; 
int p; 

{ 

double r, Is, lr, dcore, ds, dr, wr, g, w, wa, 1, i2rr, ila, ilb, 
rl, r2, xbelt, xzz, xss, xse, xl, xrd, xrs, r2saax, ttaax, 
xa, vag, vt, via, rel, xel, taax, radical, sb, tainpwr, 
x2, va, ph, pe, i2r, vol, wt, effcy, ew, 1e, cw, siv, scv, 
sa, sc, riv, rev, loa, doa, jr, js, 
ler, fl, f2, f3, slip, saax, il, 12, xslot, wtjron; 
extern double sqrtO, randoed, cos(), sin(), coshO, sinhO, vratO, 
besipO, beskpO, tan(), find sivO, swf(); 
int d, i=1, ccg, ccr, ccs, nr=71, k, cslip; 
f1=0.; 
f2=0.; 
f3=0.; 

for (k=l; k < 100; ++k) 

{ 

fl = (l./(p + klnr)) I (l./(p + ktnrl); 

<2 = (l./(p - ktnr)) t Cl./Cp - ktnrl); 
f3 += fl + f2; 

) 

while (i <= 101 It ten steps around the design point tl 

{ 

It read in the walk around the design point tl 
w=2IPItfreq; 

r=h[0H41t(1 + stepsizettrandoad - 0.5)1; 
if ((wtr/p) > (1AX.TIP.SPEED) 

continue; It go to next 1-loop if violated tl 

g=h[0K5]t(l + stepsizet(randoa() - 0.5)); 
if (g < SHIN) 
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g=6HIN; /» reset to the lieit */ 

dcore=(BRIr)/(BSATIp); /» *ost efficient use of iron 1/ 

ds=ht0][7]t(1 + stepsizet(randoa() - 0.5)); 

dr=htO][B]l(l <• stepsizel(randoel) - 0.5)); 

if (dr > r/3.0) 

dr = r/3.0; It reset to the li«it 1/ 

ls=h[0][?]»(i + stepsizeKrandoeO - 0.5)); 
if (Is > 0.75) 

ls=0.75; It reset to the lint »/ 

if (Is < 0.25) 

1s=0.25; 

lr=h[0][10]t(l + stepsizetlrandoeO - 0.5)); 
if (lr > 0.75) 

lr=0.75; It reset to the liiit t / 

if (lr < 0.25) 
lr=0.25; 

It »ech angular velocity in rad per sec I/ 
It length starting point tl 
It pull-out torque tl 


**•=**/p; 

1 = ht01C19]; 
tiax = htOl1221; 


ccs = 1; 
ccr = 1; 

while (ccr <= 5) 

{ 

if (ccs > 5) 

{ 

r 1= jr/JSHAX; 

if (r > (f1AX_TlP_SPEED*p/n)) 
r = HAX_TlP_SPEEDtp/w; 
if (r > RHAX) 

r = RHAX; 
lr 1= jr/JSHAX; 
if (lr > 0.75) 

lr = 0.75; 

} 

wr = (BtPltlrKr - dr/2.))/(4tnr + Pltlr); It rotor slot width tl 
ler = 2tPlt(r - wr/4.0 - dr/2.0); It end ring length tl 

ccs = 1; 

while (ccs <= 5) It start convergence loop on ds and drt/ 

{ 

ccg = 1; 

while (ccg <= 10) It start convergence loop on length 1/ 

{ 

rl = (6tCRH0t(l + 2.3094IP1I(r+g)tCP/p))/(SSFtdstPII(r+g)tls); 
r2 = (6»xstl]tCRH0tU(l + (2trtrtlr)/(lerlptp)))/(Pltrtlrtdr); 
xbelt = <A*wlltHUtrt(xs[5] + xst7])/(PUg»p»p)); 
xzz = (6twtl»HU»r»(xstll] + xsU3] + xs[17] + xstl91 + xst231 + 
xsC29) + xst311 + xst371 + xst411)) / (PUgtptp); 
xss = (3tHUtw»llds»(l+CP))/(lsl2»Pll(r+g)); 
le = (Pl»rlCPltan(PI/3))/(3»p); 

xse = (27twmHUtrlrt(-beskp(p, ptPltr/le)»besip(p, ptPItr/le)))/ 
(Pltlelplp); 
xl = xbelt + xzz + xss + xse; 


It first tiae thru, no adj in r tl 

It start convergence loop on r tl 
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xrd = (6IHUMl$xsmirtf3)/(PIl9); 

xrs = (6<X5tl]lHm!1Ut(wr+0.5tdrJ)/(lrIPIt(r - (dr/2.0) - (Mr/8.0))); 
x2 = xrd + xrs; 

xffl = (6IMUIwtl Ir IxeH 3) / (Pltglptp); 

vag = 2lrlBRtwlllswf(l)/p; /) air gap voltage )/ 

rel * (xalxatrl)/(rltrl + (xl+xa)t(xl+xa)); 

ft thev equiv resistance 1/ 

xel = (xiMrllrl + xltxl + xllxa)) / (r 1 trl + (x 1+xa) I (x 1+xa)); 

ft thev equiv inductance 1/ 

r2siax = sqrt(rellrel + (xel+x2)I(xel+x2)); ft r2 at saax 1/ 
vt = vag/vrat(r2siax, x2, xl, xa, r1); 

ft terminal voltage at saax 1/ 
via = (vtlx*)/sqrt((xl + xa)t(xl + xa) + (rltrl)); 

it thevenin equivalent voltage tl 
ttaax = (l.Stvlatvla)/(wat(rel+sqrt(reltrel+(xel+x2)I(xel+x2)))); 

it test eaxiaua torque 1/ 
if (abs(<taax-ttaax»/taax) <= 0.005) 

{ 

ccq = 1; 

break; It we have convergence 1/ 

> 

1 1= t»ax/ttsax; It reset 1, reiterate 1/ 

++ccg; 

) It end convergence loop on length 1/ 

if (ccs > 5) 
break; 

s»ax = r2/sqrt(reltrel + (xel+x2)t(xel+x2)); 
slip = siax/3.0; 

It starting point for slip converge I/ 

cslip = 1; 
while (cslip <= 20) 

£ It start convergence loop on slip tl 

tiinpwr = (3tvlatvlalr2/slip)/((rel+r2/slip)t(rel+r2/slip) + 
(xel+x2)t(xel+x2>); 

if (abs((«inpwr-tainpwr)/«inpwr) <= 0.005) 

break; It we have convergence tl 

slip 1= tiinpwr/unpwr; 

++cslip; It no »ore than 20 tries tl 

} It end convergence loop on slip tl 

i2 = vagtsqrt((r2/slip)»(r2/slip) + x2tx2)/((r2/slip)t(r2/slip) + x2tx2); 
jr = (3ti2)/(PItlrtdrtrtRSF); It rotor current density tl 

if (jr > JSMfiX) 

£ 

dr *= jr/JSNAX; 
if (dr > (r/3)) 

dr = r/3; It reset to liiit tl 

) 

ila = Ir2/slip)t(x»+x2) - x2tr2/slip; 
ilb = x2t(xi+x2) + (r2/slip)t(r2/slip); 

il = vaglsqrt(ilalila + ilbtilb)/(xit(x2tx2 + (r2/slip)t(r2/slip))); 
js = (3.0ti1)/(Pit(r+g)IdstlstSSF); 
if (js > JSttfiX) 
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ds 1= js/JSMAX; 

if (Ijs <= JSflAX) M (jr <= JSI1AX)) 
break; 

++ccs; 

} It end convergence loop on ds and dr I/ 

if (Ices <= 5) !! (ccg > 10)) 
break; 

++ccr; 

) It end convergence loop on r tl 


if ((ccg >10) 1! (ccs >5) I! (cslip >20) !! (ccr > 5)) 

{ It no design convergence tl 

for (d=i; d <=10; t+d) 

h[dKl8] = 10000000.0; 
printf("burp ■); 
break; 

} 


It calculations tl 

saax = r2/sqrt(rel»rel + (xel+x2)!(xel+x2>); 


siv = find_siv(l, r, g, ds, dcore, Is); It 
riv = HPIlrKr - 21dr!lr>; It 
scv = 21PIl(r+g)idsilsi(l + 2.3094!PIi(r+g)iCP/p); 

It 

rev = wrtdrilnrll + 2!ler); It 
cw = (rev + scvlIDCU; It 
loa = 1 + 4!(r+g>; It 
doa = 2*(r+g+ds+dcore); It 
vol = V0LALL!(loalPIldoaidoa/4); It 
wt = UTAH! (cm + DtlBRNBSIriv + sivl); It 
wt_iron = Dt(riv + siv); It 
ph’= 31.86225!BETAifreqtBRllHCltwt_iron/D; It 
cachine tl 


stator iron voluae tl 
rotor iron voluae tl 

stator copper voluae tl 
rotor copper voluae tl 
total copper weight tl 
length-over-all !/ 
over-all-diaaeter tl 
aachine envelope voluae !/ 
aachine weight in kg tl 
iron weight only tl 

hysteresis loss in watts, uses iron weight of 


pe = (106236.9!NU1BSAT!BSAT!T1!T1Ifreqlfreq!wt_iron)/(RH0!D); 

It eddy current loss in watts, uses iron weight of aachine !/ 
i2r = 3.0tiliil!rl; It stator copper loss in watts tl 

i2rr = sliptainpwr; It rotor copper loss in watts tl 

vt = sqrt((vag+il!rl)!(vag+il!rl) + (il!xlH(il!xl)); It rated voltage tl 
va = 3!vttil; It VA rating at the terainals tl 

effcy=(ainpwr)/(ainpwr + ph + pe + i2r + i2rr); ew=wt + ket(l-effcy) + kvtvol; It Ef¬ 

fective weight tl 


hti][l]=il; 

hli][6]=dcore; 

hliIt 113=vol; 

h[i][161=va; 

hti][21]=saaxj 

h[iH26]=rl) 

hti](31]=1oa; 

It this 

++ i; 

) 


h[iH2]=freq; hli]13]=w; 
hti]t7]=ds; hli][B]=dr; 
hliIt12]=wt; hli3113]=ph; 
htiH17]=effcy; hli][10]=ew; 
hti][22]*taax; htilt23]=vt; 
hti][27]=xl; h[iH28]=xa; 

hti 31321=doa; 

section just changed all the 


hli][4]=r; 
hti][9]=ls; 
hliDt14]=pe; 
hti11191=1; 
h[iH241=slip; 
hti11291=x2; 

variables in the 


hti115]=g; 
hli][10]=lr; 
htilt153=i2r; 
hti][203=i2rr; 
hli312S3=vag; 
hti3t303=r2; 

hold" array tl 


return; 

> 
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print out (best, p, rnnpwr, ke, kv, rp») 
int best, p; 

double Binpwr, ke, kv, rpo; 

{ 

char outfilet 143; 

FILE Ifpo, JfopenO; 
int i; 

printfCXnWhat is the na«e of the file where you want the output? °) 
scanf(*^s*, outfile); 

fpo=fopen(outfi1e, V); 


fprintf(fpo,-%d a , p); 
fprintf(fpo, B \nXlf", iinpwr); 
fprintf(fpo,*\nZlf, ke); 
fprintf(fpo,“XnZlf", kv); 
fprintf(fpo,“\nZlf", rpi); 
for (i=l; 1 <= 32; ++i) 

fprintf (fpo,"\n%e“,bibest Hi]); 
fprintf (fpo, 
fdose(fpo); 


double findsiv(l, r, g, ds, dcore, Is) 
double 1, r, g, ds, dcore, Is; 

/ 

i 

double one, two, three, four; 

one = (r+g+ds+dcore)l(r+g+ds+dcore) - (r+g)l(r+g); 
two = 2IPII(r+g)ldslls; 

three = (r+g+ds+dcore)I(r+g+ds+dcore) - (r+g+ds)I(r+g+ds); 

four = lKPIlone - two) + FI 141 (r+g)Ithree; 
return(four); 

} 

double vrat(r2ssax, x2, xl, x«, rl) It ratio of vag/vt tl 
double r2saax, x2, xl, xb, rl; 

{ 

double a, b, c, d, vrat; 
extern double sqrtO; 

a = r2saaxl(xa+xl) + rll(x2+xi); 
b = rllr2s«ax - x«lx2 - xll(x2+xa); 
c = alr2siax - blx2; 
d = blr2saax + alx2; 


vrat = (xalsqrtklc + did)) / (ala + bib); 
return(vrat);} 
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Table 43. Listing of induction e-f-ficiency program 


♦include "stdio.h" 
♦include "def.h" 


It progras naae: ieff.c to find efficiency of induction ®achines tl 

It works with a single aachine tl 

tainl) 

{ 

FILE ifopenl), Ifp; 

double r, Is, lr, dcore, ds, dr, g, w, 1, i2rr, va, ph, pe, i2r, vol, 

wt, effcy, ew, siv, riv, Ht_iron, find_siv(), p®eff, pgeff, teff, 
p®rp®, geff =0.98, pceff=0.99, gr, dhp, rp®, ainpwr, ke, kv, freq, 
sgeff=1.0, wr, ns, rl, r2, x1 f x®, vag, vt, via, rel, xel, radical, 
x2, sa, sb, sc, slip, slipl, slip2, snax, il, i2, taax; 
extern double sqrtl); 
int e-0, f, p, i; 
char infilet 141; 

printf("\nCalculates efficiency for a single ®achine.\n’); 


while (e != 1) 


{ 

f = 0; 

printffWhat is the na®e of the input file? "); 
scanfl’Xs", infile); 
fp = f open (infile, "r"); 

fscanflfp, "Xd", ip); It input nu®ber of pole pairs 1/ 

fscant(fp, *Xlf", i«inpwr); 

unpwr t- 746.0; It now in watts tl 

fscanflfp, "Xlf, ike); 

fscanflfp, •Xlf, ikv); 

fscanflfp, ■Xlf", irpat; 

fscanflfp, "Xlf, ill); 

fscanflfp, "Xlf", ifreq); 

fscanflfp, ’Xlf’, in); 

fscanflfp, “Xlf, ir); 

fscanflfp, “Xlf*, ig); 

fscanflfp, "Xlf, idcore); 

fscanflfp, “Xlf, ids); 

fscanflfp, “Xlf", idr); 

fscanflfp, "Xlf, ils); 

fscanflfp, "Xlf", Idr); 

fscanflfp, "Xlf", Ivol); 

fscanflfp, "Xlf", Iwt); 

fscanflfp, "Xlf, iph); 

fscanflfp, "Xlf", ipe); 

fscanflfp, "Xlf, ii2r); 

fscanflfp, "Xlf", iva); 

fscanflfp, "Xlf, ieffcy); 

fscanflfp, "Xlf", Lew); 

fscanflfp, "Xlf", il); 

fscanflfp, "Xlf", ii2rr); 

fscanflfp, "Xlf", isaax); 
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fscanflfp, “Xlftctaax); 
fscanfIfp, "Xlf% 4vt); 
fscanfIfp, "Xlf", 4slip); 

^scant(fp f "Xlf", &vag); 
fscanf Ifp, "Xlf", 4rl); 
fscanfIfp, “Xlf", ixl); 
fscanfIfp, ‘Ilf’, 4x«); 
fscanfIfp, "Xlf", 4x2); 
fscanf Ifp, "Xlf, 4r2); 
fcloseIfp); 

while If != 1) 

{ 

printf(“\nNhat is the sustained speed aachine horseponer? "); 
scanf("Xlf”, 4dhp); 
dhp t= 746.0; 

printf1"What is the sustained speed aachine rpa? '); 
scanfI"XIf", iparpa); 

freq = parpslp/60.0; It PH frequency 1/ 

via = (vtlxe)/sqrt((xl + xalllxl + xa) + Irllrll); 

It thevinin equivalent voltage tl 

rel = lxelxatrl)/(rUrl + (xl+xa)tlxl+xa)l; 

/I thev equiv resistance */ 

xel = (xattrltrl + xllxl + xltxad/lrltrl + lxl+xa)t(xl+xa)); 

/I thev equiv inductance tl 

sa - reltrel + Ixel+x2)t(xel+x2); It pieces of slip quadratic tl 
sb = 2trel*r2 - (3lllvla*vlalr2/dhp); 
sc = r2tr2; 

radical = sblsb - 4<satsc; 
if (radical < 0.0) 

abort I"\n6ot a negative radical in the slip eqn.‘); 
slipl - l-sb + sqrt(radical))/I2*sa); 
slip2 - l-sb - sqrtlradical))/l2*sa); 
if llslipl < 0.0) 44 Islip2 < 0.0)) ’ 

abort("\nTHO negative slips."); 

It noH we Hill use the saallest positive slip II 
else if llslipl < 0.0) 44 (slip2 > 0.0)) 
slip - slip2; 

else if llslipl > 0.0) 44 (slip2 < 0.01) 
slip = slipl; 

else 

slip = (slipl > slip2l ? slip2 : slipl; 

It tslip-ain(slipl, slip2) (/ 

if (slip > saax) 

abort I *\n slip is aore than saax."); 

siv = find_siv11, r, q, ds, dcore, Is); It stator iron voluae tl 

riv = llF'llrllr - 2ldrllr); It rotor iron voluae tl 

ntjron = Dllriv + siv); It iron Height only 1/ 

ph = 31.86225tBETA*freqtBRltHCllHt_iron/D; 

It hysteresis loss in natts, uses iron Height of aachine tl 
pe = 1106236.9*NU*BSATIBSAT*Tl*Tl*freq»freq*Ht.iron)/(RH0*D); 

It eddy current loss in natts, uses iron Height of aachine tl 


192 


> 


i2 = sqrt((siipldhp)/(3tr2l(l-slip))); 
il = i2 + vag/'xa; 


It load current tl 
It i2 plus lagnetizing current tl 


i2r=3.01iltiltrl; /I stator copper loss in watts tl 

i2rr = sliptdhp; It rotor copper loss in watts tl 

pieff = dhp/(dhp + ph + pe + i2r + i 2rr); 

printf(“\n Sustained efficiency is XIf“, pieff); 

printf(“\n\nWhat is the endurance speed lachine horsepower? *); 
scanfCW, &dhp); 
dhp 1= 746.0; 

printf("What is the endurance speed sachine rpi? B ); 
scanf( B Zlf B , ipirpa); 

freq = parpatp/60.0; It PH frequency tl 

sb = 2trellr2 - (3lltvlatvlatr2/dhp); 
radical = sbtsb - 4tsatsc; 
if (radical < 0.0) 

abort(“\n6ot a negative radical in the slip eqn. B ); 
slipi = (-sb + sqrt(radical))/(2tsa); 
slip2 = (-sb - sqrt(radical))/(2tsa); 
if ((slipi < 0.0) tt (slip2 < 0.0)) ' 

abort("\nTwo negative slips. 8 ); 
else if ((slipi < 0.0) M (slip2 > 0.0)) 
slip = slip2; 

else if ((slipi > 0.0) U (siip2 < 0.0)) 
slip = slipi; 

else 

slip = (slipi > slip2) ? slip2 : slipi; 
if (slip > siax) 

abort( B \n slip is lore than siax. B ); 
ph = 31.86225tBETAIfreqlBRltHCitwt_iron/B; 
pe = (106236.91NUJBSATlBSATtTllTllfreqlfreqtwtJron)/(RHGtD); 
i2 = sqrt((sliptdhp)/(3tr2f(1-slip))); It load current tl 

i1 = i2 + vag/xt; It i 2 plus magnetizing current tl 

i2r=3.0tiHi 1lr1; It stator copper loss in watts tl 

i2rr = sliptdhp; It rotor copper loss in watts tl 

pieff = dhp/(dhp + ph + pe + i2r + i2rr); 
printf( a \n Endurance efficiency is XIf B , pieff); 

printf( B \nSaae aachine? B ); 
scanf( B Xd B , fcf); 
if (f == 0) 

continue; 
else if (f == 2) 

{ 

e = 1; 
break; 

} 

} It end of f-loop tl 

} It end of e-loop tl 

} It end of tain prograi 1/ 
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double find_siv(l, r, g, ds, dcore, Is) 































































































double 1 , r, 9 , ds, dcore, Is; 

{ 

double one, two, three, four; 

one = <r+q+ds+dcore)t(r+g+ds+dcore) - (r+g)l(r+g); 
tno = 2IFII(r+g)tdstIs; 

three = (r+q+ds+dcore)I(r+q+ds+dcore) - (r+g+ds)l(r+g+ds); 

four = IMPUone - two) + PlHJ(r+g)tthree; 
return(four )5 
} 
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Appendix E. Weight and volume algorithms 


The ability to correctly characterize a new naval ship 
technology depends in part on the ability to calculate the 
weight and volume associated with that technology. Algo¬ 
rithms for this purpose were taken from a variety of 
sources, including the Advanced Surface Ship Evaluation Tool 
(ASSET) theory manuals. The ASSET algorithms are the result 
of data analysis for naval ships that have been constructed, 
as well as studies for other ship designs. 

The Ship Work Breakdown Structure (SWBS) categorizes 
all ship weights. The general categories are: 


W100 

Ship structures 

W200 

Propulsion plant 

W300 

Electric generation plant 

W400 

Command and control equipment 

W500 

Auxiliaries and distributed 
systems 

W600 

Outfit and furnishings 

W700 

Ship armament 

WFOO 

Ship Loads 


Within W200 are several sub-groups that pertain to 
electric propulsion. They are: 


W235 
W235.1 
W235.2 
W235.3 

W235.4 
W235.5 
W241.1 

W242 

W243 

W244 

W245 


Electric propulsion devices 
Propulsion motors 
Propulsion generators 
Transmission lines and 
propulsion cables 
Cooling systems 
Switchgear 

Locked-train-double-reduction 
reduction gears 
Propulsion clutches and 
couplings 
Shafting 

Propulsion shaft bearings 
Propeller weight 
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ASSET allows SWBS groups to be adjusted in weight, 
which allows technology sensitivity analyses, such as this 
thesis, to be performed. Only a few of the W235 sub-groups 
need to be calculated outside of ASSET and adjusted within 
ASSET. These are W235.1, W235.2, W235.4, W235.5, W241.1, 
and W243. The W235.1 and W235.2 weights and volumes were 
calculated as part of the machine design. The rest of the 
needed W235 weights were calculated in "wt.c", a copy of 
which follows. These algorithms are the U. S. Navy 
standard, and were verified against actual ships and com¬ 
ponents . 

Shafting and transmission line weights are dependent on 
motor and generator positions within a ship. Their weights 
and volumes were calculated from ASSET equations, using the 
layouts of the baseline and variant ships. 

Where no volume equation was found in the ASSET 
documentation, or where the result of such an equation was 
unrealistic, a literature search generally found enough ac¬ 
tual equipment to permit a relationship to be empirically 
determined. A linear scaling of such volumes provided 
adequate results. 
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Table 44. Listing of off-line weight and volume program 


tindude "stdio.h" 

tdefine K 150.75 It gear hardness factor 1/ 

It prograa naae: wt.c, to find the weights of SHBS groups tl 


aainO 

{ 

FILE Ifp, tfopen(); 

double xg, xa, zg, z», pa, ng, na, gra, grg, np, dhp, gn, ds, q, ne, 
w235, w2353, w2354, w2355, w24l=0.0, w243, w298=0.0, xprop, 
wpc, vpc, ppc, wexc, vexc, wbrk, vbrk, pa, v235, 
v2353=0.0, v2354=0.0, v2355, v241, v243; ’ 
extern double sqrtl), powl); 


printf(“\nThis prograa calculates 5HB5 200 weights for screen output."); 
printf(’\n\nHhat is the LCG of the propulsion generator(s)? "); 
scanfCXlf", &xg); 

printf("What is the LCG of the propulsion aotor(s)? '); 
scanfCXlf, ixai; 

printf("What is the VC6 of the propulsion generator<s)? "); 
scanfCXlf", izg); 

printf(“What is the VC6 of the propulsion aotorls)? "); 
scanfCXlf", 4za); 

printf("Hhat is the nuaber of propulsion generator Is)? "); 
scanfCXlf", ing); 

printf("Hhat is the nuaber of propulsion aotorls)? "); 
scanfCXlf", 4na); 

printf("Hhat is the nuaber of gas turbines aboard? "); 
scanfCXlf", (ne); 

printfI"What is the rated horsepower of each gas turbine? ’); 
scanfCXlf", 4pa); 

printf("What is the propeller rpa? '); 
scanfCXlf", 4np); 

printf("How auch horsepower is delivered to each propeller? "); 
scanfCXlf", Mhp); 

printfI“What is the LCG of the propeller? "); 
scanfCXlf", ixprop); 

printfI"Hhat is the gear ratio at the propulsion aotor? "); 
scanfCXlf", igra); 

printf("What is the gear ratio at the propulsion generator? '); 
scanfCXlf", tgrg); 


It H235.3 Transaission lines tl 

w2353 = 0.0000091(xa-xg+zg-za+27)I(pal746/30000>; It LT, enhanced tl 
v2353 = 0.06545*(xa-xg+zg-za+27); It cubic feet, enhanced tl 


It swag tl 


It N235.4 Cooling systeas tl 
w2354 = 0.26tpalng/2240; 
v2354 = lOO.OIng; 


197 


It LT 1/ 



/• H235.5 Switch gear 1/ 

*2355 = 0.26t(ng + 2lni>; It enhanced ASSET, LT 1/ 

v2355 = 45.0<w2355; It switchgear volume, ft A 3, enhanced I/ 


ppc r 0.000939523lpa/m; 
wpc - ni(3.55tppc/12; 
vpc = nit540lppc/12; 


It power converter rating, HH, enhanced ASSET I/ 
It weight pwr conv, LT, enhanced ASSET 1/ 
It vol pwr conv, LT, enhanced ASSET 1/ 


pi = pa!746/1000000; It *otor rating in HN 1/ 

wexc = (ne * ni)lpow((pi/30), 0.3); It weight of exciters, LT, enhanced I/ 

vexc = 79twexc; It vol of exciters, ft A 3 tl 

wbrk = 0.26lpi<m; It weight of braking resistors, LT, enhanced t/ 

vbrk = 37lpilm; It vol of resistors, ft A 3, enhanced 1/ 


w235 = w2353 + w2354 + w2355 + wpc + wexc + wbrk; 
v235 = v2353 + v2354 + v2355 + vpc + vexc + vbrk; 

It H241 locked-train double reduction gears 1/ 

if (grg != 1.0) It there are pg gears, btw gt and pg 1/ 

w241 = (l.57lpatpow((grg+l), 3.0)lng)/(3600lgrglgrglgrg4K); 
if (gri != 1.0) It there are pi gears, btw ps and propeller tl 

w241 += (1.57tdhptpow((gri+l), 3.0Uni)/(nplgrilgriiK); It LT 1/ 

v24l = w24H34.612; It cubic feet, ratioed fi FFG7 tl 

It N243 Shafting tl 

ds = 2.152Ipow((4.22347tdhp/np>, 0.333); It shaft dial, sq inches tl 
w243 = 1.5708lds*dsl2l(xprop - xi - 1.01/2240; It LT 1/ 

v243 = 0.01091ldsldsl(xprop - xi - 6.0); It ft A 3 1/ 


It W298.1 LTDR operating fluids, additon to wt in ASSET tl 
if (grg != 1.0) 

w298 = 0.27lpalng/2240; 
if (gri != 1.0) 

w298 += 0.27lpalm/2240; It LT 1/ 


printf C\n\nH235.3 Transiission lines 
printf(‘\nN235.4 Cooling systeis 
printff"\nW235.5 Switchgear 
printff*\n Power converters 
printf( , \n Exciters 
printf("\n Braking resistors 
printf("\n\nH235 Electric propulsion 
printf C\nH241 Reduction gears 
printf C\nH243 Shafting 
printf C\nH298 6ear operating fluid 
} 


I7.21f LT 17.21 f ft A 3*, *2353, v2353); 

Z7.21f LT X7.2H ft A 3’, w2354, v2354); 

17.21f LT Z7.21f ft A 3", w2355, v2355); 

17.21f LT 17.21f ft A 3‘, wpc, vpc); 

17.21f LT 17.21f ft A 3‘, wexc, vexc); 

17.21 f LT 17.2H ft A 3‘, wbrk, vbrk); 

17.21f LT 17.21f ft A 3, less PGs and PMs 1 , w235, v235) 
17.21f LT 17.21f ft A 3‘, w241, v241>; 

17.21f LT 17.21f ft A 3", w243, v243); 

17.21f LT\n\n‘, w298); 
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Appendix F. Advanced Surface Ship Evaluation Tool output 

The output of ASSET is in text and graphic form. The 
total text output for any particular synthesis run is more 
than thirty pages. Following are several graphic outputs of 
ASSET, showing the mechanical and electrical transmission 
ships used in this thesis. 
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Figure 49. Hull isometric view of all ships 
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Figure 50. Body plan of all ships 
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Figure 51. Plan view of subdivision in mechanical baseline 
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Figure 52. Plan view of subdivision in rearranged electrical 
ship 
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